Population-level exploration of alternative splicing and its unique role in controlling agronomic traits of rice.

Alternative splicing (AS) plays crucial roles in regulating various biological processes in plants. However, the genetic mechanisms underlying AS and its role in controlling important agronomic traits in rice (Oryza sativa) remain poorly understood. In this study, we explored AS in rice leaves and panicles using the rice minicore collection. Our analysis revealed a high level of transcript isoform diversity, with approximately one fifth of potential isoforms acting as major transcripts in both tissues. Regarding the genetic mechanism of AS, we found that the splicing of 833 genes in the leaf and 1,230 genes in the panicle was affected by cis-genetic variation. Twenty-one percent of these AS events could only be explained by large structural variations. Approximately 77.5% of genes with significant splicing quantitative trait loci (sGenes) exhibited tissue-specific regulation, and AS can cause 26.9% (leaf) and 23.6% (panicle) of sGenes to have altered, lost or gained functional domains. Additionally, through splicing-phenotype association analysis, we identified phosphate-starvation induced RING-type E3 ligase (OsPIE1; LOC_Os01g72480), whose splicing ratio was significantly associated with plant height. In summary, this study provides an understanding of AS in rice and its contribution to the regulation of important agronomic traits.

Gibberellin signaling regulates lignin biosynthesis to modulate rice seed shattering.

The elimination of seed shattering was a key step in rice (Oryza sativa) domestication. In this paper, we show that increasing the gibberellic acid (GA) content or response in the abscission region enhanced seed shattering in rice. We demonstrate that SLENDER RICE1 (SLR1), the key repressor of GA signaling, could physically interact with the rice seed shattering-related transcription factors quantitative trait locus of seed shattering on chromosome 1 (qSH1), O. sativa HOMEOBOX 15 (OSH15), and SUPERNUMERARY BRACT (SNB). Importantly, these physical interactions interfered with the direct binding of these three regulators to the lignin biosynthesis gene 4-COUMARATE: COENZYME A LIGASE 3 (4CL3), thereby derepressing its expression. Derepression of 4CL3 led to increased lignin deposition in the abscission region, causing reduced rice seed shattering. Importantly, we also show that modulating GA content could alter the degree of seed shattering to increase harvest efficiency. Our results reveal that the "Green Revolution" phytohormone GA is important for regulating rice seed shattering, and we provide an applicable breeding strategy for high-efficiency rice harvesting.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice.

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

A rice variation map derived from 10 548 rice accessions reveals the importance of rare variants.

Detailed knowledge of the genetic variations in diverse crop populations forms the basis for genetic crop improvement and gene functional studies. In the present study, we analyzed a large rice population with a total of 10 548 accessions to construct a rice super-population variation map (RSPVM), consisting of 54 378 986 single nucleotide polymorphisms, 11 119 947 insertion/deletion mutations and 184 736 presence/absence variations. Assessment of variation detection efficiency for different population sizes revealed a sharp increase of all types of variation as the population size increased and a gradual saturation of that after the population size reached 10 000. Variant frequency analysis indicated that ∼90% of the obtained variants were rare, and would therefore likely be difficult to detect in a relatively small population. Among the rare variants, only 2.7% were predicted to be deleterious. Population structure, genetic diversity and gene functional polymorphism of this large population were evaluated based on different subsets of RSPVM, demonstrating the great potential of RSPVM for use in downstream applications. Our study provides both a rich genetic basis for understanding natural rice variations and a powerful tool for exploiting great potential of rare variants in future rice research, including population genetics and functional genomics.

Programmable RNA N6-methyladenosine editing with CRISPR/dCas13a in plants.

N6-methyladenonsine (m6A) is the most prevalent internal modification of messenger RNA (mRNA) and plays critical roles in mRNA processing and metabolism. However, perturbation of individual m6A modification to reveal its function and the phenotypic effects is still lacking in plants. Here, we describe the construction and characterization of programmable m6A editing tools by fusing the m6A writers, the core catalytic domain of the MTA and MTB complex, and the AlkB homologue 5 (ALKBH5) eraser, to catalytically dead Cas13a (dCas13a) to edit individual m6A sites on mRNAs. We demonstrated that our m6A editors could efficiently and specifically deposit and remove m6A modifications on specific RNA transcripts in both Nicotiana benthamiana and Arabidopsis thaliana. Moreover, we found that targeting SHORT-ROOT (SHR) transcripts with a methylation editor could significantly increase its m6A levels with limited off-target effects and promote its degradation. This leads to a boost in plant growth with enlarged leaves and roots, increased plant height, plant biomass, and total grain weight in Arabidopsis. Collectively, these findings suggest that our programmable m6A editing tools can be applied to study the functions of individual m6A modifications in plants, and may also have potential applications for future crop improvement.

BTA2 regulates tiller angle and the shoot gravity response through controlling auxin content and distribution in rice.

Tiller angle is a key agricultural trait that establishes plant architecture, which in turn strongly affects grain yield by influencing planting density in rice. The shoot gravity response plays a crucial role in the regulation of tiller angle in rice, but the underlying molecular mechanism is largely unknown. Here, we report the identification of the BIG TILLER ANGLE2 (BTA2), which regulates tiller angle by controlling the shoot gravity response in rice. Loss-of-function mutation of BTA2 dramatically reduced auxin content and affected auxin distribution in rice shoot base, leading to impaired gravitropism and therefore a big tiller angle. BTA2 interacted with AUXIN RESPONSE FACTOR7 (ARF7) to modulate rice tiller angle through the gravity signaling pathway. The BTA2 protein was highly conserved during evolution. Sequence variation in the BTA2 promoter of indica cultivars harboring a less expressed BTA2 allele caused lower BTA2 expression in shoot base and thus wide tiller angle during rice domestication. Overexpression of BTA2 significantly increased grain yield in the elite rice cultivar Huanghuazhan under appropriate dense planting conditions. Our findings thus uncovered the BTA2-ARF7 module that regulates tiller angle by mediating the shoot gravity response. Our work offers a target for genetic manipulation of plant architecture and valuable information for crop improvement by producing the ideal plant type.

A centromere map based on super pan-genome highlights the structure and function of rice centromeres.

Rice (Oryza sativa) is a significant crop worldwide with a genome shaped by various evolutionary factors. Rice centromeres are crucial for chromosome segregation, and contain some unreported genes. Due to the diverse and complex centromere region, a comprehensive understanding of rice centromere structure and function at the population level is needed. We constructed a high-quality centromere map based on the rice super pan-genome consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. We showed that rice centromeres have diverse satellite repeat CentO, which vary across chromosomes and subpopulations, reflecting their distinct evolutionary patterns. We also revealed that long terminal repeats (LTRs), especially young Gypsy-type LTRs, are abundant in the peripheral CentO-enriched regions and drive rice centromere expansion and evolution. Furthermore, high-quality genome assembly and complete telomere-to-telomere (T2T) reference genome enable us to obtain more centromeric genome information despite mapping and cloning of centromere genes being challenging. We investigated the association between structural variations and gene expression in the rice centromere. A centromere gene, OsMAB, which positively regulates rice tiller number, was further confirmed by expression quantitative trait loci, haplotype analysis and clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 methods. By revealing the new insights into the evolutionary patterns and biological roles of rice centromeres, our finding will facilitate future research on centromere biology and crop improvement.

Identification of natural allelic variation in TTL1 controlling thermotolerance and grain size by a rice super pan-genome.

Continuously increasing global temperatures present great challenges to food security. Grain size, one of the critical components determining grain yield in rice (Oryza sativa L.), is a prime target for genetic breeding. Thus, there is an immediate need for genetic improvement in rice to maintain grain yield under heat stress. However, quantitative trait loci (QTLs) endowing heat stress tolerance and grain size in rice are extremely rare. Here, we identified a novel negative regulator with pleiotropic effects, Thermo-Tolerance and grain Length 1 (TTL1), from the super pan-genomic and transcriptomic data. Loss-of-function mutations in TTL1 enhanced heat tolerance, and caused an increase in grain size by coordinating cell expansion and proliferation. TTL1 was shown to function as a transcriptional regulator and localized to the nucleus and cell membrane. Furthermore, haplotype analysis showed that hapL and hapS of TTL1 were obviously correlated with variations of thermotolerance and grain size in a core collection of cultivars. Genome evolution analysis of available rice germplasms suggested that TTL1 was selected during domestication of the indica and japonica rice subspecies, but still had much breeding potential for increasing grain length and thermotolerance. These findings provide insights into TTL1 as a novel potential target for the development of high-yield and thermotolerant rice varieties.

Systematic Analysis of BELL Family Genes in Zizania latifolia and Functional Identification of ZlqSH1a/b in Rice Seed Shattering.

The loss of seed shattering is an important event in crop domestication, and elucidating the genetic mechanisms underlying seed shattering can help reduce yield loss during crop production. This study is the first to systematically identify and analyse the BELL family of transcription factor-encoding genes in Chinese wild rice (Zizania latifolia). ZlqSH1a (Zla04G033720) and ZlqSH1b (Zla02G027130) were identified as key candidate genes involved in seed shattering in Z. latifolia. These genes were involved in regulating the development of the abscission layer (AL) and were located in the nucleus of the cell. Over-expression of ZlqSH1a and ZlqSH1b resulted in a complete AL between the grain and pedicel and significantly enhanced seed shattering after grain maturation in rice. Transcriptome sequencing revealed that 172 genes were differentially expressed between the wild type (WT) and the two transgenic (ZlqSH1a and ZlqSH1b over-expressing) plants. Three of the differentially expressed genes related to seed shattering were validated using qRT-PCR analysis. These results indicate that ZlqSH1a and ZlqSH1b are involved in AL development in rice grains, thereby regulating seed shattering. Our results could facilitate the genetic improvement of seed-shattering behaviour in Z. latifolia and other cereal crops.

Natural alleles of a uridine 5'-diphospho-glucosyltransferase gene responsible for differential endosperm development between upland rice and paddy rice.

Traditional upland rice generally exhibits insufficient grains resulting from abnormal endosperm development compared to paddy rice. However, the underlying molecular mechanism of this trait is poorly understood. Here, we cloned the uridine 5'-diphospho (UDP)-glucosyltransferase gene EDR1 (Endosperm Development in Rice) responsible for differential endosperm development between upland rice and paddy rice by performing quantitative trait loci analysis and map-based cloning. EDR1 was highly expressed in developing seeds during grain filling. Natural variations in EDR1 significantly reduced the UDP-glucosyltransferase activity of EDR1YZN compared to EDR1YD1 , resulting in abnormal endosperm development in the near-isogenic line, accompanied by insufficient grains and changes in grain quality. By analyzing the distribution of the two alleles EDR1YD1 and EDR1YZN among diverse paddy rice and upland rice varieties, we discovered that EDR1 was conserved in upland rice, but segregated in paddy rice. Further analyses of grain chalkiness in the alleles of EDR1YD1 and EDR1YZN varieties indicated that rice varieties harboring EDR1YZN and EDR1YD1 preferentially showed high chalkiness, and low chalkiness, respectively. Taken together, these results suggest that the UDP-glucosyltransferase gene EDR1 is an important determinant controlling differential endosperm development between upland rice and paddy rice.

OsRLR4 binds to the OsAUX1 promoter to negatively regulate primary root development in rice.

Root architecture is one of the most important agronomic traits that determines rice crop yield. The primary root (PR) absorbs mineral nutrients and provides mechanical support; however, the molecular mechanisms of PR elongation remain unclear in rice. Here, the two loss-of-function T-DNA insertion mutants of root length regulator 4 (OsRLR4), osrlr4-1 and osrlr4-2 with longer PR, and three OsRLR4 overexpression lines, OE-OsRLR4-1/-2/-3 with shorter PR compared to the wild type/Hwayoung (WT/HY), were identified. OsRLR4 is one of five members of the PRAF subfamily of the regulator chromosome condensation 1 (RCC1) family. Phylogenetic analysis of OsRLR4 from wild and cultivated rice indicated that it is under selective sweeps, suggesting its potential role in domestication. OsRLR4 controls PR development by regulating auxin accumulation in the PR tip and thus the root apical meristem activity. A series of biochemical and genetic analyses demonstrated that OsRLR4 functions directly upstream of the auxin transporter OsAUX1. Moreover, OsRLR4 interacts with the TRITHORAX-like protein OsTrx1 to promote H3K4me3 deposition at the OsAUX1 promoter, thus altering its transcription level. This work provides insight into the cooperation of auxin and epigenetic modifications in regulating root architecture and provides a genetic resource for plant architecture breeding.

Comparative Transcriptomics of Lowland Rice Varieties Uncovers Novel Candidate Genes for Adaptive Iron Excess Tolerance.

Iron (Fe) toxicity is a major challenge for plant cultivation in acidic waterlogged soil environments, where lowland rice is a major staple food crop. Only few studies have addressed the molecular characterization of excess Fe tolerance in rice, and these highlight different mechanisms for Fe tolerance. Out of 16 lowland rice varieties, we identified a pair of contrasting lines, Fe-tolerant Lachit and -susceptible Hacha. The two lines differed in their physiological and morphological responses to excess Fe, including leaf growth, leaf rolling, reactive oxygen species generation and Fe and metal contents. These responses were likely due to genetic origin as they were mirrored by differential gene expression patterns, obtained through RNA sequencing, and corresponding gene ontology term enrichment in tolerant vs. susceptible lines. Thirty-five genes of the metal homeostasis category, mainly root expressed, showed differential transcriptomic profiles suggestive of an induced tolerance mechanism. Twenty-two out of these 35 metal homeostasis genes were present in selection sweep genomic regions, in breeding signatures, and/or differentiated during rice domestication. These findings suggest that Fe excess tolerance is an important trait in the domestication of lowland rice, and the identified genes may further serve to design the targeted Fe tolerance breeding of rice crops.

Loci and natural alleles for cadmium-mediated growth responses revealed by a genome wide association study and transcriptome analysis in rice.

BACKGROUND: Cadmium (Cd) is a toxic heavy metal that is harmful to the environment and human health. Cd pollution threatens the cultivation of rice (Oryza sativa L.) in many countries. Improving rice performance under Cd stress could potentially improve rice productivity. RESULTS: In this study, 9 growth traits of 188 different cultivated rice accessions under normal and Cd stress conditions were found to be highly variable during the seedling stage. Based on ~3.3 million single nucleotide polymorphisms (SNPs), 119 Cd-mediated growth response (CGR) quantitative trait loci (QTL) were identified by a genome-wide association study (GWAS), 55 of which have been validated by previously reported QTL and 64 were new CGR loci. Combined with the data from the GWAS, transcriptome analysis, gene annotations from the gene ontology (GO) Slim database, and annotations and functions of homologous genes, 148 CGR candidate genes were obtained. Additionally, several reported genes have been found to play certain roles in CGRs. Seven Cd-related cloned genes were found among the CGR genes. Natural elite haplotypes/alleles in these genes that increased Cd tolerance were identified by a haplotype analysis of a diverse mini core collection. More importantly, this study was the first to uncover the natural variations of 5 GST genes that play important roles in CGRs. CONCLUSION: The exploration of Cd-resistant rice germplasm resources and the identification of elite natural variations related to Cd-resistance will help improve the tolerance of current major rice varieties to Cd, as well as provide raw materials and new genes for breeding Cd-resistant varieties.

Linkage analysis, GWAS, transcriptome analysis to identify candidate genes for rice seedlings in response to high temperature stress.

BACKGROUND: Rice plants suffer from the rising temperature which is becoming more and more prominent. Mining heat-resistant genes and applying them to rice breeding is a feasible and effective way to solve the problem. RESULT: Three main biomass traits, including shoot length, dry weight, and fresh weight, changed after abnormally high-temperature treatment in the rice seedling stage of a recombinant inbred lines and the natural indica germplasm population. Based on a comparison of the results of linkage analysis and genome-wide association analysis, two loci with lengths of 57 kb and 69 kb in qDW7 and qFW6, respectively, were associated with the rice response to abnormally high temperatures at the seedling stage. Meanwhile, based on integrated transcriptome analysis, some genes are considered as important candidate genes. Combining with known genes and analysis of homologous genes, it was found that there are eight genes in candidate intervals that need to be focused on in subsequent research. CONCLUSIONS: The results indicated several relevant loci, which would help researchers to further discover beneficial heat-resistant genes that can be applied to rice heat-resistant breeding.

OsbHLH98 regulates leaf angle in rice through transcriptional repression of OsBUL1.

Leaf angle is an important agronomic trait in cereals that helps determine plant yield by affecting planting density. However, the regulation mechanism of leaf angle remained elusive. Here, we show that OsbHLH98, a rice bHLH transcription factor, negatively regulates leaf angle. osbhlh98 mutant leaves formed a larger leaf angle, whereas transgenic plants overexpressing OsbHLH98 exhibited a slight reduction in leaf angle. We determined that the changes in leaf angle resulted from increased number and size of parenchyma cells on the adaxial side of the lamina joint in osbhlh98 mutants. Experiments using reporter constructs showed that OsbHLH98 is expressed on the adaxial side of lamina joints, consistent with its proposed function in regulating leaf angle. Furthermore, we established by chromatin immunoprecipitation and CUT&RUN that OsBUL1 is a direct downstream target of OsbHLH98. Transactivation assays and reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis indicated that OsbHLH98 represses OsBUL1 transcription. Our results demonstrate that OsbHLH98 negatively regulates leaf angle by counteracting brassinosteroid-induced cell elongation via the repression of OsBUL1 transcription. The characterization of OsbHLH98 and its role in determining leaf angle will lay the foundation to develop the ideal plant architecture for adaptation to high planting density.

A genome-wide survey of copy number variations reveals an asymmetric evolution of duplicated genes in rice.

BACKGROUND: Copy number variations (CNVs) are an important type of structural variations in the genome that usually affect gene expression levels by gene dosage effect. Understanding CNVs as part of genome evolution may provide insights into the genetic basis of important agricultural traits and contribute to the crop breeding in the future. While available methods to detect CNVs utilizing next-generation sequencing technology have helped shed light on prevalence and effects of CNVs, the complexity of crop genomes poses a major challenge and requires development of additional tools. RESULTS: Here, we generated genomic and transcriptomic data of 93 rice (Oryza sativa L.) accessions and developed a comprehensive pipeline to call CNVs in this large-scale dataset. We analyzed the correlation between CNVs and gene expression levels and found that approximately 13% of the identified genes showed a significant correlation between their expression levels and copy numbers. Further analysis showed that about 36% of duplicate pairs were involved in pseudogenetic events while only 5% of them showed functional differentiation. Moreover, the offspring copy mainly contributed to the expression levels and seemed more likely to become a pseudogene, whereas the parent copy tended to maintain the function of ancestral gene. CONCLUSION: We provide a high-accuracy CNV dataset that will contribute to functional genomics studies and molecular breeding in rice. We also showed that gene dosage effect of CNVs in rice is not exponential or linear. Our work demonstrates that the evolution of duplicated genes is asymmetric in both expression levels and gene fates, shedding a new insight into the evolution of duplicated genes.

Natural variation in the promoter of TGW2 determines grain width and weight in rice.

Understanding the genetic basis of natural variation in grain size among diverse rice varieties can help breeders develop high-yielding rice cultivars. Here, we report the discovery of qTGW2, a new semidominant quantitative trait locus for grain width and weight. The corresponding gene, TGW2, encodes CELL NUMBER REGULATOR 1 (OsCNR1) localized to the plasma membrane. A single nucleotide polymorphism (SNP) variation 1818 bp upstream of TGW2 is responsible for its different expression, leading to alteration in grain width and weight by influencing cell proliferation and expansion in glumes. TGW2 interacts with KRP1, a regulator of cell cycle in plants, to negatively regulate grain width and weight. Genetic diversity analysis of TGW2 in 141 rice accessions revealed it as a breeding target in a selective sweep region. Our findings provide new insights into the genetic mechanism underlying grain morphology and grain weight, and uncover a promising gene for improving rice yield.

Natural variation in the promoter of OsHMA3 contributes to differential grain cadmium accumulation between Indica and Japonica rice.

Rice is a major source of cadmium (Cd) intake for Asian people. Indica rice usually accumulates more Cd in shoots and grains than Japonica rice. However, underlying genetic bases for differential Cd accumulation between Indica and Japonica rice are still unknown. In this study, we cloned a quantitative trait locus (QTL) grain Cd concentration on chromosome 7 (GCC7) responsible for differential grain Cd accumulation between two rice varieties by performing QTL analysis and map-based cloning. We found that the two GCC7 alleles, GCC7PA64s and GCC793-11 , had different promoter activity of OsHMA3, leading to different OsHMA3 expression and different shoot and grain Cd concentrations. By analyzing the distribution of different haplotypes of GCC7 among diverse rice accessions, we discovered that the high and low Cd accumulation alleles, namely GCC793-11 and GCC7PA64s , were preferentially distributed in Indica and Japonica rice, respectively. We further showed that the GCC7PA64s allele can be used to replace the GCC793-11 allele in the super cultivar 93-11 to reduce grain Cd concentration without adverse effect on agronomic traits. Our results thus reveal that the QTL GCC7 with sequence variation in the OsHMA3 promoter is an important determinant controlling differential grain Cd accumulation between Indica and Japonica rice.