Induction of IL-10-producing type 2 innate lymphoid cells by allergen immunotherapy is associated with clinical response.

The role of innate immune cells in allergen immunotherapy that confers immune tolerance to the sensitizing allergen is unclear. Here, we report a role of interleukin-10-producing type 2 innate lymphoid cells (IL-10+ ILC2s) in modulating grass-pollen allergy. We demonstrate that KLRG1+ but not KLRG1- ILC2 produced IL-10 upon activation with IL-33 and retinoic acid. These cells attenuated Th responses and maintained epithelial cell integrity. IL-10+ KLRG1+ ILC2s were lower in patients with grass-pollen allergy when compared to healthy subjects. In a prospective, double-blind, placebo-controlled trial, we demonstrated that the competence of ILC2 to produce IL-10 was restored in patients who received grass-pollen sublingual immunotherapy. The underpinning mechanisms were associated with the modification of retinol metabolic pathway, cytokine-cytokine receptor interaction, and JAK-STAT signaling pathways in the ILCs. Altogether, our findings underscore the contribution of IL-10+ ILC2s in the disease-modifying effect by allergen immunotherapy.

Developmental and biophysical determinants of grass leaf size worldwide.

One of the most notable ecological trends-described more than 2,300  years ago by Theophrastus-is the association of small leaves with dry and cold climates, which has recently been recognized for eudicotyledonous plants at a global scale1-3. For eudicotyledons, this pattern has been attributed to the fact that small leaves have a thinner boundary layer that helps to avoid extreme leaf temperatures4 and their leaf development results in vein traits that improve water transport under cold or dry climates5,6. However, the global distribution of leaf size and its adaptive basis have not been tested in the grasses, which represent a diverse lineage that is distinct in leaf morphology and that contributes 33% of terrestrial primary productivity (including the bulk of crop production)7. Here we demonstrate that grasses have shorter and narrower leaves under colder and drier climates worldwide. We show that small grass leaves have thermal advantages and vein development that contrast with those of eudicotyledons, but that also explain the abundance of small leaves in cold and dry climates. The worldwide distribution of leaf size in grasses exemplifies how biophysical and developmental processes result in convergence across major lineages in adaptation to climate globally, and highlights the importance of leaf size and venation architecture for grass performance in past, present and future ecosystems.

Investigating the cis-regulatory basis of C3 and C4 photosynthesis in grasses at single-cell resolution.

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-dependent malic enzyme), Panicum miliaceum (NAD-dependent malic enzyme), Urochloa fusca (phosphoenolpyruvate carboxykinase), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of cis-regulatory elements evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Genome evolution and initial breeding of the Triticeae grass Leymus chinensis dominating the Eurasian Steppe.

Leymus chinensis, a dominant perennial grass in the Eurasian Steppe, is well known for its remarkable adaptability and forage quality. Hardly any breeding has been done on the grass, limiting its potential in ecological restoration and forage productivity. To enable genetic improvement of the untapped, important species, we obtained a 7.85-Gb high-quality genome of L. chinensis with a particularly long contig N50 (318.49 Mb). Its allotetraploid genome is estimated to originate 5.29 million years ago (MYA) from a cross between the Ns-subgenome relating to Psathyrostachys and the unknown Xm-subgenome. Multiple bursts of transposons during 0.433-1.842 MYA after genome allopolyploidization, which involved predominantly the Tekay and Angela of LTR retrotransposons, contributed to its genome expansion and complexity. With the genome resource available, we successfully developed a genetic transformation system as well as the gene-editing pipeline in L. chinensis. We knocked out the monocot-specific miR528 using CRISPR/Cas9, resulting in the improvement of yield-related traits with increases in the tiller number and growth rate. Our research provides valuable genomic resources for Triticeae evolutionary studies and presents a conceptual framework illustrating the utilization of genomic information and genome editing to accelerate the improvement of wild L. chinensis with features such as polyploidization and self-incompatibility.

Paleolakes and socioecological implications of last glacial "greening" of the South African interior.

Determining the timing and drivers of Pleistocene hydrological change in the interior of South Africa is critical for testing hypotheses regarding the presence, dynamics, and resilience of human populations. Combining geological data and physically based distributed hydrological modeling, we demonstrate the presence of large paleolakes in South Africa's central interior during the last glacial period, and infer a regional-scale invigoration of hydrological networks, particularly during marine isotope stages 3 and 2, most notably 55 to 39 ka and 34 to 31 ka. The resulting hydrological reconstructions further permit investigation of regional floral and fauna responses using a modern analog approach. These suggest that the climate change required to sustain these water bodies would have replaced xeric shrubland with more productive, eutrophic grassland or higher grass-cover vegetation, capable of supporting a substantial increase in ungulate diversity and biomass. The existence of such resource-rich landscapes for protracted phases within the last glacial period likely exerted a recurrent draw on human societies, evidenced by extensive pan-side artifact assemblages. Thus, rather than representing a perennially uninhabited hinterland, the central interior's underrepresentation in late Pleistocene archeological narratives likely reflects taphonomic biases stemming from a dearth of rockshelters and regional geomorphic controls. These findings suggest that South Africa's central interior experienced greater climatic, ecological, and cultural dynamism than previously appreciated and potential to host human populations whose archaeological signatures deserve systematic investigation.

Mutations in NAKED-ENDOSPERM IDD genes reveal functional interactions with SCARECROW during leaf patterning in C4 grasses.

Leaves comprise a number of different cell-types that are patterned in the context of either the epidermal or inner cell layers. In grass leaves, two distinct anatomies develop in the inner leaf tissues depending on whether the leaf carries out C3 or C4 photosynthesis. In both cases a series of parallel veins develops that extends from the leaf base to the tip but in ancestral C3 species veins are separated by a greater number of intervening mesophyll cells than in derived C4 species. We have previously demonstrated that the GRAS transcription factor SCARECROW (SCR) regulates the number of photosynthetic mesophyll cells that form between veins in the leaves of the C4 species maize, whereas it regulates the formation of stomata in the epidermal leaf layer in the C3 species rice. Here we show that SCR is required for inner leaf patterning in the C4 species Setaria viridis but in this species the presumed ancestral stomatal patterning role is also retained. Through a comparative mutant analysis between maize, setaria and rice we further demonstrate that loss of NAKED-ENDOSPERM (NKD) INDETERMINATE DOMAIN (IDD) protein function exacerbates loss of function scr phenotypes in the inner leaf tissues of maize and setaria but not rice. Specifically, in both setaria and maize, scr;nkd mutants exhibit an increased proportion of fused veins with no intervening mesophyll cells. Thus, combined action of SCR and NKD may control how many mesophyll cells are specified between veins in the leaves of C4 but not C3 grasses. Together our results provide insight into the evolution of cell patterning in grass leaves and demonstrate a novel patterning role for IDD genes in C4 leaves.

Standing genetic variation fuels rapid evolution of herbicide resistance in blackgrass.

Repeated herbicide applications in agricultural fields exert strong selection on weeds such as blackgrass (Alopecurus myosuroides), which is a major threat for temperate climate cereal crops. This inadvertent selection pressure provides an opportunity for investigating the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much target-site resistance (TSR) relies on de novo mutations vs. standing variation is important for developing strategies to manage herbicide resistance. We first generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Next, through empirical data in the form of highly accurate long-read amplicons of alleles encoding acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) variants, we showed that most populations with resistance due to TSR mutations-23 out of 27 and six out of nine populations for ACCase and ALS, respectively-contained at least two TSR haplotypes, indicating that soft sweeps are the norm. Finally, through forward-in-time simulations, we inferred that TSR is likely to mainly result from standing genetic variation, with only a minor role for de novo mutations.

Cross-species predictive modeling reveals conserved drought responses between maize and sorghum.

Drought tolerance is a highly complex trait controlled by numerous interconnected pathways with substantial variation within and across plant species. This complexity makes it difficult to distill individual genetic loci underlying tolerance, and to identify core or conserved drought-responsive pathways. Here, we collected drought physiology and gene expression datasets across diverse genotypes of the C4 cereals sorghum and maize and searched for signatures defining water-deficit responses. Differential gene expression identified few overlapping drought-associated genes across sorghum genotypes, but using a predictive modeling approach, we found a shared core drought response across development, genotype, and stress severity. Our model had similar robustness when applied to datasets in maize, reflecting a conserved drought response between sorghum and maize. The top predictors are enriched in functions associated with various abiotic stress-responsive pathways as well as core cellular functions. These conserved drought response genes were less likely to contain deleterious mutations than other gene sets, suggesting that core drought-responsive genes are under evolutionary and functional constraints. Our findings support a broad evolutionary conservation of drought responses in C4 grasses regardless of innate stress tolerance, which could have important implications for developing climate resilient cereals.

Similar photosynthetic but different yield responses of C3 and C4 crops to elevated O3.

The deleterious effects of ozone (O3) pollution on crop physiology, yield, and productivity are widely acknowledged. It has also been assumed that C4 crops with a carbon concentrating mechanism and greater water use efficiency are less sensitive to O3 pollution than C3 crops. This assumption has not been widely tested. Therefore, we compiled 46 journal articles and unpublished datasets that reported leaf photosynthetic and biochemical traits, plant biomass, and yield in five C3 crops (chickpea, rice, snap bean, soybean, and wheat) and four C4 crops (sorghum, maize, Miscanthus × giganteus, and switchgrass) grown under ambient and elevated O3 concentration ([O3]) in the field at free-air O3 concentration enrichment (O3-FACE) facilities over the past 20 y. When normalized by O3 exposure, C3 and C4 crops showed a similar response of leaf photosynthesis, but the reduction in chlorophyll content, fluorescence, and yield was greater in C3 crops compared with C4 crops. Additionally, inbred and hybrid lines of rice and maize showed different sensitivities to O3 exposure. This study quantitatively demonstrates that C4 crops respond less to elevated [O3] than C3 crops. This understanding could help maintain cropland productivity in an increasingly polluted atmosphere.

Reinventing metabolic pathways: Independent evolution of benzoxazinoids in flowering plants.

Benzoxazinoids (BXDs) form a class of indole-derived specialized plant metabolites with broad antimicrobial and antifeedant properties. Unlike most specialized metabolites, which are typically lineage-specific, BXDs occur sporadically in a number of distantly related plant orders. This observation suggests that BXD biosynthesis arose independently numerous times in the plant kingdom. However, although decades of research in the grasses have led to the elucidation of the BXD pathway in the monocots, the biosynthesis of BXDs in eudicots is unknown. Here, we used a metabolomic and transcriptomic-guided approach, in combination with pathway reconstitution in Nicotiana benthamiana, to identify and characterize the BXD biosynthetic pathways from both Aphelandra squarrosa and Lamium galeobdolon, two phylogenetically distant eudicot species. We show that BXD biosynthesis in A. squarrosa and L. galeobdolon utilize a dual-function flavin-containing monooxygenase in place of two distinct cytochrome P450s, as is the case in the grasses. In addition, we identified evolutionarily unrelated cytochrome P450s, a 2-oxoglutarate-dependent dioxygenase, a UDP-glucosyltransferase, and a methyltransferase that were also recruited into these BXD biosynthetic pathways. Our findings constitute the discovery of BXD pathways in eudicots. Moreover, the biosynthetic enzymes of these pathways clearly demonstrate that BXDs independently arose in the plant kingdom at least three times. The heterogeneous pool of identified BXD enzymes represents a remarkable example of metabolic plasticity, in which BXDs are synthesized according to a similar chemical logic, but with an entirely different set of metabolic enzymes.

Wild grass-derived alleles represent a genetic architecture for the resilience of modern common wheat to stresses.

This review explores the integration of wild grass-derived alleles into modern bread wheat breeding to tackle the challenges of climate change and increasing food demand. With a focus on synthetic hexaploid wheat, this review highlights the potential of genetic variability in wheat wild relatives, particularly Aegilops tauschii, for improving resilience to multifactorial stresses like drought, heat, and salinity. The evolutionary journey of wheat (Triticum spp.) from diploid to hexaploid species is examined, revealing significant genetic contributions from wild grasses. We also emphasize the importance of understanding incomplete lineage sorting in the genomic evolution of wheat. Grasping this information is crucial as it can guide breeders in selecting the appropriate alleles from the gene pool of wild relatives to incorporate into modern wheat varieties. This approach improves the precision of phylogenetic relationships and increases the overall effectiveness of breeding strategies. This review also addresses the challenges in utilizing the wheat wild genetic resources, such as the linkage drag and cross-compatibility issues. Finally, we culminate the review with future perspectives, advocating for a combined approach of high-throughput phenotyping tools and advanced genomic techniques to comprehensively understand the genetic and regulatory architectures of wheat under stress conditions, paving the way for more precise and efficient breeding strategies.

Short day promotes gall swelling by a CONSTANS-FLOWERING LOCUS T pathway in Zizania latifolia.

"Jiaobai" is a symbiont of Zizania latifolia and Ustilago esculenta, producing fleshy galls as a popular vegetable in South and East Asia. Current "Jiaobai" cultivars exhibit abundant variation in their gall formation date; however, the underlying mechanism is not clear. In this study, a strict short-day (SD) "Jiaobai" line "YD-3" was used. Plants were treated with two day-length regimes [14 h/10 h (day/night) (control) and 8 h/16 h (day/night) (SD)] from 100 to 130 days after planting. The gall swelling rate of the two treatments and another early SD treatment (from 60 to 90 days after planting), together with the contingent flowering plants in the experiment population, revealed that SD can improve both gall enlargement and flowering of "Jiaobai" plants. Comparison of RNA sequencing data among control, SD swelling, and SD flowering treatments of leaves and meristems indicated that SD promotion of "Jiaobai" swelling is conducted by the CONSTANS (CO)-FLOWERING LOCUS T (FT) pathway, similar but not identical to the SD-induced flowering pathway in Z latifolia and rice. "Virus-induced gene silencing", "Yeast one-hybrid assay" and "Dual-luciferase assay" showed that a FT gene, ZlGsd1, is critical in SD promotion of gall formation and is positively regulated by a CO gene, ZlCOL1. Our study elucidated how photoperiod affects the formation of a unique organ produced by plant-fungus symbiosis. The difference in SD response between "Jiaobai" and rice, as well as their potential applications in breeding of "Jiaobai" and rice, were also discussed.

Xylan-directed cell wall assembly in grasses.

Xylan is the most abundant hemicellulosic polysaccharide in the cell walls of grasses and is pivotal for the assembly of distinct cell wall structures that govern various cellular functions. Xylan also plays a crucial role in regulating biomass recalcitrance, ultimately affecting the utilization potential of lignocellulosic materials. Over the past decades, our understanding of the xylan biosynthetic machinery and cell wall organization has substantially improved due to the innovative application of multiple state-of-the-art techniques. Notably, novel xylan-based nanostructures have been revealed in the cell walls of xylem vessels, promoting a more extensive exploration of the role of xylan in the formation of cell wall structures. This Update summarizes recent achievements in understanding xylan biosynthesis, modification, modeling, and compartmentalization in grasses, providing a brief overview of cell wall assembly regarding xylan. We also discuss the potential for tailoring xylan to facilitate the breeding of elite energy and feed crops.

Chromosome-specific painting reveals the Y genome origin and chromosome rearrangements of the St genome in Triticeae.

Karyotypes provide key cytogenetic information on phylogenetic relationships and evolutionary origins in related plant species. The St genome of Pseudoroegneria contributes to 8 alloploid genera, representing over half of the species that are highly valuable for wheat (Triticum aestivum) breeding and for understanding Triticeae species evolution. However, St chromosome characterization is challenging due to limited cytogenetic markers and DNA information. We developed a complete set of St genome-specific chromosome painting probes for identification of the individual chromosomes 1St to 7St based on the genome sequences of Pseudoroegneria libanotica and wheat. We revealed the conservation of St chromosomes in St-containing species by chromosome painting, including Pseudoroegneria, Roegneria, Elymus, and Campeiostachys. Notably, the Y genome showed hybridization signals, albeit weaker than those of the St genome. The awnless species harboring the Y genome exhibited more intense hybridization signals compare to the awned species in Roegneria and Campeiostachys, yet weaker than the hybridization signals of the St genome in autotetraploid Pseudoroegneria strigosa. Although awnless species were morphologically more similar to each other, phenotypic divergence progressively increased from awnless to awned species. Our results indicate that the Y genome originated from the St genome and shed light on the possible origin of the Roegneria and Campeiostachys species, enhancing our understanding of St-genome-containing species evolution.

Poaceae plants transfer cyclobutane pyrimidine dimer photolyase to chloroplasts for ultraviolet-B resistance.

Photoreactivation enzyme that repairs cyclobutane pyrimidine dimer (CPD) induced by ultraviolet-B radiation, commonly called CPD photolyase (PHR) is essential for plants living under sunlight. Rice (Oryza sativa) PHR (OsPHR) is a unique triple-targeting protein. The signal sequences required for its translocation to the nucleus or mitochondria are located in the C-terminal region but have yet to be identified for chloroplasts. Here, we identified sequences located in the N-terminal region, including the serine-phosphorylation site at position 7 of OsPHR, and found that OsPHR is transported/localized to chloroplasts via a vesicle transport system under the control of serine-phosphorylation. However, the sequence identified in this study is only conserved in some Poaceae species, and in many other plants, PHR is not localized to the chloroplasts. Therefore, we reasoned that Poaceae species need the ability to repair CPD in the chloroplast genome to survive under sunlight and have uniquely acquired this mechanism for PHR chloroplast translocation.

Variation in leaf dark respiration among C3 and C4 grasses is associated with use of different substrates.

Measurements of respiratory properties have often been made at a single time point either during daytime using dark-adapted leaves or during nighttime. The influence of the day-night cycle on respiratory metabolism has received less attention but is crucial to understand photosynthesis and photorespiration. Here, we examined how CO2- and O2-based rates of leaf dark respiration (Rdark) differed between midday (after 30-min dark adaptation) and midnight in 8 C3 and C4 grasses. We used these data to calculate the respiratory quotient (RQ; ratio of CO2 release to O2 uptake), and assessed relationships between Rdark and leaf metabolome. Rdark was higher at midday than midnight, especially in C4 species. The day-night difference in Rdark was more evident when expressed on a CO2 than O2 basis, with the RQ being higher at midday than midnight in all species, except in rice (Oryza sativa). Metabolomic analyses showed little correlation of Rdark or RQ with leaf carbohydrates (sucrose, glucose, fructose, or starch) but strong multivariate relationships with other metabolites. The results suggest that rates of Rdark and differences in RQ were determined by several concurrent CO2-producing and O2-consuming metabolic pathways, not only the tricarboxylic acid cycle (organic acids utilization) but also the pentose phosphate pathway, galactose metabolism, and secondary metabolism. As such, Rdark was time-, type- (C3/C4) and species-dependent, due to the use of different substrates.

Rapid growth of Moso bamboo (Phyllostachys edulis): Cellular roadmaps, transcriptome dynamics, and environmental factors.

Moso bamboo (Phyllostachys edulis) shows remarkably rapid growth (114.5 cm/day), but the underlying biological mechanisms remain unclear. After examining more than 12,750 internodes from more than 510 culms from 17 Moso populations, we identified internode 18 as a representative internode for rapid growth. This internode includes a 2-cm cell division zone (DZ), a cell elongation zone up to 12 cm, and a secondary cell wall (SCW) thickening zone. These zones elongated 11.8 cm, produced approximately 570,000,000 cells, and deposited ∼28 mg g-1 dry weight (DW) lignin and ∼44 mg g-1 DW cellulose daily, far exceeding vegetative growth observed in other plants. We used anatomical, mathematical, physiological, and genomic data to characterize development and transcriptional networks during rapid growth in internode 18. Our results suggest that (1) gibberellin may directly trigger the rapid growth of Moso shoots, (2) decreased cytokinin and increased auxin accumulation may trigger cell DZ elongation, and (3) abscisic acid and mechanical pressure may stimulate rapid SCW thickening via MYB83L. We conclude that internode length involves a possible tradeoff mediated by mechanical pressure caused by rapid growth, possibly influenced by environmental temperature and regulated by genes related to cell division and elongation. Our results provide insight into the rapid growth of Moso bamboo.

Genetic control of branching patterns in grass inflorescences.

Inflorescence branching in the grasses controls the number of florets and hence the number of seeds. Recent data on the underlying genetics come primarily from rice and maize, although new data are accumulating in other systems as well. This review focuses on a window in developmental time from the production of primary branches by the inflorescence meristem through to the production of glumes, which indicate the transition to producing a spikelet. Several major developmental regulatory modules appear to be conserved among most or all grasses. Placement and development of primary branches are controlled by conserved auxin regulatory genes. Subtending bracts are repressed by a network including TASSELSHEATH4, and axillary branch meristems are regulated largely by signaling centers that are adjacent to but not within the meristems themselves. Gradients of SQUAMOSA-PROMOTER BINDING-like and APETALA2-like proteins and their microRNA regulators extend along the inflorescence axis and the branches, governing the transition from production of branches to production of spikelets. The relative speed of this transition determines the extent of secondary and higher order branching. This inflorescence regulatory network is modified within individual species, particularly as regards formation of secondary branches. Differences between species are caused both by modifications of gene expression and regulators and by presence or absence of critical genes. The unified networks described here may provide tools for investigating orphan crops and grasses other than the well-studied maize and rice.

The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata.

The unique morphology of grass stomata enables rapid responses to environmental changes. Deciphering the basis for these responses is critical for improving food security. We have developed a planta platform of single-nucleus RNA-sequencing by combined fluorescence-activated nuclei flow sorting, and used it to identify cell types in mature and developing stomata from 33,098 nuclei of the maize epidermis-enriched tissues. Guard cells (GCs) and subsidiary cells (SCs) displayed differential expression of genes, besides those encoding transporters, involved in the abscisic acid, CO2, Ca2+, starch metabolism, and blue light signaling pathways, implicating coordinated signal integration in speedy stomatal responses, and of genes affecting cell wall plasticity, implying a more sophisticated relationship between GCs and SCs in stomatal development and dumbbell-shaped guard cell formation. The trajectory of stomatal development identified in young tissues, and by comparison to the bulk RNA-seq data of the MUTE defective mutant in stomatal development, confirmed known features, and shed light on key participants in stomatal development. Our study provides a valuable, comprehensive, and fundamental foundation for further insights into grass stomatal function.

Alpine Extremophytes in Evolutionary Turmoil: Complex Diversification Patterns and Demographic Responses of a Halophilic Grass in a Central Asian Biodiversity Hotspot.

Diversification and demographic responses are key processes shaping species evolutionary history. Yet we still lack a full understanding of ecological mechanisms that shape genetic diversity at different spatial scales upon rapid environmental changes. In this study, we examined genetic differentiation in an extremophilic grass Puccinellia pamirica and factors affecting its population dynamics among the occupied hypersaline alpine wetlands on the arid Pamir Plateau in Central Asia. Using genomic data, we found evidence of fine-scale population structure and gene flow among the localities established across the high-elevation plateau as well as fingerprints of historical demographic expansion. We showed that an increase in the effective population size could coincide with the Last Glacial Period, which was followed by the species demographic decline during the Holocene. Geographic distance plays a vital role in shaping the spatial genetic structure of P. pamirica alongside with isolation-by-environment and habitat fragmentation. Our results highlight a complex history of divergence and gene flow in this species-poor alpine region during the Late Quaternary. We demonstrate that regional climate specificity and a shortage of nonclimate data largely impede predictions of future range changes of the alpine extremophile using ecological niche modeling. This study emphasizes the importance of fine-scale environmental heterogeneity for population dynamics and species distribution shifts.