Dynamics of extrachromosomal circular DNA in rice.

The genome's dynamic nature, exemplified by elements like extrachromosomal circular DNA (eccDNA), is crucial for biodiversity and adaptation. Yet, the role of eccDNA in plants, particularly rice, remains underexplored. Here, we identify 25,598 eccDNAs, unveiling the widespread presence of eccDNA across six rice tissues and revealing its formation as a universal and random process. Interestingly, we discover that direct repeats play a pivotal role in eccDNA formation, pointing to a unique origin mechanism. Despite eccDNA's prevalence in coding sequences, its impact on gene expression is minimal, implying its roles beyond gene regulation. We also observe the association between eccDNA's formation and minor chromosomal deletions, providing insights of its possible function in regulating genome stability. Further, we discover eccDNA specifically accumulated in rice leaves, which may be associated with DNA damage caused by environmental stressors like intense light. In summary, our research advances understanding of eccDNA's role in the genomic architecture and offers valuable insights for rice cultivation and breeding.

Technical strategy for monozygotic twin discrimination by single-nucleotide variants.

Monozygotic (MZ) twins are theoretically genetically identical. Although they are revealed to accumulate mutations after the zygote splits, discriminating between twin genomes remains a formidable challenge in the field of forensic genetics. Single-nucleotide variants (SNVs) are responsible for a substantial portion of genetic variation, thus potentially serving as promising biomarkers for the identification of MZ twins. In this study, we sequenced the whole genome of a pair of female MZ twins when they were 27 and 33 years old to approximately 30 × coverage using peripheral blood on an Illumina NovaSeq 6000 Sequencing System. Potentially discordant SNVs supported by whole-genome sequencing were validated extensively by amplicon-based targeted deep sequencing and Sanger sequencing. In total, we found nine bona fide post-twinning SNVs, all of which were identified in the younger genomes and found in the older genomes. None of the SNVs occurred within coding exons, three of which were observed in introns, supported by whole-exome sequencing results. A double-blind test was employed, and the reliability of MZ twin discrimination by discordant SNVs was endorsed. All SNVs were successfully detected when input DNA amounts decreased to 0.25 ng, and reliable detection was limited to seven SNVs below 0.075 ng input. This comprehensive analysis confirms that SNVs could serve as cost-effective biomarkers for MZ twin discrimination.

Structure and function of rice hybrid genomes reveal genetic basis and optimal performance of heterosis.

Exploitation of crop heterosis is crucial for increasing global agriculture production. However, the quantitative genomic analysis of heterosis was lacking, and there is currently no effective prediction tool to optimize cross-combinations. Here 2,839 rice hybrid cultivars and 9,839 segregation individuals were resequenced and phenotyped. Our findings demonstrated that indica-indica hybrid-improving breeding was a process that broadened genetic resources, pyramided breeding-favorable alleles through combinatorial selection and collaboratively improved both parents by eliminating the inferior alleles at negative dominant loci. Furthermore, we revealed that widespread genetic complementarity contributed to indica-japonica intersubspecific heterosis in yield traits, with dominance effect loci making a greater contribution to phenotypic variance than overdominance effect loci. On the basis of the comprehensive dataset, a genomic model applicable to diverse rice varieties was developed and optimized to predict the performance of hybrid combinations. Our data offer a valuable resource for advancing the understanding and facilitating the utilization of heterosis in rice.

The WD40 domain-containing protein Ehd5 positively regulates flowering in rice (Oryza sativa).

Heading date (flowering time), which greatly influences regional and seasonal adaptability in rice (Oryza sativa), is regulated by many genes in different photoperiod pathways. Here, we characterized a heading date gene, Early heading date 5 (Ehd5), using a modified bulked segregant analysis method. The ehd5 mutant showed late flowering under both short-day and long-day conditions, as well as reduced yield, compared to the wild type. Ehd5, which encodes a WD40 domain-containing protein, is induced by light and follows a circadian rhythm expression pattern. Transcriptome analysis revealed that Ehd5 acts upstream of the flowering genes Early heading date 1 (Ehd1), RICE FLOWERING LOCUS T 1 (RFT1), and Heading date 3a (Hd3a). Functional analysis showed that Ehd5 directly interacts with Rice outermost cell-specific gene 4 (Roc4) and Grain number, plant height, and heading date 8 (Ghd8), which might affect the formation of Ghd7-Ghd8 complexes, resulting in increased expression of Ehd1, Hd3a, and RFT1. In a nutshell, these results demonstrate that Ehd5 functions as a positive regulator of rice flowering and provide insight into the molecular mechanisms underlying heading date.

Armadillo repeat only protein GS10 negatively regulates brassinosteroid signaling to control rice grain size.

Grain yield and grain quality are major determinants in modern breeding controlled by many quantitative traits loci (QTLs) in rice (Oryza sativa). However, the mechanisms underlying grain shape and quality are poorly understood. Here, we characterize a QTL for grain size and grain quality via map-based cloning from wild rice (W1943), GS10 (Grain Size on Chromosome 10), which encodes a protein with 6 tandem armadillo repeats. The null mutant gs10 shows slender and narrow grains with altered cell size, which has a pleiotropic effect on other agronomical traits. Functional analysis reveals that GS10 interacts with TUD1 (Taihu Dwarf1) and is epistatic to OsGSK2 (glycogen synthase kinase 2) through regulating grain shape and lamina joint inclination, indicating it is negatively involved in brassinosteroid (BR) signaling. Pyramiding gs10 and the grain size gene GW5 into cultivar GLA4 substantially improved grain shape and appearance quality. Natural variation analysis revealed that gs10 from the wild rice Oryza rufipogon W1943 is a rare allele across the rice population. Collectively, these findings advance our understanding of the underlying mechanism of grain shape and provide the beneficial allele of gs10 for future rice breeding and genetic improvement.

Chromosome-scale assembly and analysis of biomass crop Miscanthus lutarioriparius genome.

Miscanthus, a rhizomatous perennial plant, has great potential for bioenergy production for its high biomass and stress tolerance. We report a chromosome-scale assembly of Miscanthus lutarioriparius genome by combining Oxford Nanopore sequencing and Hi-C technologies. The 2.07-Gb assembly covers 96.64% of the genome, with contig N50 of 1.71 Mb. The centromere and telomere sequences are assembled for all 19 chromosomes and chromosome 10, respectively. Allotetraploid origin of the M. lutarioriparius is confirmed using centromeric satellite repeats. The tetraploid genome structure and several chromosomal rearrangements relative to sorghum are clearly demonstrated. Tandem duplicate genes of M. lutarioriparius are functional enriched not only in terms related to stress response, but cell wall biosynthesis. Gene families related to disease resistance, cell wall biosynthesis and metal ion transport are greatly expanded and evolved. The expansion of these families may be an important genomic basis for the enhancement of remarkable traits of M. lutarioriparius.

Genome-wide identification and analysis of heterotic loci in three maize hybrids.

Heterosis, or hybrid vigour, is a predominant phenomenon in plant genetics, serving as the basis of crop hybrid breeding, but the causative loci and genes underlying heterosis remain unclear in many crops. Here, we present a large-scale genetic analysis using 5360 offsprings from three elite maize hybrids, which identifies 628 loci underlying 19 yield-related traits with relatively high mapping resolutions. Heterotic pattern investigations of the 628 loci show that numerous loci, mostly with complete-incomplete dominance (the major one) or overdominance effects (the secondary one) for heterozygous genotypes and nearly equal proportion of advantageous alleles from both parental lines, are the major causes of strong heterosis in these hybrids. Follow-up studies for 17 heterotic loci in an independent experiment using 2225 F2 individuals suggest most heterotic effects are roughly stable between environments with a small variation. Candidate gene analysis for one major heterotic locus (ub3) in maize implies that there may exist some common genes contributing to crop heterosis. These results provide a community resource for genetics studies in maize and new implications for heterosis in plants.

A noncanonical vacuolar sugar transferase required for biosynthesis of antimicrobial defense compounds in oat.

Plants produce an array of natural products with important ecological functions. These compounds are often decorated with oligosaccharide groups that influence bioactivity, but the biosynthesis of such sugar chains is not well understood. Triterpene glycosides (saponins) are a large family of plant natural products that determine important agronomic traits, as exemplified by avenacins, antimicrobial defense compounds produced by oats. Avenacins have a branched trisaccharide moiety consisting of l-arabinose linked to 2 d-glucose molecules that is critical for antifungal activity. Plant natural product glycosylation is usually performed by uridine diphosphate-dependent glycosyltransferases (UGTs). We previously characterized the arabinosyltransferase that initiates the avenacin sugar chain; however, the enzymes that add the 2 remaining d-glucose molecules have remained elusive. Here we characterize the enzymes that catalyze these last 2 glucosylation steps. AsUGT91G16 is a classical cytosolic UGT that adds a 1,2-linked d-glucose molecule to l-arabinose. Unexpectedly, the enzyme that adds the final 1,4-linked d-glucose (AsTG1) is not a UGT, but rather a sugar transferase belonging to Glycosyl Hydrolase family 1 (GH1). Unlike classical UGTs, AsTG1 is vacuolar. Analysis of oat mutants reveals that AsTG1 corresponds to Sad3, a previously uncharacterized locus shown by mutation to be required for avenacin biosynthesis. AsTG1 and AsUGT91G16 form part of the avenacin biosynthetic gene cluster. Our demonstration that a vacuolar transglucosidase family member plays a critical role in triterpene biosynthesis highlights the importance of considering other classes of carbohydrate-active enzymes in addition to UGTs as candidates when elucidating pathways for the biosynthesis of glycosylated natural products in plants.

The PLATZ Transcription Factor GL6 Affects Grain Length and Number in Rice.

Grain size is one of the key determinants of grain yield. Although a number of genes that control grain size in rice (Oryza sativa) have been identified, the overall regulatory networks behind this process remain poorly understood. Here, we report the map-based cloning and functional characterization of the quantitative trait locus GL6, which encodes a plant-specific plant AT-rich sequence- and zinc-binding transcription factor that regulates rice grain length and spikelet number. GL6 positively controls grain length by promoting cell proliferation in young panicles and grains. The null gl6 mutant possesses short grains, whereas overexpression of GL6 results in large grains and decreased grain number per panicle. We demonstrate that GL6 participates in RNA polymerase III transcription machinery by interacting with RNA polymerase III subunit C53 and transcription factor class C1 to regulate the expression of genes involved in rice grain development. Our findings reveal a further player involved in the regulation of rice grain size that may be exploited in future rice breeding.

Publisher Correction: Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice.

When published, this article did not initially appear open access. This error has been corrected, and the open access status of the paper is noted in all versions of the paper.

Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice.

The rich genetic diversity in Oryza sativa and Oryza rufipogon serves as the main sources in rice breeding. Large-scale resequencing has been undertaken to discover allelic variants in rice, but much of the information for genetic variation is often lost by direct mapping of short sequence reads onto the O. sativa japonica Nipponbare reference genome. Here we constructed a pan-genome dataset of the O. sativa-O. rufipogon species complex through deep sequencing and de novo assembly of 66 divergent accessions. Intergenomic comparisons identified 23 million sequence variants in the rice genome. This catalog of sequence variations includes many known quantitative trait nucleotides and will be helpful in pinpointing new causal variants that underlie complex traits. In particular, we systemically investigated the whole set of coding genes using this pan-genome data, which revealed extensive presence and absence of variation among rice accessions. This pan-genome resource will further promote evolutionary and functional studies in rice.

Genomic architecture of heterosis for yield traits in rice.

Increasing grain yield is a long-term goal in crop breeding to meet the demand for global food security. Heterosis, when a hybrid shows higher performance for a trait than both parents, offers an important strategy for crop breeding. To examine the genetic basis of heterosis for yield in rice, here we generate, sequence and record the phenotypes of 10,074 F2 lines from 17 representative hybrid rice crosses. We classify modern hybrid rice varieties into three groups, representing different hybrid breeding systems. Although we do not find any heterosis-associated loci shared across all lines, within each group, a small number of genomic loci from female parents explain a large proportion of the yield advantage of hybrids over their male parents. For some of these loci, we find support for partial dominance of heterozygous locus for yield-related traits and better-parent heterosis for overall performance when all of the grain-yield traits are considered together. These results inform on the genomic architecture of heterosis and rice hybrid breeding.

OsSPL13 controls grain size in cultivated rice.

Although genetic diversity has a cardinal role in domestication, abundant natural allelic variations across the rice genome that cause agronomically important differences between diverse varieties have not been fully explored. Here we implement an approach integrating genome-wide association testing with functional analysis on grain size in a diverse rice population. We report that a major quantitative trait locus, GLW7, encoding the plant-specific transcription factor OsSPL13, positively regulates cell size in the grain hull, resulting in enhanced rice grain length and yield. We determine that a tandem-repeat sequence in the 5' UTR of OsSPL13 alters its expression by affecting transcription and translation and that high expression of OsSPL13 is associated with large grains in tropical japonica rice. Further analysis indicates that the large-grain allele of GLW7 in tropical japonica rice was introgressed from indica varieties under artificial selection. Our study demonstrates that new genes can be effectively identified on the basis of genome-wide association data.

Genetic discovery for oil production and quality in sesame.

Oilseed crops are used to produce vegetable oil. Sesame (Sesamum indicum), an oilseed crop grown worldwide, has high oil content and a small diploid genome, but the genetic basis of oil production and quality is unclear. Here we sequence 705 diverse sesame varieties to construct a haplotype map of the sesame genome and de novo assemble two representative varieties to identify sequence variations. We investigate 56 agronomic traits in four environments and identify 549 associated loci. Examination of the major loci identifies 46 candidate causative genes, including genes related to oil content, fatty acid biosynthesis and yield. Several of the candidate genes for oil content encode enzymes involved in oil metabolism. Two major genes associated with lignification and black pigmentation in the seed coat are also associated with large variation in oil content. These findings may inform breeding and improvement strategies for a broad range of oilseed crops.

Transcriptome-wide investigation of circular RNAs in rice.

Various stable circular RNAs (circRNAs) are newly identified to be the abundance of noncoding RNAs in Archaea, Caenorhabditis elegans, mice, and humans through high-throughput deep sequencing coupled with analysis of massive transcriptional data. CircRNAs play important roles in miRNA function and transcriptional controlling by acting as competing endogenous RNAs or positive regulators on their parent coding genes. However, little is known regarding circRNAs in plants. Here, we report 2354 rice circRNAs that were identified through deep sequencing and computational analysis of ssRNA-seq data. Among them, 1356 are exonic circRNAs. Some circRNAs exhibit tissue-specific expression. Rice circRNAs have a considerable number of isoforms, including alternative backsplicing and alternative splicing circularization patterns. Parental genes with multiple exons are preferentially circularized. Only 484 circRNAs have backsplices derived from known splice sites. In addition, only 92 circRNAs were found to be enriched for miniature inverted-repeat transposable elements (MITEs) in flanking sequences or to be complementary to at least 18-bp flanking intronic sequences, indicating that there are some other production mechanisms in addition to direct backsplicing in rice. Rice circRNAs have no significant enrichment for miRNA target sites. A transgenic study showed that overexpression of a circRNA construct could reduce the expression level of its parental gene in transgenic plants compared with empty-vector control plants. This suggested that circRNA and its linear form might act as a negative regulator of its parental gene. Overall, these analyses reveal the prevalence of circRNAs in rice and provide new biological insights into rice circRNAs.

The draft genome of the grass carp (Ctenopharyngodon idellus) provides insights into its evolution and vegetarian adaptation.

The grass carp is an important farmed fish, accounting for ∼16% of global freshwater aquaculture, and has a vegetarian diet. Here we report a 0.9-Gb draft genome of a gynogenetic female adult and a 1.07-Gb genome of a wild male adult. Genome annotation identified 27,263 protein-coding gene models in the female genome. A total of 114 scaffolds consisting of 573 Mb are anchored on 24 linkage groups. Divergence between grass carp and zebrafish is estimated to have occurred 49-54 million years ago. We identify a chromosome fusion in grass carp relative to zebrafish and report frequent crossovers between the grass carp X and Y chromosomes. We find that transcriptional activation of the mevalonate pathway and steroid biosynthesis in liver is associated with the grass carp's adaptation from a carnivorous to an herbivorous diet. We believe that the grass carp genome could serve as an initial platform for breeding better-quality fish using a genomic approach.

Genomic analysis of hybrid rice varieties reveals numerous superior alleles that contribute to heterosis.

Exploitation of heterosis is one of the most important applications of genetics in agriculture. However, the genetic mechanisms of heterosis are only partly understood, and a global view of heterosis from a representative number of hybrid combinations is lacking. Here we develop an integrated genomic approach to construct a genome map for 1,495 elite hybrid rice varieties and their inbred parental lines. We investigate 38 agronomic traits and identify 130 associated loci. In-depth analyses of the effects of heterozygous genotypes reveal that there are only a few loci with strong overdominance effects in hybrids, but a strong correlation is observed between the yield and the number of superior alleles. While most parental inbred lines have only a small number of superior alleles, high-yielding hybrid varieties have several. We conclude that the accumulation of numerous rare superior alleles with positive dominance is an important contributor to the heterotic phenomena.

A haplotype map of genomic variations and genome-wide association studies of agronomic traits in foxtail millet (Setaria italica).

Foxtail millet (Setaria italica) is an important grain crop that is grown in arid regions. Here we sequenced 916 diverse foxtail millet varieties, identified 2.58 million SNPs and used 0.8 million common SNPs to construct a haplotype map of the foxtail millet genome. We classified the foxtail millet varieties into two divergent groups that are strongly correlated with early and late flowering times. We phenotyped the 916 varieties under five different environments and identified 512 loci associated with 47 agronomic traits by genome-wide association studies. We performed a de novo assembly of deeply sequenced genomes of a Setaria viridis accession (the wild progenitor of S. italica) and an S. italica variety and identified complex interspecies and intraspecies variants. We also identified 36 selective sweeps that seem to have occurred during modern breeding. This study provides fundamental resources for genetics research and genetic improvement in foxtail millet.

The draft genome of the fast-growing non-timber forest species moso bamboo (Phyllostachys heterocycla).

Bamboo represents the only major lineage of grasses that is native to forests and is one of the most important non-timber forest products in the world. However, no species in the Bambusoideae subfamily has been sequenced. Here, we report a high-quality draft genome sequence of moso bamboo (P. heterocycla var. pubescens). The 2.05-Gb assembly covers 95% of the genomic region. Gene prediction modeling identified 31,987 genes, most of which are supported by cDNA and deep RNA sequencing data. Analyses of clustered gene families and gene collinearity show that bamboo underwent whole-genome duplication 7-12 million years ago. Identification of gene families that are key in cell wall biosynthesis suggests that the whole-genome duplication event generated more gene duplicates involved in bamboo shoot development. RNA sequencing analysis of bamboo flowering tissues suggests a potential connection between drought-responsive and flowering genes.