MXenes with ordered triatomic-layer borate polyanion terminations.

Surface terminations profoundly influence the intrinsic properties of MXenes, but existing terminations are limited to monoatomic layers or simple groups, showing disordered arrangements and inferior stability. Here we present the synthesis of MXenes with triatomic-layer borate polyanion terminations (OBO terminations) through a flux-assisted eutectic molten etching approach. During the synthesis, Lewis acidic salts act as the etching agent to obtain the MXene backbone, while borax generates BO2- species, which cap the MXene surface with an O-B-O configuration. In contrast to conventional chlorine/oxygen-terminated Nb2C with localized charge transport, OBO-terminated Nb2C features band transport described by the Drude model, exhibiting a 15-fold increase in electrical conductivity and a 10-fold improvement in charge mobility at the d.c. limit. This transition is attributed to surface ordering that effectively mitigates charge carrier backscattering and trapping. Additionally, OBO terminations provide Ti3C2 MXene with substantially enriched Li+-hosting sites and thereby a large charge-storage capacity of 420 mAh g-1. Our findings illustrate the potential of intricate termination configurations in MXenes and their applications for (opto)electronics and energy storage.

Upland rice genomic signatures of adaptation to drought resistance and navigation to molecular design breeding.

Upland rice is a distinctive drought-aerobic ecotype of cultivated rice highly resistant to drought stress. However, the genetic and genomic basis for the drought-aerobic adaptation of upland rice remains largely unclear due to the lack of genomic resources. In this study, we identified 25 typical upland rice accessions and assembled a high-quality genome of one of the typical upland rice varieties, IRAT109, comprising 384 Mb with a contig N50 of 19.6 Mb. Phylogenetic analysis revealed upland and lowland rice have distinct ecotype differentiation within the japonica subgroup. Comparative genomic analyses revealed that adaptive differentiation of lowland and upland rice is likely attributable to the natural variation of many genes in promoter regions, formation of specific genes in upland rice, and expansion of gene families. We revealed differentiated gene expression patterns in the leaves and roots of the two ecotypes and found that lignin synthesis mediated by the phenylpropane pathway plays an important role in the adaptive differentiation of upland and lowland rice. We identified 28 selective sweeps that occurred during domestication and validated that the qRT9 gene in selective regions can positively regulate drought resistance in rice. Eighty key genes closely associated with drought resistance were appraised for their appreciable potential in drought resistance breeding. Our study enhances the understanding of the adaptation of upland rice and provides a genome navigation map of drought resistance breeding, which will facilitate the breeding of drought-resistant rice and the "blue revolution" in agriculture.

The wall-associated kinase GWN1 controls grain weight and grain number in rice.

Grain size is a crucial agronomic trait that determines grain weight and final yield. Although several genes have been reported to regulate grain size in rice (Oryza sativa), the function of Wall-Associated Kinase family genes affecting grain size is still largely unknown. In this study, we identified GRAIN WEIGHT AND NUMBER 1 (GWN1) using map-based cloning. GWN1 encodes the OsWAK74 protein kinase, which is conserved in plants. GWN1 negatively regulates grain length and weight by regulating cell proliferation in spikelet hulls. We also found that GWN1 negatively influenced grain number by influencing secondary branch numbers and finally increased plant grain yield. The GWN1 gene was highly expressed in inflorescences and its encoded protein is located at the cell membrane and cell wall. Moreover, we identified three haplotypes of GWN1 in the germplasm. GWN1hap1 showing longer grain, has not been widely utilized in modern rice varieties. In summary, GWN1 played a very important role in regulating grain length, weight and number, thereby exhibiting application potential in molecular breeding for longer grain and higher yield.

Genomic variation in 3,010 diverse accessions of Asian cultivated rice.

Here we analyse genetic variation, population structure and diversity among 3,010 diverse Asian cultivated rice (Oryza sativa L.) genomes from the 3,000 Rice Genomes Project. Our results are consistent with the five major groups previously recognized, but also suggest several unreported subpopulations that correlate with geographic location. We identified 29 million single nucleotide polymorphisms, 2.4 million small indels and over 90,000 structural variations that contribute to within- and between-population variation. Using pan-genome analyses, we identified more than 10,000 novel full-length protein-coding genes and a high number of presence-absence variations. The complex patterns of introgression observed in domestication genes are consistent with multiple independent rice domestication events. The public availability of data from the 3,000 Rice Genomes Project provides a resource for rice genomics research and breeding.

Changes in the soil fungal communities of steppe grasslands at varying degradation levels in North China.

The grasslands in North China are rich in fungal resources. However, the knowledge of the structure and function of fungal communities and the role of microbial communities in vegetation restoration and succession are limited. Thus, we used an Illumina HiSeq PE250 high-throughput sequencing platform to study the changing characteristics of soil fungal communities in degraded grasslands, which were categorized as non-degraded (ND), lightly degraded, moderately degraded, and severely degraded (SD). Moreover, a correlation analysis between soil physical and chemical properties and fungal communities was completed. The results showed that the number of plant species, vegetation coverage, aboveground biomass, and diversity index decreased significantly with increasing degradation, and there were significant differences in the physical and chemical properties of the soil among the different degraded grasslands. The dominant fungal phyla in the degraded grassland were as follows: Ascomycota, 44.88%-65.03%; Basidiomycota, 12.68%-29.91%; and unclassified, 5.51%-16.91%. The dominant fungi were as follows: Mortierella, 6.50%-11.41%; Chaetomium, 6.71%-11.58%; others, 25.95%-36.14%; and unclassified, 25.56%-53.0%. There were significant differences in the microbial Shannon-Wiener and Chao1 indices between the ND and degraded meadows, and the composition and diversity of the soil fungal community differed significantly as the meadows continued to deteriorate. The results showed that pH was the most critical factor affecting soil microbial and fungal communities in SD grasslands, whereas soil microbial and fungal communities in ND grasslands were mainly affected by water content and other environmental factors.

Natural variation in MORE GRAINS 1 regulates grain number and grain weight in rice.

Grain yield is determined mainly by grain number and grain weight. In this study, we identified and characterized MORE GRAINS1 (MOG1), a gene associated with grain number and grain weight in rice (Oryza sativa L.), through map-based cloning. Overexpression of MOG1 increased grain yield by 18.6%-22.3% under field conditions. We determined that MOG1, a bHLH transcription factor, interacts with OsbHLH107 and directly activates the expression of LONELY GUY (LOG), which encodes a cytokinin-activating enzyme and the cell expansion gene EXPANSIN-LIKE1 (EXPLA1), positively regulating grain number per panicle and grain weight. Natural variations in the promoter and coding regions of MOG1 between Hap-LNW and Hap-HNW alleles resulted in changes in MOG1 expression level and transcriptional activation, leading to functional differences. Haplotype analysis revealed that Hap-HNW, which results in a greater number and heavier grains, has undergone strong selection but has been poorly utilized in modern lowland rice breeding. In summary, the MOG1-OsbHLH107 complex activates LOG and EXPLA1 expression to promote cell expansion and division of young panicles through the cytokinin pathway, thereby increasing grain number and grain weight. These findings suggest that Hap-HNW could be used in strategies to breed high-yielding temperate japonica lowland rice.

Cold-adaptive evolution at the reproductive stage in Geng/japonica subspecies reveals the role of OsMAPK3 and OsLEA9.

Cold stress at the reproductive stage severely affects the production and geographic distribution of rice. The Geng/japonica subpopulation gradually developed stronger cold adaptation than the Xian/indica subpopulation during the long-term domestication of cultivated rice. However, the evolutionary path and natural alleles underlying the cold adaptability of intra-Geng subspecies remain largely unknown. Here, we identified MITOGEN-ACTIVATED PROTEIN KINASE 3 (OsMAPK3) and LATE EMBRYOGENESIS ABUNDANT PROTEIN 9 (OsLEA9) as two important regulators for the cold adaptation of Geng subspecies from a combination of transcriptome analysis and genome-wide association study. Transgenic validation showed that OsMAPK3 and OsLEA9 confer cold tolerance at the reproductive stage. Selection and evolution analysis suggested that the Geng version of OsMAPK3 (OsMAPK3Geng ) directly evolved from Chinese Oryza rufipogon III and was largely retained in high-latitude and high-altitude regions with low temperatures during domestication. Later, the functional nucleotide polymorphism (FNP-776) in the Kunmingxiaobaigu and Lijiangxiaoheigu version of the OsLEA9 (OsLEA9KL ) promoter originated from novel variation of intra-Geng was selected and predominantly retained in temperate Geng to improve the adaptation of Geng together with OsMAPK3Geng to colder climatic conditions in high-latitude areas. Breeding potential analysis suggested that pyramiding of OsMAPK3Geng and OsLEA9KL enhanced the cold tolerance of Geng and promotes the expansion of cultivated rice to colder regions. This study not only highlights the evolutionary path taken by the cold-adaptive differentiation of intra-Geng, but also provides new genetic resources for rice molecular breeding in low-temperature areas.

Semiconducting Conjugated Coordination Polymer with High Charge Mobility Enabled by "4 + 2" Phenyl Ligands.

Electrically conductive coordination polymers and metal-organic frameworks are attractive emerging electroactive materials for (opto-)electronics. However, developing semiconducting coordination polymers with high charge carrier mobility for devices remains a major challenge, urgently requiring the rational design of ligands and topological networks with desired electronic structures. Herein, we demonstrate a strategy for synthesizing high-mobility semiconducting conjugated coordination polymers (c-CPs) utilizing novel conjugated ligands with D2h symmetry, namely, "4 + 2" phenyl ligands. Compared with the conventional phenyl ligands with C6h symmetry, the reduced symmetry of the "4 + 2" ligands leads to anisotropic coordination in the formation of c-CPs. Consequently, we successfully achieve a single-crystalline three-dimensional (3D) c-CP Cu4DHTTB (DHTTB = 2,5-dihydroxy-1,3,4,6-tetrathiolbenzene), containing orthogonal ribbon-like π-d conjugated chains rather than 2D conjugated layers. DFT calculation suggests that the resulting Cu4DHTTB exhibits a small band gap (∼0.2 eV), strongly dispersive energy bands near the Fermi level with a low electron-hole reduced effective mass (∼0.2m0*). Furthermore, the four-probe method reveals a semiconducting behavior with a decent conductivity of 0.2 S/cm. Thermopower measurement suggests that it is a p-type semiconductor. Ultrafast terahertz photoconductivity measurements confirm Cu4DHTTB's semiconducting nature and demonstrate the Drude-type transport with high charge carrier mobilities up to 88 ± 15 cm2 V-1 s-1, outperforming the conductive 3D coordination polymers reported till date. This molecular design strategy for constructing high-mobility semiconducting c-CPs lays the foundation for achieving high-performance c-CP-based (opto-)electronics.

Genetic-based dissection of resistance to bacterial leaf streak in rice by GWAS.

BACKGROUND: Rice is the second-largest food crop in the world and vulnerable to bacterial leaf streak disease. A thorough comprehension of the genetic foundation of agronomic traits was essential for effective implementation of molecular marker-assisted selection. RESULTS: Our study aimed to evaluate the vulnerability of rice to bacterial leaf streak disease (BLS) induced by the gram-negative bacterium Xanthomonas oryzae pv. oryzicola (Xoc). In order to accomplish this, we first analyzed the population structure of 747 accessions and subsequently assessed their phenotypes 20 days after inoculation with a strain of Xoc, GX01. We conducted genome-wide association studies (GWAS) on a population of 747 rice accessions, consisting of both indica and japonica subpopulations, utilizing phenotypic data on resistance to bacterial leaf streak (RBLS) and sequence data. We identified a total of 20 QTLs associated with RBLS in our analysis. Through the integration of linkage mapping, sequence analysis, haplotype analysis, and transcriptome analysis, we were able to identify five potential candidate genes (OsRBLS1-OsRBLS5) that possess the potential to regulate RBLS in rice. In order to gain a more comprehensive understanding of the genetic mechanism behind resistance to bacterial leaf streak, we conducted tests on these genes in both the indica and japonica subpopulations, ultimately identifying superior haplotypes that suggest the potential utilization of these genes in breeding disease-resistant rice varieties. CONCLUSIONS: The findings of our study broaden our comprehension of the genetic mechanisms underlying RBLS in rice and offer significant insights that can be applied towards genetic improvement and breeding of disease-resistant rice in rapidly evolving environmental conditions.

Natural variation in Tiller Number 1 affects its interaction with TIF1 to regulate tillering in rice.

Tiller number per plant-a cardinal component of ideal plant architecture-affects grain yield potential. Thus, alleles positively affecting tillering must be mined to promote genetic improvement. Here, we report a Tiller Number 1 (TN1) protein harbouring a bromo-adjacent homology domain and RNA recognition motifs, identified through genome-wide association study of tiller numbers. Natural variation in TN1 affects its interaction with TIF1 (TN1 interaction factor 1) to affect DWARF14 expression and negatively regulate tiller number in rice. Further analysis of variations in TN1 among indica genotypes according to geographical distribution revealed that low-tillering varieties with TN1-hapL are concentrated in Southeast Asia and East Asia, whereas high-tillering varieties with TN1-hapH are concentrated in South Asia. Taken together, these results indicate that TN1 is a tillering regulatory factor whose alleles present apparent preferential utilization across geographical regions. Our findings advance the molecular understanding of tiller development.

RRS1 shapes robust root system to enhance drought resistance in rice.

A strong root system facilitates the absorption of water and nutrients from the soil, to improve the growth of crops. However, to date, there are still very few root development regulatory genes that can be used in crop breeding for agriculture. In this study, we cloned a negative regulator gene of root development, Robust Root System 1 (RRS1), which encodes an R2R3-type MYB family transcription factor. RRS1 knockout plants showed enhanced root growth, including longer root length, longer lateral root length, and larger lateral root density. RRS1 represses root development by directly activating the expression of OsIAA3 which is involved in the auxin signaling pathway. A natural variation in the coding region of RRS1 changes the transcriptional activity of its protein. RRS1T allele, originating from wild rice, possibly increases root length by means of weakening regulation of OsIAA3. Knockout of RRS1 enhances drought resistance by promoting water absorption and improving water use efficiency. This study provides a new gene resource for improving root systems and cultivating drought-resistant rice varieties with important values in agricultural applications.

Variations in OsSPL10 confer drought tolerance by directly regulating OsNAC2 expression and ROS production in rice.

Drought is a major factor restricting the production of rice (Oryza sativa L.). The identification of natural variants for drought stress-related genes is an important step toward developing genetically improved rice varieties. Here, we characterized a member of the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) family, OsSPL10, as a transcription factor involved in the regulation of drought tolerance in rice. OsSPL10 appears to play a vital role in drought tolerance by controlling reactive oxygen species (ROS) production and stomatal movements. Haplotype and allele frequency analyses of OsSPL10 indicated that most upland rice and improved lowland rice varieties harbor the OsSPL10Hap1 allele, whereas the OsSPL10Hap2 allele was mainly present in lowland and landrace rice varieties. Importantly, we demonstrated that the varieties with the OsSPL10Hap1 allele showed low expression levels of OsSPL10 and its downstream gene, OsNAC2, which decreases the expression of OsAP37 and increases the expression of OsCOX11, thus preventing ROS accumulation and programmed cell death (PCD). Furthermore, the knockdown or knockout of OsSPL10 induced fast stomatal closure and prevented water loss, thereby improving drought tolerance in rice. Based on these observations, we propose that OsSPL10 confers drought tolerance by regulating OsNAC2 expression and that OsSPL10Hap1 could be a valuable haplotype for the genetic improvement of drought tolerance in rice.

RGN1 controls grain number and shapes panicle architecture in rice.

Yield in rice is determined mainly by panicle architecture. Using map-based cloning, we identified an R2R3 MYB transcription factor REGULATOR OF GRAIN NUMBER1 (RGN1) affecting grain number and panicle architecture. Mutation of RGN1 caused an absence of lateral grains on secondary branches. We demonstrated that RGN1 controls lateral grain formation by regulation of LONELY GUY (LOG) expression, thus controlling grain number and shaping panicle architecture. A novel favourable allele, RGN1C , derived from the Or-I group in wild rice affected panicle architecture by means longer panicles. Identification of RGN1 provides a theoretical basis for understanding the molecular mechanism of lateral grain formation in rice; RGN1 will be an important gene resource for molecular breeding for higher yield.

Uncovering Hierarchical Regulation among MYB-bHLH-WD40 Proteins and Manipulating Anthocyanin Pigmentation in Rice.

Anthocyanins accumulate in various organs of rice, and the regulatory genes involved in pigmentation of specific organs, such as pericarp, hull, leaf, apiculus, and stigma have been elucidated. However, the corresponding gene for rice culm pigmentation has not been clarified. The well-known MYB-bHLH-WD40 (MBW) complex plays vital role in regulating the anthocyanin biosynthesis pathway in plants. However, the core members of MBW and the hierarchical regulation between these members are not fully elucidated in rice. Here, by map-based cloning, we identified the culm-specific pigmentation gene S1 whose alleles are also known for hull/pericarp pigmentation. We also clarified that one WD40 protein encoding gene, WA1, is indispensable for anthocyanin biosynthesis in rice. In the cascading regulation among MBW members, S1 (bHLH) acts as the master gene by activating the expression of C1 (MYB), and then C1 activates the expression of WA1 (WD40), which is unique in plant species. This enables MBW members to be coordinated in a common way to efficiently regulate anthocyanin biosynthesis genes. Based on these studies, we explored the minimal gene set required for anthocyanin biosynthesis in rice. These findings will help us design new rice varieties with anthocyanin accumulation in specific organs as needed.

Effects of air pollutants and temperature on the number of asthma outpatient visits in Hohhot, China.

Although numerous studies have linked asthma to air temperature and pollution, few studies have examined their interactive effects on asthma outpatient visits. This study investigated how air pollutants and their interactions with temperature affect asthma outpatient visits in a city in northern Chinaduring the time period 2018 - 2020 . . As the results, 24,163 asthma outpatients were recorded, a 10-µg/m3 increase in PM10, PM2.5, SO2, and NO2 concentrations was associated with significant increases in visits of 3.47% (95% CI: 2.35%-4.60%), 0.83% (95% CI: 0.36%-1.30%), 3.17% (95% CI: 1.47%-4.90%), and 8.90% (95% CI: 6.09%-11.79%), respectively. The effect was stronger in females than males, and stronger in the elderly (≥65 years) than among the young. The interaction between low temperatures and high air pollution levels significantly increased the number of asthma outpatient visits. This study emphasizesthe importance of reducing air pollution in order to lessen the effects of cold.

Genetic architecture to cause dynamic change in tiller and panicle numbers revealed by genome-wide association study and transcriptome profile in rice.

Panicle number (PN) is one of the three yield components in rice. As one of the most unstable traits, the dynamic change in tiller number (DCTN) may determine the final PN. However, the genetic basis of DCTN and its relationship with PN remain unclear. Here, 377 deeply re-sequenced rice accessions were used to perform genome-wide association studies (GWAS) for tiller/PN. It was found that the DCTN pattern rather than maximum tiller number or effective tiller ratio is the determinant factor of high PN. The DCTN pattern that affords more panicles exhibits a period of stable tillering peak between 30 and 45 days after transplant (called DT30 and DT45, respectively), which was believed as an ideal pattern contributing to the steady transition from tiller development to panicle development (ST-TtP). Consistently, quantitative trait loci (QTL) expressed near DT30-DT45 were especially critical to the rice DCTN and in supporting the ST-TtP. The spatio-temporal expression analysis showed that the expression pattern of keeping relatively high expression in root at 24:00 (R24-P2) from about DT30 to DT45 is a typical expression pattern of cloned tiller genes, and the candidate genes with R24-P2 can facilitate the prediction of PN. Moreover, gene OsSAUR27 was identified by an integrated approach combining GWAS, bi-parental QTL mapping and transcription. These findings related to the genetic basis underlying the DCTN will provide the genetic theory in making appropriate decisions on field management, and in developing new varieties with high PN and ideal dynamic plant architecture.

Defective Nanographenes Containing Seven-Five-Seven (7-5-7)-Membered Rings.

Defects have been observed in graphene and are expected to play a key role in its optical, electronic, and magnetic properties. However, because most of the studies focused on the structural characterization, the implications of topological defects on the physicochemical properties of graphene remain poorly understood. Here, we demonstrate a bottom-up synthesis of three novel nanographenes (1-3) with well-defined defects in which seven-five-seven (7-5-7)-membered rings were introduced to their sp2 carbon frameworks. From the X-ray crystallographic analysis, compound 1 adopts a nearly planar structure. Compound 2, with an additional five-membered ring compared to 1, possesses a slightly saddle-shaped geometry. Compound 3, which can be regarded as the "head-to-head" fusion of 1 with two bonds, features two saddles connected together. The resultant defective nanographenes 1-3 were well-investigated by UV-vis absorption, cyclic voltammetry, and time-resolved absorption spectra and further corroborated by density functional theory (DFT) calculations. Detailed experimental and theoretical investigations elucidate that these three nanographenes 1-3 exhibit an anti-aromatic character in their ground states and display a high stability under ambient conditions, which contrast with the reported unstable biradicaloid nanographenes that contain heptagons. Our work reported herein offers insights into the understanding of structure-related properties and enables the control of the electronic structures of expanded nanographenes with atomically precise defects.

Dual-Redox-Sites Enable Two-Dimensional Conjugated Metal-Organic Frameworks with Large Pseudocapacitance and Wide Potential Window.

Advanced supercapacitor electrodes require the development of materials with dense redox sites embedded into conductive and porous skeletons. Two-dimensional (2D) conjugated metal-organic frameworks (c-MOFs) are attractive supercapacitor electrode materials due to their high intrinsic electrical conductivities, large specific surface areas, and quasi-one-dimensional aligned pore arrays. However, the reported 2D c-MOFs still suffer from unsatisfying specific capacitances and narrow potential windows because large and redox-inactive building blocks lead to low redox-site densities of 2D c-MOFs. Herein, we demonstrate the dual-redox-site 2D c-MOFs with copper phthalocyanine building blocks linked by metal-bis(iminobenzosemiquinoid) (M2[CuPc(NH)8], M = Ni or Cu), which depict both large specific capacitances and wide potential windows. Experimental results accompanied by theoretical calculations verify that phthalocyanine monomers and metal-bis(iminobenzosemiquinoid) linkages serve as respective redox sites for pseudocapacitive cation (Na+) and anion (SO42-) storage, enabling the continuous Faradaic reactions of M2[CuPc(NH)8] occurring in a large potential window of -0.8 to 0.8 V vs Ag/AgCl (3 M KCl). The decent conductivity (0.8 S m-1) and high active-site density further endow the Ni2[CuPc(NH)8] with a remarkable specific capacitance (400 F g-1 at 0.5 A g-1) and excellent rate capability (183 F g-1 at 20 A g-1). Quasi-solid-state symmetric supercapacitors are further assembled to demonstrate the practical application of Ni2[CuPc(NH)8] electrode, which deliver a state-of-the-art energy density of 51.6 Wh kg-1 and a peak power density of 32.1 kW kg-1.

Migration Kinetics of Surface Ions in Oxygen-Deficient Perovskite During Topotactic Transitions.

Oxygen diffusivity and surface exchange kinetics underpin the ionic, electronic, and catalytic functionalities of complex multivalent oxides. Towards understanding and controlling the kinetics of oxygen transport in emerging technologies, it is highly desirable to reveal the underlying lattice dynamics and ionic activities related to oxygen variation. In this study, the evolution of oxygen content is identified in real-time during the progress of a topotactic phase transition in La0.7 Sr0.3 MnO3-δ epitaxial thin films, both at the surface and throughout the bulk. Using polarized neutron reflectometry, a quantitative depth profile of the oxygen content gradient is achieved, which, alongside atomic-resolution scanning transmission electron microscopy, uniquely reveals the formation of a novel structural phase near the surface. Surface-sensitive X-ray spectroscopies further confirm a significant change of the electronic structure accompanying the transition. The anisotropic features of this novel phase enable a distinct oxygen diffusion pathway in contrast to conventional observation of oxygen motion at moderate temperatures. The results provide insights furthering the design of solid oxygen ion conductors within the framework of topotactic phase transitions.

Stepwise selection of natural variations at CTB2 and CTB4a improves cold adaptation during domestication of japonica rice.

The improvement of cold adaptation has contributed to the increased growing area of rice. Standing variation and de novo mutation are distinct natural sources of beneficial alleles in plant adaptation. However, the genetic mechanisms and evolutionary patterns underlying these sources in a single population during crop domestication remain elusive. Here we cloned the CTB2 gene, encoding a UDP-glucose sterol glucosyltransferase, for cold tolerance in rice at the booting stage. A single standing variation (I408V) in the conserved UDPGT domain of CTB2 originated from Chinese Oryza rufipogon and contributed to the cold adaptation of Oryza sativa ssp. japonica. CTB2 is located in a 56.8 kb region, including the previously reported gene CTB4a in which de novo mutation arose c. 3200 yr BP in Yunnan province, China, conferring cold tolerance. Standing variation of CTB2 and de novo mutation of CTB4a underwent stepwise selection to facilitate cold adaptation to expand rice cultivation from high-altitude to high-latitude regions. These results provide an example of stepwise selection on two kinds of variation and describe a new molecular mechanism of cold adaptation in japonica rice.