MORE FLORET1 Encodes a MYB Transcription Factor That Regulates Spikelet Development in Rice

MORE FLORET1 Interacts with C-type Replication Protein A Complex and Regulates Male Meiosis in Rice

Meiosis plays a pivotal role in plant reproduction, which is also crucial for enhancing genetic diversity. Although the impact of MOF1 on floral organ development and its negative regulation of the key tapetal gene PKS2 have been established, the specific function of MOF1 in male meiotic process remains elusive. In this study, we identified two mutant lines of MOF1 in Nipponbare background. Compared to the wild-type controls, MOF1 mutations resulted in significant reductions in seed setting rate and pollen fertility, partially attributed to its defects in the formation of male meiotic bivalents. RNA-seq analyses and RT-qPCR assays revealed that loss-of-function mutation of MOF1 didn't alter expression levels of 60 known meiotic-regulated genes, suggesting that MOF1 may not function as a transcriptional factor in its meiotic regulation. Yeast two-hybrid and bimolecular fluorescence complementation assays demonstrated the protein-protein interactions among MOF1, RPA2c, RPA1c, as well as FAR1, among which RPA1c and RPA2c involved in meiotic bivalent formation. Furthermore, gene expression pattern analyses and subcellular localization studies indicated the co-expression among above interacted proteins in nucleus during anther development. Our findings provide a mechanistic insight into how MOF1 modulate male meiosis possibly through interactions with key meiotic proteins, facilitating a better understanding of male reproductive regulation.

Chromosomal inversion at the DG1 promoter drives double-grain spikelets and enhances grain yield in sorghum

The phenomenon of multiple-grain spikelets is frequently observed in gramineous crops. In the case of dual-floret spikelets, the upper fertile floret develops normally to form a single grain, while the lower sterile floret undergoes abortion. Here we elucidate the role of Double-Grain 1 (DG1), a gene encoding a homeobox-domain-containing protein, in regulating the lower floret meristem activity and double-grain spikelet trait in sorghum. A 35.7-kb paracentric inversion in the DG1 promoter region leads to increased DG1 expression, probably by reducing repressive histone modifications. This increase in DG1 expression transforms the degenerated lower floret into a fertile one. The use of the superior DG1 allele results in an increase of approximately 40.7% to 46.1% in grain number per panicle and a 10.1% to 14.3% increase in overall grain yield. Our findings shed light on the sorghum double-grain spikelet characteristic, offering valuable insights for high-yield breeding designs in cereals.

MORE FLORET1 controls anther development by negatively regulating key tapetal genes in both diploid and tetraploid rice

Polyploid hybrid rice (Oryza sativa) has great potential for increasing yields. However, hybrid rice depends on male fertility and its regulation, which is less well studied in polyploid rice than in diploid rice. We previously identified an MYB transcription factor, MORE FLORET1 (MOF1), whose mutation causes male sterility in neo-tetraploid rice. MOF1 expression in anthers peaks at anther Stage 7 (S7) and progressively decreases to low levels at S10. However, it remains unclear how the dynamics of MOF1 expression contribute to male fertility. Here, we carefully examined anther development in both diploid and tetraploid mof1 rice mutants, as well as lines ectopically expressing MOF1 in a temporal manner. MOF1 mutations caused delayed degeneration of the tapetum and middle layer of anthers and aberrant pollen wall organization. Ectopic MOF1 expression at later stages of anther development led to retarded cytoplasmic reorganization of tapetal cells. In both cases, pollen grains were aborted and seed production was abolished, indicating that precise control of MOF1 expression is essential for male reproduction. We demonstrated that 5 key tapetal genes, CYP703A3 (CYTOCHROME P450 HYDROXYLASE 703A3), OsABCG26 (O. sativa ATP BINDING CASSETTE G26), PTC1 (PERSISTENT TAPETAL CELL1), PKS2 (POLYKETIDE SYNTHASE 2), and OsABCG15 (O. sativa ATP BINDING CASSETTE G15), exhibit expression patterns opposite to those of MOF1 and are negatively regulated by MOF1. Moreover, DNA affinity purification sequencing (DAP-seq), luciferase activity assays, and electrophoretic mobility shift assays indicated that MOF1 binds directly to the PKS2 promoter for transcriptional repression. Our results provide a mechanistic basis for the regulation of male reproduction by MOF1 in both diploid and tetraploid rice. This study will facilitate the development of polyploid male sterile lines, which are useful for breeding of polyploid hybrid rice.

The effect of exogenous melatonin on waterlogging stress in Clematis

Clematis is the queen of the vines, being an ornamental plant with high economic value. Waterlogging stress reduces the ornamental value of the plant and limits its application. Melatonin plays an important role in plant resistance to abiotic stresses. In this study, the physiological responses and gene expression levels of two wild species, namely, Clematis tientaiensis and Clematis lanuginosa, and two horticultural varieties, namely, 'Sen-No-Kaze' and 'Viva Polonia,' under waterlogging stress were analyzed to determine the effect of melatonin on waterlogging tolerance. The results showed that the waterlogging tolerances of C. lanuginosa and 'Sen-No-Kaze' were relatively poor, but were significantly improved by concentrations of 100 μmol·L-1 and 50 μmol·L-1 melatonin. C. tientaiensis and 'Viva Polonia' had relatively strong tolerance to waterlogging, and this was significantly improved by 200 μmol·L-1 melatonin. Under waterlogging stress, the relative conductivity and H2O2 content of Clematis increased significantly; the photosynthetic parameters and chlorophyll contents were significantly decreased; photosynthesis was inhibited; the contents of soluble protein and soluble sugars were decreased. Effective improvement of waterlogging tolerance after exogenous melatonin spraying, the relative conductivity was decreased by 4.05%-27.44%; the H2O2 content was decreased by 3.84%-23.28%; the chlorophyll content was increased by 35.59%-103.36%; the photosynthetic efficiency was increased by 25.42%-45.86%; the antioxidant enzyme activities of APX, POD, SOD, and CAT were increased by 28.03%-158.61%; the contents of proline, soluble protein, and soluble sugars were enhanced, and cell homeostasis was improved. Transcription sequencing was performed on wild Clematis with differences in waterlogging tolerance, and nine transcription factors were selected that were highly correlated with melatonin and that had the potential to improve waterlogging tolerance, among which LBD4, and MYB4 were significantly positively correlated with the antioxidant enzyme system, and bHLH36, DOF36, and WRKY4 were significantly negatively correlated. Photosynthetic capacity was positively correlated with DOF36 and WRKY4 while being significantly negatively correlated with MYB4, MOF1, DOF47, REV1 and ABR1. Melatonin could enhance the flooding tolerance of Clematis by improving photosynthetic efficiency and antioxidant enzyme activity. This study provides an important basis and reference for the application of melatonin in waterlogging-resistant breeding of Clematis.

OsMADS6-OsMADS32 and REP1 control palea cellular heterogeneity and morphogenesis in rice

Precise regulation of cell proliferation and differentiation is vital for organ morphology. Rice palea, serving as sepal, comprises two distinct regions: the marginal region (MRP) and body of palea (BOP), housing heterogeneous cell populations, which makes it an ideal system for studying organ morphogenesis. We report that the transcription factor (TF) REP1 promotes epidermal cell proliferation and differentiation in the BOP, resulting in hard silicified protrusion cells, by regulating the cyclin-dependent kinase gene, OsCDKB1;1. Conversely, TFs OsMADS6 and OsMADS32 are expressed exclusively in the MRP, where they limit cell division rates by inhibiting OsCDKB2;1 expression and promote endoreduplication, yielding elongated epidermal cells. Furthermore, reciprocal inhibition between the OsMADS6-OsMADS32 complex and REP1 fine-tunes the balance between cell division and differentiation during palea morphogenesis. We further show the functional conservation of these organ identity genes in heterogeneous cell growth in Arabidopsis, emphasizing a critical framework for controlling cellular heterogeneity in organ morphogenesis.

OsCRLK2, a Receptor-Like Kinase Identified by QTL Analysis, is Involved in the Regulation of Rice Quality

The quality of rice (Oryza sativa L) is determined by a combination of appearance, flavor, aroma, texture, storage characteristics, and nutritional composition. Rice quality directly influences acceptance by consumers and commercial value. The genetic mechanism underlying rice quality is highly complex, and is influenced by genotype, environment, and chemical factors such as starch type, protein content, and amino acid composition. Minor variations in these chemical components may lead to substantial differences in rice quality. Among these components, starch is the most crucial and influential factor in determining rice quality. In this study, quantitative trait loci (QTLs) associated with eight physicochemical properties related to the rapid viscosity analysis (RVA) profile were identified using a high-density sequence map constructed using recombinant inbred lines (RILs). Fifty-nine QTLs were identified across three environments, among which qGT6.4 was a novel locus co-located across all three environments. By integrating RNA-seq data, we identified the differentially expressed candidate gene OsCRLK2 within the qGT6.4 interval. osclrk2 mutants exhibited decreased gelatinization temperature (GT), apparent amylose content (AAC) and viscosity, and increased chalkiness. Furthermore, osclrk2 mutants exhibited downregulated expression of the majority of starch biosynthesis-related genes compared to wild type (WT) plants. In summary, OsCRLK2, which encodes a receptor-like protein kinase, appears to consistently influence rice quality across different environments. This discovery provides a new genetic resource for use in the molecular breeding of rice cultivars with improved quality.

Identification of Genomic Regions for Traits Associated with Flowering in Cassava (Manihot esculenta Crantz)

Flowering in cassava (Manihot esculenta Crantz) is crucial for the generation of botanical seed for breeding. However, genotypes preferred by most farmers are erect and poor at flowering or never flower. To elucidate the genetic basis of flowering, 293 diverse cassava accessions were evaluated for flowering-associated traits at two locations and seasons in Uganda. Genotyping using the Diversity Array Technology Pty Ltd. (DArTseq) platform identified 24,040 single-nucleotide polymorphisms (SNPs) distributed on the 18 cassava chromosomes. Population structure analysis using principal components (PCs) and kinships showed three clusters; the first five PCs accounted for 49.2% of the observed genetic variation. Linkage disequilibrium (LD) estimation averaged 0.32 at a distance of ~2850 kb (kilo base pairs). Polymorphism information content (PIC) and minor allele frequency (MAF) were 0.25 and 0.23, respectively. A genome-wide association study (GWAS) analysis uncovered 53 significant marker-trait associations (MTAs) with flowering-associated traits involving 27 loci. Two loci, SNPs S5_29309724 and S15_11747301, were associated with all the traits. Using five of the 27 SNPs with a Phenotype_Variance_Explained (PVE) ≥ 5%, 44 candidate genes were identified in the peak SNP sites located within 50 kb upstream or downstream, with most associated with branching traits. Eight of the genes, orthologous to Arabidopsis and other plant species, had known functional annotations related to flowering, e.g., eukaryotic translation initiation factor and myb family transcription factor. This study identified genomic regions associated with flowering-associated traits in cassava, and the identified SNPs can be useful in marker-assisted selection to overcome hybridization challenges, like unsynchronized flowering, and candidate gene validation.

Boosting Triticeae crop grain yield by manipulating molecular modules to regulate inflorescence architecture: insights and knowledge from other cereal crops

One of the challenges for global food security is to reliably and sustainably improve the grain yield of cereal crops. One solution is to modify the architecture of the grain-bearing inflorescence to optimize for grain number and size. Cereal inflorescences are complex structures, with determinacy, branching patterns, and spikelet/floret growth patterns that vary by species. Recent decades have witnessed rapid advancements in our understanding of the genetic regulation of inflorescence architecture in rice, maize, wheat, and barley. Here, we summarize current knowledge on key genetic factors underlying the different inflorescence morphologies of these crops and model plants (Arabidopsis and tomato), focusing particularly on the regulation of inflorescence meristem determinacy and spikelet meristem identity and determinacy. We also discuss strategies to identify and utilize these superior alleles to optimize inflorescence architecture and, ultimately, improve crop grain yield.

Evolution of cereal floral architecture and threshability

Hulled grains, while providing natural protection for seeds, pose a challenge to manual threshing due to the pair of glumes tightly encasing them. Based on natural evolution and artificial domestication, gramineous crops evolved various hull-like floral organs. Recently, progress has been made in uncovering novel domesticated genes associated with cereal threshability and deciphering common regulatory modules pertinent to the specification of hull-like floral organs. Here we review morphological similarities, principal regulators, and common mechanisms implicated in the easy-threshing traits of crops. Understanding the shared and unique features in the developmental process of cereal threshability may not only shed light on the convergent evolution of cereals but also facilitate the de novo domestication of wild cereal germplasm resources through genome-editing technologies.

Comparative cytological and transcriptome analyses of ny2 mutant delayed degeneration of tapetal cells and promotes abnormal microspore development in neo-tetraploid rice

We aimed to investigate the genetic defects related to pollen development and infertility in NY2, a novel tetraploid rice germplasm known as Neo-tetraploid rice. This rice variety was created through the crossbreeding and selective breeding of various autotetraploid rice lines and has previously shown high fertility. Our previous research has revealed that the NY2 gene, encoding a eukaryotic translation initiation factor 3 subunit E, regulates pollen fertility. However, the underlying mechanism behind this fertility is yet to be understood. To shed light on this matter, we performed a combined cytological and transcriptome analysis of the NY2 gene. Cytological analysis indicated that ny2 underwent abnormal tapetal cells, microspore, and middle layer development, which led to pollen abortion and ultimately to male sterility. Genetic analysis revealed that the F₁ plants showed normal fertility and an obvious advantage for seed setting compared to ny2. Global gene expression analysis in ny2 revealed a total of 7545 genes were detected at the meiosis stage, and 3925 and 3620 displayed upregulation and downregulation, respectively. The genes were significantly enriched for the gene ontology (GO) term "carbohydrate metabolic process. Moreover, 9 genes related to tapetum or pollen fertility showed down-regulation, such as OsABCG26 (ATP Binding Cassette G26), TMS9-1 (Thermosensitive Male Sterility), EAT1 (Programmed cell death regulatory), KIN14M (Kinesin Motor), OsMT1a (Metallothionein), and OsSTRL2 (Atypical strictosidine synthase), which were validated by qRT-PCR. Further analyses of DEGs identified nine down-regulated transcription factor genes related to pollen development. NY2 is an important regulator of the development of tapetum and microspore. The regulatory gene network described in this study may offer important understandings into the molecular processes that underlie fertility control in tetraploid rice.

A novel transcriptional repressor complex MYB22-TOPLESS-HDAC1 promotes rice resistance to brown planthopper by repressing F3'H expression

The brown planthopper (BPH) is the most destructive pest of rice. The MYB transcription factors are vital for rice immunity, but most are activators. Although MYB22 positively regulates rice resistance to BPH and has an EAR motif associated with active repression, it remains unclear whether it is a transcriptional repressor affecting rice-BPH interaction. Genetic analyses revealed that MYB22 regulates rice resistance to BPH via its EAR motif. Several biochemical experiments (e.g. transient transcription assay, Y2H, LCA, and BiFC) indicated that MYB22 is a transcriptional repressor that interacts with the corepressor TOPLESS via its EAR motif and recruits HDAC1 to form a tripartite complex. Flavonoid-3'-hydroxylase (F3'H) is a flavonoid biosynthesis pathway-related gene that negatively regulates rice resistance to BPH. Based on a bioinformatics analysis and the results of EMSA and transient transcription assays, MYB22 can bind directly to the F3'H promoter and repress gene expression along with TOPLESS and HDAC1. We revealed a transcriptional regulatory mechanism influencing the rice-BPH interaction that differs from previously reported mechanisms. Specifically, MYB22-TOPLESS-HDAC1 is a novel transcriptional repressor complex with components that synergistically and positively regulate rice resistance to BPH through the transcriptional repression of F3'H.

SiMYBS3, Encoding a Setaria italica Heterosis-Related MYB Transcription Factor, Confers Drought Tolerance in Arabidopsis

Drought is a major limiting factor affecting grain production. Drought-tolerant crop varieties are required to ensure future grain production. Here, 5597 DEGs were identified using transcriptome data before and after drought stress in foxtail millet (Setaria italica) hybrid Zhangza 19 and its parents. A total of 607 drought-tolerant genes were screened through WGCNA, and 286 heterotic genes were screened according to the expression level. Among them, 18 genes overlapped. One gene, Seita.9G321800, encoded MYBS3 transcription factor and showed upregulated expression after drought stress. It is highly homologous with MYBS3 in maize, rice, and sorghum and was named SiMYBS3. Subcellular localization analysis showed that the SiMYBS3 protein was located in the nucleus and cytoplasm, and transactivation assay showed SiMYBS3 had transcriptional activation activity in yeast cells. Overexpression of SiMYBS3 in Arabidopsis thaliana conferred drought tolerance, insensitivity to ABA, and earlier flowering. Our results demonstrate that SiMYBS3 is a drought-related heterotic gene and it can be used for enhancing drought resistance in agricultural crop breeding.

The upstream regulatory mechanism of BplMYB46 and the function of upstream regulatory factors that mediate resistance to stress in Betula platyphylla

Previously, we have shown that the transcription factor BplMYB46 in Betula platyphylla can enhance tolerance to salt and osmotic stress and promote secondary cell wall deposition, and we characterized its downstream regulatory mechanism. However, its upstream regulatory mechanism remains unclear. Here, the promoter activity and upstream regulatory factors of BplMYB46 were studied. Analyses of β-glucuronidase (GUS) staining and activity indicated that BplMYB46 promoter was specific temporal and spatial expression, and its expression can be induced by salt and osmotic stress. We identified three upstream regulatory factors of BplMYB46: BpDof1, BpWRKY3, and BpbZIP3. Yeast-one hybrid assays, GUS activity, chromatin immunoprecipitation, and quantitative real-time polymerase chain reaction revealed that BpDof1, BpWRKY3, and BpbZIP3 can directly regulate the expression of BplMYB46 by specifically binding to Dof, W-box, and ABRE elements in the BplMYB46 promoter, respectively. BpDof1, BpWRKY3, and BpbZIP3 were all localized to the nucleus, and their expressions can be induced by stress. Overexpression of BpDof1, BpWRKY3, and BpbZIP3 conferred the resistance of transgenic birch plants to salt and osmotic stress. Our findings provide new insights into the upstream regulatory mechanism of BplMYB46 and reveal new upstream regulatory genes that mediate resistance to adverse environments. The genes identified in our study provide novel targets for the breeding of forest tree species.

EgMYB108 regulates very long-chain fatty acid (VLCFA) anabolism in the mesocarp of oil palm

EgMYB108 regulates VLCFA anabolism in oil palm. Very long-chain fatty acids (VLCFAs), which are fatty acids with more than 18 C, can not only be used as a form of triglyceride (TAG) but also provide precursors for the biosynthesis of cuticle wax, and they exist in plant epidermal cells in the form of wax in higher plants. However, which and how transcriptional factors (TFs) regulate this process is largely unknown in oil palm. In this study, a MYB transcription factor (EgMYB108) with high expression in the mesocarp of oil palm fruit was characterized. Overexpression of EgMYB108 promoted not only total lipid content but also VLCFA accumulation in oil palm embryoids. Subsequently, transient transformation in protoplasts and qRT-PCR analysis indicated that the EgKCS5 and EgLACS4 genes were significantly increased with the overexpression of EgMYB108. Furthermore, yeast one‑hybrid assays, dual-luciferase assays and EMSAs demonstrated that EgMYB108 binds to the promoters of EgKCS5 and EgLACS4 and regulates their transcription. Finally, EgMYB108 interacts with the promoters of EgLACS and EgKCS simultaneously and finally improves the VLCFA and total lipid contents; a pathway summarizing this interaction was depicted.. The results provide new insight into the mechanism by which EgMYB108 regulates lipid and VLCFA accumulation in oil palm.

Identification and Analysis of MYB Gene Family for Discovering Potential Regulators Responding to Abiotic Stresses in Curcuma wenyujin

MYB superfamily is one of the most abundant families in plants, and plays critical role in plant growth, development, metabolism regulation, and stress response. Curcuma wenyujin is the main source plant of three traditional Chinese medicines, which are widely used in clinical treatment due to its diverse pharmacological activities. In present study, 88 CwMYBs were identified and analyzed in C. wenyujin, including 43 MYB-related genes, 42 R2R3-MYB genes, two 3R-MYB genes, and one 4R-MYB gene. Forty-three MYB-related proteins were classified into several types based on conserved domains and specific motifs, including CCA1-like type, R-R type, Myb-CC type, GARP-like type, and TBR-like type. The analysis of motifs in MYB DBD and no-MYB regions revealed the relevance of protein structure and function. Comparative phylogeny analysis divided 42 R2R3-MYB proteins into 19 subgroups and provided a reference for understanding the functions of some CwMYBs based on orthologs of previously characterized MYBs. Expression profile analysis of CwMYB genes revealed the differentially expressed genes responding to various abiotic stresses. Four candidate MYB genes were identified by combining the results of phylogeny analysis and expression analysis. CwMYB10, CwMYB18, CwMYB39, and CwMYB41 were significantly induced by cold, NaCl, and MeJA stress treatments. CwMYB18 and CwMYB41 were proved as regulators with activity of transcriptional activation, whereas CwMYB39 and CwMYB10 were not. They may participate in the response to abiotic stresses through different mechanisms in C. wenyujin. This study was the first step toward understanding the CwMYB family and the response to abiotic stresses in C. wenyujin.

What are the regulatory targets for intervention in assimilate partitioning to improve crop yield and resilience?

Sucrose utilisation for the synthesis of cellular components involved in growth and development and the accumulation of biomass determines diversity in the plant kingdom; sucrose utilisation and partitioning also underpin crop yields. As a complex process the use of sucrose for the partitioning of plant products for yield is decided by the interaction of several regulatory hubs and the integration of metabolism and development. Understanding the regulation of assimilate partitioning has been a grand challenge in plant and crop science. There are emerging examples of genes and processes that appear important for assimilate partitioning that underpin yield in crops and which are amenable to intervention. Enzymes of carbon metabolism were some of the first targets in attempts to modify assimilate partitioning at the beginning (source) and end (sink) of the whole plant assimilate partitioning process. Metabolic enzymes are subject to regulatory and homeostatic mechanisms, a key factor to consider in modifying assimilate partitioning. Trehalose 6-phosphate, as a sucrose signal, may represent a special case in its ability to regulate and coordinate source and sink processes. This review summarises recent progress in understanding the underlying regulators of assimilate partitioning and the current and potentially most promising routes to crop yield enhancement with a main focus on cereals. A framework for how source-sink may regulate whole plant assimilate partitioning involving a few key elements and the central importance of reproductive development is presented.

Modification of cereal plant architecture by genome editing to improve yields

KEY MESSAGE

We summarize recent genome editing studies that have focused on the examination (or reexamination) of plant architectural phenotypes in cereals and the modification of these traits for crop improvement. Plant architecture is defined as the three-dimensional organization of the entire plant. Shoot architecture refers to the structure and organization of the aboveground components of a plant, reflecting the developmental patterning of stems, branches, leaves and inflorescences/flowers. Root system architecture is essentially determined by four major shape parameters-growth, branching, surface area and angle. Interest in plant architecture has arisen from the profound impact of many architectural traits on agronomic performance, and the genetic and hormonal regulation of these traits which makes them sensitive to both selective breeding and agronomic practices. This is particularly important in staple crops, and a large body of literature has, therefore, accumulated on the control of architectural phenotypes in cereals, particularly rice due to its twin role as one of the world's most important food crops as well as a model organism in plant biology and biotechnology. These studies have revealed many of the molecular mechanisms involved in the regulation of tiller/axillary branching, stem height, leaf and flower development, root architecture and the grain characteristics that ultimately help to determine yield. The advent of genome editing has made it possible, for the first time, to introduce precise mutations into cereal crops to optimize their architecture and close in on the concept of the ideotype. In this review, we consider recent genome editing studies that have focused on the examination (or reexamination) of plant architectural phenotypes in cereals and the modification of these traits for crop improvement.

Transcriptome and Gene Editing Analyses Reveal MOF1a Defect Alters the Expression of Genes Associated with Tapetum Development and Chromosome Behavior at Meiosis Stage Resulting in Low Pollen Fertility of Tetraploid Rice

Autotetraploid rice is a useful rice germplasm for polyploid rice breeding. However, low fertility limits its commercial production. A neo-tetraploid rice with high fertility was developed from the progenies of crossing between autotetraploid lines by our research group. Our previous study showed that a myeloblastosis (MYB) transcription factor, MOF1, might be associated with the pollen development in tetraploid rice. However, little information is available about its role in pollen development in tetraploid rice. Here, we identified a new haplotype of MOF1 from neo-tetraploid rice and marked it as MOF1a. Transcriptome and qRT-PCR analysis demonstrated that MOF1a highly expressed in anthers, and displayed differential expression in neo-tetraploid rice compared to tetraploid rice line with low pollen fertility. The mutant (mof1a) of MOF1a, which was generated by CRISPR/Cas9 knockout in neo-tetraploid rice, showed low pollen fertility, and also exhibited abnormal tapetum and middle layer development, and defective chromosome behaviors during meiosis. A total of 13 tapetal related genes were found to be up-regulated in meiotic anthers of MOF1a compared with wild type plants by RNA-seq analysis, including CYP703A3, PTC1, and OsABCG26, which had been demonstrated to affect tapetal development. Moreover, 335 meiosis-related genes displayed differential expression patterns at same stage, including nine important meiosis-related genes, such as metallothionein OsMT1a. These results demonstrated that MOF1a plays an important role in pollen development and provides a foundation for understanding the molecular mechanism underlying MOF1a in reproduction of tetraploid rice.