Lost genome segments associate with trait diversity during rice domestication.

BACKGROUND: DNA mutations of diverse types provide the raw material required for phenotypic variation and evolution. In the case of crop species, previous research aimed to elucidate the changing patterns of repetitive sequences, single-nucleotide polymorphisms (SNPs), and small InDels during domestication to explain morphological evolution and adaptation to different environments. Additionally, structural variations (SVs) encompassing larger stretches of DNA are more likely to alter gene expression levels leading to phenotypic variation affecting plant phenotypes and stress resistance. Previous studies on SVs in rice were hampered by reliance on short-read sequencing limiting the quantity and quality of SV identification, while SV data are currently only available for cultivated rice, with wild rice largely uncharacterized. Here, we generated two genome assemblies for O. rufipogon using long-read sequencing and provide insights on the evolutionary pattern and effect of SVs on morphological traits during rice domestication. RESULTS: In this study, we identified 318,589 SVs in cultivated and wild rice populations through a comprehensive analysis of 13 high-quality rice genomes and found that wild rice genomes contain 49% of unique SVs and an average of 1.76% of genes were lost during rice domestication. These SVs were further genotyped for 649 rice accessions, their evolutionary pattern during rice domestication and potential association with the diversity of important agronomic traits were examined. Genome-wide association studies between these SVs and nine agronomic traits identified 413 candidate causal variants, which together affect 361 genes. An 824-bp deletion in japonica rice, which encodes a serine carboxypeptidase family protein, is shown to be associated with grain length. CONCLUSIONS: We provide relatively accurate and complete SV datasets for cultivated and wild rice accessions, especially in TE-rich regions, by comparing long-read sequencing data for 13 representative varieties. The integrated rice SV map and the identified candidate genes and variants represent valuable resources for future genomic research and breeding in rice.

LncPheDB: a genome-wide lncRNAs regulated phenotypes database in plants.

LncPheDB (https://www.lncphedb.com/) is a systematic resource of genome-wide long non-coding RNAs (lncRNAs)-phenotypes associations for multiple species. It was established to display the genome-wide lncRNA annotations, target genes prediction, variant-trait associations, gene-phenotype correlations, lncRNA-phenotype correlations, and the similar non-coding regions of the queried sequence in multiple species. LncPheDB sorted out a total of 203,391 lncRNA sequences, 2000 phenotypes, and 120,271 variants of nine species (Zea mays L., Gossypium barbadense L., Triticum aestivum L., Lycopersicon esculentum Mille, Oryza sativa L., Hordeum vulgare L., Sorghum bicolor L., Glycine max L., and Cucumis sativus L.). By exploring the relationship between lncRNAs and the genomic position of variants in genome-wide association analysis, a total of 68,862 lncRNAs were found to be related to the diversity of agronomic traits. More importantly, to facilitate the study of the functions of lncRNAs, we analyzed the possible target genes of lncRNAs, constructed a blast tool for performing similar fragmentation studies in all species, linked the pages of phenotypic studies related to lncRNAs that possess similar fragments and constructed their regulatory networks. In addition, LncPheDB also provides a user-friendly interface, a genome visualization platform, and multi-level and multi-modal convenient data search engine. We believe that LncPheDB plays a crucial role in mining lncRNA-related plant data. Supplementary Information: The online version contains supplementary material available at 10.1007/s42994-022-00084-3.

Identification of candidate genes and clarification of the maintenance of the green pericarp of weedy rice grains.

The weedy rice (Oryza sativa f. spontanea) pericarp has diverse colors (e.g., purple, red, light-red, and white). However, research on pericarp colors has focused on red and purple, but not green. Unlike many other common weedy rice resources, LM8 has a green pericarp at maturity. In this study, the coloration of the LM8 pericarp was evaluated at the cellular and genetic levels. First, an examination of their ultrastructure indicated that LM8 chloroplasts were normal regarding plastid development and they contained many plastoglobules from the early immature stage to maturity. Analyses of transcriptome profiles and differentially expressed genes revealed that most chlorophyll (Chl) degradation-related genes in LM8 were expressed at lower levels than Chl a/b cycle-related genes in mature pericarps, suggesting that the green LM8 pericarp was associated with inhibited Chl degradation in intact chloroplasts. Second, the F2 generation derived from a cross between LM8 (green pericarp) and SLG (white pericarp) had a pericarp color segregation ratio of 9:3:4 (green:brown:white). The bulked segregant analysis of the F2 populations resulted in the identification of 12 known genes in the chromosome 3 and 4 hotspot regions as candidate genes related to Chl metabolism in the rice pericarp. The RNA-seq and sqRT-PCR assays indicated that the expression of the Chl a/b cycle-related structural gene DVR (encoding divinyl reductase) was sharply up-regulated. Moreover, genes encoding magnesium-chelatase subunit D and the light-harvesting Chl a/b-binding protein were transcriptionally active in the fully ripened dry pericarp. Regarding the ethylene signal transduction pathway, the CTR (encoding an ethylene-responsive protein kinase) and ERF (encoding an ethylene-responsive factor) genes expression profiles were determined. The findings of this study highlight the regulatory roles of Chl biosynthesis- and degradation-related genes influencing Chl accumulation during the maturation of the LM8 pericarp.

Genomic Analysis Provides Insights Into the Plant Architecture Variations in in situ Conserved Chinese Wild Rice (Oryza rufipogon Griff.).

In situ conserved wild rice (Oryza rufipogon Griff.) is a promising source of alleles for improving rice production worldwide. In this study, we conducted a genomic analysis of an in situ conserved wild rice population (Guiping wild rice) growing at the center of wild rice genetic diversity in South China. Differences in the plant architecture in this population were investigated. An analysis using molecular markers revealed the substantial genetic diversity in this population, which was divided into subgroups according to the plant architecture. After resequencing representative individuals, the Guiping wild rice population was compared with other O. rufipogon and Oryza sativa populations. The results indicated that this in situ conserved wild rice population has a unique genetic structure, with genes that were introgressed from aromatic and O. sativa ssp. indica and japonica populations. The QTLs associated with plant architecture in this population were detected via a pair-wise comparison analysis of the sequencing data for multiple DNA pools. The results suggested that a heading date-related gene (DHD1) might be associated with variations in plant architecture and may have originated in cultivated rice. Our findings provide researchers with useful insights for future genomic analyses of in situ conserved wild rice populations.

Genome-wide selection and introgression of Chinese rice varieties during breeding.

China is the largest rice-producing country, but the genomic landscape of rice diversity has not yet been clarified. In this study, we re-sequence 1070 rice varieties collected from China (400) and other regions in Asia (670). Among the six major rice groups (aus, indica-I, indica-II, aromatic, temperate japonica, and tropical japonica), almost all Chinese varieties belong to the indica-II or temperate japonica group. Most Chinese indica varieties belong to indica-II, which consists of two subgroups developed during different phases of rice breeding. The genomic segments underlying the differences between these subgroups span 36.32 Mb. The Chinese japonica rice varieties fall into the temperate japonica group, consisting of two subgroups based on their geographical distribution. The genomic segments underlying the differences between these subgroups span 27.69 Mb. These differentiated segments in the Chinese indica varieties span 45 genes with nonsynonymous mutations that are closely related to variations in plant height and grain width. Fifty-four genes with nonsynonymous mutations are associated with the differences in heading date between the two Chinese japonica subgroups. These findings provide new insights into rice diversity in China that will facilitate the molecular breeding.

Global whole-genome comparison and analysis to classify subpopulations and identify resistance genes in weedy rice relevant for improving crops.

Common weedy rice plants are important genetic resources for modern breeding programs because they are the closest relatives to rice cultivars and their genomes contain elite genes. Determining the utility and copy numbers of WRKY and nucleotide-binding site (NBS) resistance-related genes may help to clarify their variation patterns and lead to crop improvements. In this study, the weedy rice line LM8 was examined at the whole-genome level. To identify the Oryza sativa japonica subpopulation that LM8 belongs to, the single nucleotide polymorphisms (SNPs) of 180 cultivated and 23 weedy rice varieties were used to construct a phylogenetic tree and a principal component analysis and STRUCTURE analysis were performed. The results indicated that LM8 with admixture components from japonica (GJ) and indica (XI) belonged to GJ-admixture (GJ-adm), with more than 60% of its genetic background derived from XI-2 (22.98%), GJ-tropical (22.86%), and GJ-subtropical (17.76%). Less than 9% of its genetic background was introgressed from weedy rice. Our results also suggested LM8 may have originated in a subtropical or tropical geographic region. Moreover, the comparisons with Nipponbare (NIP) and Shuhui498 (R498) revealed many specific structure variations (SVs) in the LM8 genome and fewer SVs between LM8 and NIP than between LM8 and R498. Next, 96 WRKY and 464 NBS genes were identified and mapped on LM8 chromosomes to eliminate redundancies. Three WRKY genes (ORUFILM02g002693, ORUFILM05g002725, and ORUFILM05g001757) in group III and one RNL [including the resistance to powdery mildew 8 (RPW8) domain, NBS, and leucine rich repeats (LRRs)] type NBS gene (ORUFILM12g000772) were detected in LM8. Among the NBS genes, the RPW8 domain was detected only in ORUFILM12g000772. This gene may improve plant resistance to pathogens as previously reported. Its classification and potential utility imply LM8 should be considered as a germplasm resource relevant for rice breeding programs.

High-Quality Genomes and High-Density Genetic Map Facilitate the Identification of Genes From a Weedy Rice.

Genes have been lost or weakened from cultivated rice during rice domestication and breeding. Weedy rice (Oryza sativa f. spontanea) is usually recognized as the progeny between cultivated rice and wild rice and is also known to harbor an gene pool for rice breeding. Therefore, identifying genes from weedy rice germplasms is an important way to break the bottleneck of rice breeding. To discover genes from weedy rice germplasms, we constructed a genetic map based on w-hole-genome sequencing of a F2 population derived from the cross between LM8 and a cultivated rice variety. We further identified 31 QTLs associated with 12 important agronomic traits and revealed that ORUFILM03g000095 gene may play an important role in grain length regulation and participate in grain formation. To clarify the genomic characteristics from weedy rice germplasms of LM8, we generated a high-quality genome assembly using single-molecule sequencing, Bionano optical mapping, and Hi-C technologies. The genome harbored a total size of 375.8 Mb, a scaffold N50 of 24.1 Mb, and originated approximately 0.32 million years ago (Mya) and was more closely related to Oryza sativa ssp. japonica. and contained 672 unique genes. It is related to the formation of grain shape, heading date and tillering. This study generated a high-quality reference genome of weedy rice and high-density genetic map that would benefit the analysis of genome evolution for related species and suggested an effective way to identify genes related to important agronomic traits for further rice breeding.