High light intensity plays a major role in emergence of population level variation in Arabidopsis thaliana along an altitudinal gradient

Genomic Insights into High-Altitude Adaptation: A Comparative Analysis of Roscoea alpina and R. purpurea in the Himalayas

Environmental stress at high altitudes drives the development of distinct adaptive mechanisms in plants. However, studies exploring the genetic adaptive mechanisms of high-altitude plant species are scarce. In the present study, we explored the high-altitude adaptive mechanisms of plants in the Himalayas through whole-genome resequencing. We studied two widespread members of the Himalayan endemic alpine genus Roscoea (Zingiberaceae): R. alpina (a selfing species) and R. purpurea (an outcrossing species). These species are distributed widely in the Himalayas with distinct non-overlapping altitude distributions; R. alpina is distributed at higher elevations, and R. purpurea occurs at lower elevations. Compared to R. purpurea, R. alpina exhibited higher levels of linkage disequilibrium, Tajima's D, and inbreeding coefficient, as well as lower recombination rates and genetic diversity. Approximately 96.3% of the genes in the reference genome underwent significant genetic divergence (FST ≥ 0.25). We reported 58 completely divergent genes (FST = 1), of which only 17 genes were annotated with specific functions. The functions of these genes were primarily related to adapting to the specific characteristics of high-altitude environments. Our findings provide novel insights into how evolutionary innovations promote the adaptation of mountain alpine species to high altitudes and harsh habitats.

Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na⁺ -K⁺ /H⁺ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

Polymorphisms in alternative oxidase genes from ecotypes of Arabidopsis and rice revealed an environment-induced linkage to altitude and rainfall

We investigated SNPs in alternative oxidase (AOX) genes and their connection to ecotype origins (climate, altitude, and rainfall) by using genomic data sets of Arabidopsis and rice populations from 1190 and 90 ecotypes, respectively. Parameters were defined to detect non-synonymous SNPs in the AOX ORF, which revealed amino acid (AA) changes in AOX1c, AOX1d, and AOX2 from Arabidopsis and AOX1c from rice in comparison to AOX references from Columbia-0 and Japonica ecotypes, respectively. Among these AA changes, Arabidopsis AOX1c_A161E&G165R and AOX1c_R242S revealed a link to high rainfall and high altitude, respectively, while all other changes in Arabidopsis and rice AOX was connected to high altitude and rainfall. Comparative 3D modeling showed that all mutant AOX presented structural differences in relation to the respective references. Molecular docking analysis uncovered lower binding affinity values between AOX and the substrate ubiquinol for most of the identified structures compared to their reference, indicating better enzyme-substrate binding affinities. Thus, our in silico data suggest that the majority of the AA changes found in the available ecotypes will confer better enzyme-subtract interactions and thus indicate environment-related, more efficient AOX activity.

Geometric morphometric analysis of the plateau zokor (Eospalax baileyi) revealed significant effects of environmental factors on skull variations

The plateau zokor (Eospalax baileyi) is employed as an ideal model for examining the relationships between phenotypic and ecological adaptations to the underground conditions in which the skull morphology evolves to adapt to tunnel environment. We evaluated the influence of environmental factors (altitude, temperature, and precipitation) and geographical distance on the variations in skull morphology of a native subterranean rodent plateau zokor population. Thin-plate spline showed that the trend of morphological changes along the CV1 axis was as follows: the two zygomatic arch and the two postorbital processes moved down, the two mastoid processes and the tooth row moved upward, and the tympanic bulla grew longer. The changes along the CV2 axis were as follows: the nasal bone and the tooth row became longer, the distance between the two anterior tips of zygomatic arch lengthened, the infraorbital foramen became smaller, the whole posterior part of the skull became shorter, the zygomatic bone and the two posterior tips of zygomatic arch moved down, and the foramen magnum became bigger. Thus we found significant differences in the skull shape among the seven populations studied. Along with the reduction in the altitude and increase in the mean annual temperature and mean annual precipitation, the nasal bone became shorter, the distance between the two anterior tips of the zygomatic arch became shorter, the whole posterior part of the skull lengthened, the infraorbital foramen became smaller, the two mastoid processes moved upward, and the occipital bone moved down on the dorsal surface of the skull. On the ventral surface of the skull, with an increase in the altitude, mean annual temperature, and mean annual precipitation, the tympanic bulla became shorter, the tooth row moved down, and the foramen magnum became smaller. The morphological changes in the skull were significantly positively correlated with environmental factors. Finally, there was a significant positive correlation between the Procrustes distance matrix of the skull and the geographic distance matrix, which indicates that the evolution of the plateau zokor follows the distance isolation model, but it needs to be further explored from genetic perspectives.

Modulation of miRNA expression in natural populations of A. thaliana along a wide altitudinal gradient of Indian Himalayas

Plant populations growing along an altitudinal gradient are exposed to different environmental conditions. They are excellent resources to study regulatory mechanisms adopted by plants to respond to different environmental stresses. Regulation by miRNA is one of such strategies. Here, we report how different miRNAs are preferentially expressed in the three natural populations of A. thaliana originating from a wide altitudinal range. The expression level of miRNAs was mostly governed by temperature and radiation. Majority of the identified miRNAs expressed commonly in the three populations. However, 30 miRNAs expressed significantly at different level between the low and the high altitude populations. Most of these miRNAs regulate the genes associated with different developmental processes, abiotic stresses including UV, cold, secondary metabolites, etc. Further, the expression of miR397 and miR858 involved in lignin biosynthesis and regulation of secondary metabolites respectively, may be regulated by light intensity. A few miRNAs expressed at increasing level with the increase in the altitude of the site indicating environment driven tight regulation of these miRNAs. Further, several novel miRNAs and isomiR diversity specific to the Himalayas are reported which might have an adaptive advantage. To the best of our knowledge, this is the first report on miRNA expression from natural plant populations.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome-wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome-wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome-wide associations identified known (Twin Sister of FT, FRIGIDA-like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage-Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.

High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition

BACKGROUND

Population differentiation and their adaptation to a particular environment depend on their ability to respond to a new environment. This, in turn is governed to an extent, by the degree of phenotypic plasticity exhibited by the populations. The populations of same species inhabiting different climatic conditions may differ in their phenotypic plasticity. Himalayan populations of Arabidopsis thaliana originating from a steep altitude are exposed to different climatic conditions ranging from sub-tropical to temperate. Thus they might have experienced different selection pressures during evolution and may respond differently under common environmental condition.

RESULTS

Phenotypic plasticity and differentiation of natural populations of A. thaliana grown under common garden and controlled conditions were determined. A total of seventeen morphological traits, their plasticity, association between traits and environment were performed using 45 accessions from three populations. Plants from different altitudes differed in phenotypes, their selection and fitness under two conditions. Under both the conditions lower altitude population was characterized by higher leaf count and larger silique than higher and middle altitude population. Flowering time of high altitude population increased while that of low and medium altitude decreased under controlled condition compared to open field. An increase in seed weight and germination was observed for all the population under open field than controlled. Rosette area was under divergent selection in both the condition. The heritability of lower altitude population was the highest under both the conditions, where as it was the least for higher altitude population further indicating that the high altitude populations are more responsive towards phenotypic changes under new environmental conditions. Ninety-nine percent of variability in traits and their plasticity co-varied with the altitude of their origin. The population of high altitude was more plastic and differentiated as compared to the lower altitude one.

CONCLUSIONS

Arabidopsis thaliana population native to different altitudes of the west Himalaya responds differently when grown under common environments. The success of high altitude population is more in common garden than the controlled conditions. The significant variability in phenotype and its association with altitude of origin predicts for non-random genetic differentiation among the populations.