Changes in gene expression patterns can lead to the variation of morphological traits. This phenomenon is particularly evident in recent evolution events such as crop domestication and responses to environmental stress, where alterations in expression levels can efficiently give rise to domesticated syndromes and adaptive phenotypes. Rice (Oryza sativa L.), one of the world’s most crucial cereal crops, comprises two morphologically distinct subspecies, Indica and Japonica. To investigate the morphological divergence between these two rice subspecies, this study planted a total of 315 landrace individuals of both Indica and Japonica under identical cultivation conditions. Out of the 16 quantitative traits measured in this study, 12 exhibited significant differences between the subspecies. To determine the genetic divergence between Indica and Japonica at the whole-genome sequence level, we constructed a phylogenetic tree using a resequencing dataset encompassing 95 rice landrace accessions. The samples formed two major groups that neatly corresponded to the two subspecies, Indica and Japonica. Furthermore, neighbor-joining (NJ) trees based on the expression quantity of effectively expressed genes (EEGs) across five different tissues categorized 12 representative samples into two major clades aligning with the two subspecies. These results imply that divergence in genome-wide expression levels undergoes stabilizing selection under non-stressful conditions, with evolutionary trends in expression levels mirroring sequence variation levels. This study further supports the pivotal role of changes in genome-wide expression regulation in the divergence of the two rice subspecies, Indica and Japonica.

Kingdonia uniflora is an endangered alpine herb that is distributed along an altitudinal gradient. The unique traits and important phylogenetic position make K. uniflora an ideal model for exploring how endangered plants react to altitude variation. In this study, we sampled nine individuals from three representative locations and adopted RNA-seq technology to sequence 18 tissues, aiming to uncover how K. uniflora responded to different altitudes at the gene expression level. We revealed that genes that responded to light stimuli and circadian rhythm genes were significantly enriched in DEGs in the leaf tissue group, while genes that were related to root development and peroxidase activity or involved in the pathways of cutin, suberin, wax biosynthesis, and monoterpenoid biosynthesis were significantly enriched in DEGs in the flower bud tissue group. All of the above genes may play an important role in the response of K. uniflora to various stresses, such as low temperatures and hypoxia in high-altitude environments. Furthermore, we proved that the discrepancy in gene expression patterns between leaf and flower bud tissues varied along the altitudinal gradient. Overall, our findings provide new insights into the adaptation of endangered species to high-altitude environments and further encourage parallel research to focus on the molecular mechanisms of alpine plant evolution.