Pan-genomes from large natural populations can capture genetic diversity and reveal genomic complexity. Using de novo long-read assembly, we generated a graph-based super pan-genome of rice consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. Our pan-genome reveals extensive structural variations (SVs) and gene presence/absence variations. Additionally, our pan-genome enables the accurate identification of nucleotide-binding leucine-rich repeat genes and characterization of their inter- and intraspecific diversity. Moreover, we uncovered grain weight-associated SVs which specify traits by affecting the expression of their nearby genes. We characterized genetic variants associated with submergence tolerance, seed shattering and plant architecture and found independent selection for a common set of genes that drove adaptation and domestication in Asian and African rice. This super pan-genome facilitates pinpointing of lineage-specific haplotypes for trait-associated genes and provides insights into the evolutionary events that have shaped the genomic architecture of various rice species.

Background: The advent of third-generation sequencing (TGS) technologies opens the door to improve genome assembly. Long reads are promising for enhancing the quality of fragmented draft assemblies constructed from next-generation sequencing (NGS) technologies. To date, a few algorithms that are capable of improving draft assemblies have released. There are SSPACE-LongRead, OPERA-LG, SMIS, npScarf, DBG2OLC, Unicycler, and LINKS. Hybrid assembly on large genomes remains challenging, however. Results: We develop a scalable and computationally efficient scaffolder, Long Reads Scaffolder (LRScaf, https://github.com/shingocat/lrscaf), that is capable of significantly boosting assembly contiguity using long reads. In this study, we summarise a comprehensive performance assessment for state-of-the-art scaffolders and LRScaf on seven organisms, i.e., E. coli, S. cerevisiae, A. thaliana, O. sativa, S. pennellii, Z. mays, and H. sapiens. LRScaf significantly improves the contiguity of draft assemblies, e.g., increasing the NGA50 value of CHM1 from 127.1 kbp to 9.4 Mbp using 20-fold coverage PacBio dataset and the NGA50 value of NA12878 from 115.3 kbp to 12.9 Mbp using 35-fold coverage Nanopore dataset. Besides, LRScaf generates the best contiguous NGA50 on A. thaliana, S. pennellii, Z. mays, and H. sapiens. Moreover, LRScaf has the shortest run time compared with other scaffolders, and the peak RAM of LRScaf remains practical for large genomes (e.g., 20.3 and 62.6 GB on CHM1 and NA12878, respectively). Conclusions: The new algorithm, LRScaf, yields the best or, at least, moderate scaffold contiguity and accuracy in the shortest run time compared with other scaffolding algorithms. Furthermore, LRScaf provides a cost-effective way to improve contiguity of draft assemblies on large genomes.

Plant-specific NAC proteins comprise one of the largest transcription factor families in plants and play important roles in plant development and the stress response. Gossypium hirsutum L. is a major source of fiber, but its growth and productivity are limited by many biotic and abiotic stresses. In this study, the NAC domain gene GhNAC79 was functionally characterized in detail, and according to information about the cotton genome sequences, it was located on scaffold42.1, containing three exons and two introns. Promoter analysis indicated that the GhNAC79 promoter contained both basic and stress-related elements, and it was especially expressed in the cotyledon of Arabidopsis. A transactivation assay in yeast demonstrated that GhNAC79 was a transcription activator, and its activation domain was located at its C-terminus. The results of qRT-PCR proved that GhNAC79 was preferentially expressed at later stages of cotyledon and fiber development, and it showed high sensitivity to ethylene and meJA treatments. Overexpression of GhNAC79 resulted in an early-flowering phenotype in Arabidopsis, and it also improved drought tolerance in both Arabidopsis and cotton. Furthermore, VIGS-induced silencing of GhNAC79 in cotton led to a drought-sensitive phenotype. In summary, GhNAC79 positively regulates drought stress, and it also responds to ethylene and meJA treatments, making it a candidate gene for stress studies in cotton.

Seed germination represents a major developmental switch in plants that is vital to agriculture, but how this process is controlled at the chromatin level remains obscure. Here we demonstrate that successful germination in Arabidopsis thaliana requires a chromatin mechanism that progressively silences 9-CIS-EPOXYCAROTENOID DIOXYGENASE 6 (NCED6), which encodes a rate-limiting enzyme in abscisic acid (ABA) biosynthesis, through the cooperative action of the RNA-binding protein RZ-1 and the polycomb repressive complex 2 (PRC2). Simultaneous inactivation of RZ-1 and PRC2 blocked germination and synergistically derepressed NCEDs and hundreds of genes. At NCED6, in part by promoting H3 deacetylation and suppressing H3K4me3, RZ-1 facilitates transcriptional silencing and also an H3K27me3 accumulation process that occurs during seed germination and early seedling growth. Genome-wide analysis revealed that RZ-1 is preferentially required for transcriptional silencing of many PRC2 targets early during seed germination, when H3K27me3 is not yet established. We propose RZ-1 confers a novel silencing mechanism to compensate for and synergize with PRC2. Our work highlights the progressive chromatin silencing of ABA biosynthesis genes via the RNA-binding protein RZ-1 and PRC2 acting in synergy, a process that is vital for seed germination.

Preharvest sprouting (PHS) is a deleterious phenotype that occurs frequently in rice-growing regions where the temperature and precipitation are high. It negatively affects yield, quality, and downstream grain processing. Seed dormancy is a trait related to PHS. Longer seed dormancy is preferred for rice production as it can prevent PHS. Here, we map QTLs associated with rice seed dormancy and clone Seed Dormancy 3.1 (SDR3.1) underlying one major QTL. SDR3.1 encodes a mediator of OsbZIP46 deactivation and degradation (MODD). We show that SDR3.1 negatively regulates seed dormancy by inhibiting the transcriptional activity of ABIs. In addition, we reveal two critical amino acids of SDR3.1 that are critical for the differences in seed dormancy between the Xian/indica and Geng/japonica cultivars. Further, SDR3.1 has been artificially selected during rice domestication. We propose a two-line model for the process of rice seed dormancy domestication from wild rice to modern cultivars. We believe the candidate gene and germplasm studied in this study would be beneficial for the genetic improvement of rice seed dormancy.