Abstract: Amid the rapid increase of the global population and the quest for sustainable agriculture, the need for enhanced rice breeding strategies has become increasingly pronounced, particularly in Heilongjiang, China’s foremost rice-producing province, renowned for its premium temperate japonica rice. Here, we conducted an extensive genomic investigation of the elite rice cultivars developed in Heilongjiang Province. Using whole-genome re-sequencing of a total of 376 representative cultivars from Heilongjiang, of which 14 were developed by a single research group, we identified 4.9 million single nucleotide polymorphisms (SNPs) and 0.98 million insertions and deletions (InDels), offering a comprehensive perspective on genetic diversity and population structure. We classified the 376 rice cultivars into five subgroups based on their breeding years. Recently bred cultivars, assigned to subgroups HLJ-IV-1 and HLJ-IV-2, showed notable genetic differentiation. Through a selective sweep analysis, significant genomic variation in genes such as OsACBP5, Os4CL5, and GFR1 was pinpointed, reflecting a concerted effort in selecting for broad-spectrum disease resistance and enhanced tillering capacity. Furthermore, to identify the strengths and areas for improvement within those series, we conducted an exhaustive analysis of aromatic compounds and their corresponding genes OsODC and OsBadh2, as well as the advantageous long-grain gene OsGL3.1 haplotype within Hagengdao7. Additionally, strategies for reducing plant height through the introduction of the sd1 gene have been elucidated. With a commitment to expediting the development of superior rice cultivars, our discoveries are poised to raise the sensory attributes and nutritional profile of rice, thereby bolstering the resilience and sustainability of global food systems.

中国黑龙江省水稻品种多样性的基因组学特征研究

[1]AlexanderDH, NovembreJ, LangeK, 2009. Fast model-based estimation of ancestry in unrelated individuals. Genome Res, 19(9):1655-1664.

[2]AmarawathiY, SinghR, SinghAK, et al., 2008. Mapping of quantitative trait loci for basmati quality traits in rice (Oryza sativa L.). Mol Breeding, 21(1):49-65.

[3]AshburnerM, BallCA, BlakeJA, et al., 2000. Gene Ontology: tool for the unification of biology. Nat Genet, 25(1):25-29.

[4]BolgerAM, LohseM, UsadelB, 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 30(15):2114-2120.

[5]ChenH, PattersonN, ReichD, 2010. Population differentiation as a test for selective sweeps. Genome Res, 20(3):393-402.

[6]ChenZ, LiXX, LuHW, et al., 2020. Genomic atlases of introgression and differentiation reveal breeding footprints in Chinese cultivated rice. J Genet Genomics, 47(10):637-649.

[7]ChenZ, BuQY, LiuGF, et al., 2023. Genomic decoding of breeding history to guide breeding-by-design in rice. Natl Sci Rev, 10(5):nwad029.

[8]CingolaniP, PlattsA, WangLL, et al., 2012. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w¹¹¹⁸; iso-2; iso-3. Fly, 6(2):80-92.

[9]CutterAD, PayseurBA, 2013. Genomic signatures of selection at linked sites: unifying the disparity among species. Nat Rev Genet, 14(4):262-274.

[10]DanecekP, AutonA, AbecasisG, et al., 2011. The variant call format and VCFtools. Bioinformatics, 27(15):2156-2158.

[11]DePristoMA, BanksE, PoplinR, et al., 2011. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet, 43(5):491-498.

[12]DingCH, LinXH, ZuoY, et al., 2021. Transcription factor OsbZIP49 controls tiller angle and plant architecture through the induction of indole-3-acetic acid-amido synthetases in rice. Plant J, 108(5):1346-1364.

[13]GuZL, GongJY, ZhuZ, et al., 2023. Structure and function of rice hybrid genomes reveal genetic basis and optimal performance of heterosis. Nat Genet, 55(10):1745-1756.

[14]GuiJS, ShenJH, LiLG, 2011. Functional characterization of evolutionarily divergent 4-coumarate: coenzyme a ligases in rice. Plant Physiol, 157(2):574-586.

[15]HuangTC, TengCS, ChangJL, et al., 2008. Biosynthetic mechanism of 2-acetyl-1-pyrroline and its relationship with Δ¹-pyrroline-5-carboxylic acid and methylglyoxal in aromatic rice (Oryza sativa L.) callus. J Agric Food Chem, 56(16):7399-7404.

[16]IkedaM, MiuraK, AyaK, et al., 2013. Genes offering the potential for designing yield-related traits in rice. Curr Opin Plant Biol, 16(2):213-220.

[17]JieY, ShiTY, ZhangZ, et al., 2021. Identification of key volatiles differentiating aromatic rice cultivars using an untargeted metabolomics approach. Metabolites, 11(8):528.

[18]KanehisaM, GotoS, 2000. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res, 28(1):27-30.

[19]KhushGS, 2001. Green revolution: the way forward. Nat Rev Genet, 2(10):815-822.

[20]LiH, DurbinR, 2009. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics, 25(14):1754-1760.

[21]LiH, DurbinR, 2010. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics, 26(5):589-595.

[22]LiH, HandsakerB, WysokerA, et al., 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics, 25(16):2078-2079.

[23]LiY, ZhangWT, LiMY, et al., 2024. Discovery of OsODC as a key enhancer of aroma and development of highly fragrant rice. Plant Commun, 6(1):101141.

[24]LiuCX, PengP, LiWG, et al., 2021. Deciphering variation of 239 elite japonica rice genomes for whole genome sequences-enabled breeding. Genomics, 113(5):3083-3091.

[25]LiuEB, ZengSY, ZhuSS, et al., 2019. Favorable alleles of GRAIN-FILLING RATE1 increase the grain-filling rate and yield of rice. Plant Physiol, 181(3):1207-1222.

[26]LiuJF, ChenJ, ZhengXM, et al., 2017. GW5 acts in the brassinosteroid signalling pathway to regulate grain width and weight in rice. Nat Plants, 3(5):17043.

[27]LiuWZ, WuC, FuYP, et al., 2009. Identification and characterization of HTD2: a novel gene negatively regulating tiller bud outgrowth in rice. Planta, 230(4):649-658.

[28]LongY, WangC, LiuC, et al., 2024. Molecular mechanisms controlling grain size and weight and their biotechnological breeding applications in maize and other cereal crops. J Adv Res, 62:27-46.

[29]LorieuxM, PetrovM, HuangN, et al., 1996. Aroma in rice: genetic analysis of a quantitative trait. Theor Appl Genet, 93(7):1145-1151.

[30]NarayananSP, LungSC, LiaoP, et al., 2020. The overexpression of OsACBP5 protects transgenic rice against necrotrophic, hemibiotrophic and biotrophic pathogens. Sci Rep, 10:14918.

[31]PachauriV, MishraV, MishraP, et al., 2014. Identification of candidate genes for rice grain aroma by combining QTL mapping and transcriptome profiling approaches. Cereal Res Commun, 42(3):376-388.

[32]PengH, WangK, ChenZ, et al., 2020. MBKbase for rice: an integrated omics knowledgebase for molecular breeding in rice. Nucleic Acids Res, 48(D1):D1085-D1092.

[33]PengJF, ZhuY, LinF, et al., 2023. Direct determination of 2-acetyl-1-pyrroline in rice by ultrasound-assisted solvent extraction coupled with ultra-performance liquid chromatography-tandem mass spectrometry. Food Anal Method, 16(5):900-908.

[34]ProdhanZH, ShuQY, 2020. Rice aroma: a natural gift comes with price and the way forward. Rice Sci, 27(2):86-100.

[35]PurcellS, NealeB, Todd-BrownK, et al., 2007. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 81(3):559-575.

[36]RajA, StephensM, PritchardJK, 2014. fastSTRUCTURE: variational inference of population structure in large SNP data sets. Genetics, 197(2):573-589.

[37]ReichD, PriceAL, PattersonN, 2008. Principal component analysis of genetic data. Nat Genet, 40(5):491-492.

[38]RetiefJD, 2000. Phylogenetic analysis using PHYLIP. In: Misener S, Krawetz SA (Eds.), Bioinformatics Methods and Protocols. Humana Press, Totowa, p.243-258.

[39]RyooN, YuC, ParkCS, et al., 2007. Knockout of a starch synthase gene OsSSIIIa/Flo5 causes white-core floury endosperm in rice (Oryza sativa L.). Plant Cell Rep, 26(7):1083-1095.

[40]SakthivelK, SundaramRM, RaniNS, et al., 2009. Genetic and molecular basis of fragrance in rice. Biotechnol Adv, 27(4):468-473.

[41]ShangLG, HeWC, WangTY, et al., 2023. A complete assembly of the rice Nipponbare reference genome. Mol Plant, 16(8):1232-1236.

[42]ShiCL, DongNQ, GuoT, et al., 2020. A quantitative trait locus GW6 controls rice grain size and yield through the gibberellin pathway. Plant J, 103(3):1174-1188.

[43]TalukdarPR, RathiS, PathakK, et al., 2017. Population structure and marker-trait association in indigenous aromatic rice. Rice Sci, 24(3):145-154.

[44]TeraoT, HiroseT, 2015. Control of grain protein contents through SEMIDWARF1 mutant alleles: sd1 increases the grain protein content in Dee-geo-woo-gen but not in Reimei. Mol Genet Genomics, 290(3):939-954.

[45]ValentB, 2021. The impact of blast disease: past, present, and future. In: Jacob S (Ed.), Magnaporthe Oryzae. Humana, New York, p.1-18.

[46]WangC, FengXM, YuanQB, et al., 2023. Upgrading the genome of an elite japonica rice variety Kongyu 131 for lodging resistance improvement. Plant Biotechnol J, 21(2):419-432.

[47]WangWS, MauleonR, HuZQ, et al., 2018. Genomic variation in 3010 diverse accessions of Asian cultivated rice. Nature, 557(7703):43-49.

[48]WangYX, HuangLY, DuFP, et al., 2021. Comparative transcriptome and metabolome profiling reveal molecular mechanisms underlying OsDRAP1-mediated salt tolerance in rice. Sci Rep, 11:5166.

[49]WuTZ, HuEQ, XuSB, et al., 2021. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (Camb), 2(3):100141.

[50]XinFF, XiaoXM, DongJW, et al., 2020. Large increases of paddy rice area, gross primary production, and grain production in Northeast China during 2000–2017. Sci Total Environ, 711:135183.

[51]YangJ, LeeSH, GoddardME, et al., 2011. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet, 88(1):76-82.

[52]YangS, LiuMS, ChuN, et al., 2022. Combined transcriptome and metabolome reveal glutathione metabolism plays a critical role in resistance to salinity in rice landraces HD961. Front Plant Sci, 13:952595.

[53]YeJH, ZhangMC, YuanXP, et al., 2022. Genomic insight into genetic changes and shaping of major inbred rice cultivars in China. New Phytol, 236(6):2311-2326.

[54]YouNS, DongJW, HuangJX, et al., 2021. The 10-m crop type maps in Northeast China during 2017–2019. Sci Data, 8:41.

[55]ZhangC, DongSS, XuJY, et al., 2019. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics, 35(10):1786-1788.

[56]ZhengXM, WeiF, ChengC, et al., 2024. A historical review of hybrid rice breeding. J Integr Plant Biol, 66(3):532-545.