An epithelial cell-derived metabolite tunes immunoglobulin A secretion by gut-resident plasma cells.

Immunoglobulin A (IgA) secretion by plasma cells, terminally differentiated B cells residing in the intestinal lamina propria, assures microbiome homeostasis and protects the host against enteric infections. Exposure to diet-derived and commensal-derived signals provides immune cells with organizing cues that instruct their effector function and dynamically shape intestinal immune responses at the mucosal barrier. Recent data have described metabolic and microbial inputs controlling T cell and innate lymphoid cell activation in the gut; however, whether IgA-secreting lamina propria plasma cells are tuned by local stimuli is completely unknown. Although antibody secretion is considered to be imprinted during B cell differentiation and therefore largely unaffected by environmental changes, a rapid modulation of IgA levels in response to intestinal fluctuations might be beneficial to the host. In the present study, we showed that dietary cholesterol absorption and commensal recognition by duodenal intestinal epithelial cells lead to the production of oxysterols, evolutionarily conserved lipids with immunomodulatory functions. Using conditional cholesterol 25-hydroxylase deleter mouse line we demonstrated that 7α,25-dihydroxycholesterol from epithelial cells is critical to restrain IgA secretion against commensal- and pathogen-derived antigens in the gut. Intestinal plasma cells sense oxysterols via the chemoattractant receptor GPR183 and couple their tissue positioning with IgA secretion. Our findings revealed a new mechanism linking dietary cholesterol and humoral immune responses centered around plasma cell localization for efficient mucosal protection.

Rescue of Fragile X Syndrome Neurons by DNA Methylation Editing of the FMR1 Gene.

Fragile X syndrome (FXS), the most common genetic form of intellectual disability in males, is caused by silencing of the FMR1 gene associated with hypermethylation of the CGG expansion mutation in the 5' UTR of FMR1 in FXS patients. Here, we applied recently developed DNA methylation editing tools to reverse this hypermethylation event. Targeted demethylation of the CGG expansion by dCas9-Tet1/single guide RNA (sgRNA) switched the heterochromatin status of the upstream FMR1 promoter to an active chromatin state, restoring a persistent expression of FMR1 in FXS iPSCs. Neurons derived from methylation-edited FXS iPSCs rescued the electrophysiological abnormalities and restored a wild-type phenotype upon the mutant neurons. FMR1 expression in edited neurons was maintained in vivo after engrafting into the mouse brain. Finally, demethylation of the CGG repeats in post-mitotic FXS neurons also reactivated FMR1. Our data establish that demethylation of the CGG expansion is sufficient for FMR1 reactivation, suggesting potential therapeutic strategies for FXS.

High-throughput PRIME-editing screens identify functional DNA variants in the human genome.

Despite tremendous progress in detecting DNA variants associated with human disease, interpreting their functional impact in a high-throughput and single-base resolution manner remains challenging. Here, we develop a pooled prime-editing screen method, PRIME, that can be applied to characterize thousands of coding and non-coding variants in a single experiment with high reproducibility. To showcase its applications, we first identified essential nucleotides for a 716 bp MYC enhancer via PRIME-mediated single-base resolution analysis. Next, we applied PRIME to functionally characterize 1,304 genome-wide association study (GWAS)-identified non-coding variants associated with breast cancer and 3,699 variants from ClinVar. We discovered that 103 non-coding variants and 156 variants of uncertain significance are functional via affecting cell fitness. Collectively, we demonstrate that PRIME is capable of characterizing genetic variants at single-base resolution and scale, advancing accurate genome annotation for disease risk prediction, diagnosis, and therapeutic target identification.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project.

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

The chemokine CCL1 triggers an AMFR-SPRY1 pathway that promotes differentiation of lung fibroblasts into myofibroblasts and drives pulmonary fibrosis.

Recruitment of immune cells to the site of inflammation by the chemokine CCL1 is important in the pathology of inflammatory diseases. Here, we examined the role of CCL1 in pulmonary fibrosis (PF). Bronchoalveolar lavage fluid from PF mouse models contained high amounts of CCL1, as did lung biopsies from PF patients. Immunofluorescence analyses revealed that alveolar macrophages and CD4+ T cells were major producers of CCL1 and targeted deletion of Ccl1 in these cells blunted pathology. Deletion of the CCL1 receptor Ccr8 in fibroblasts limited migration, but not activation, in response to CCL1. Mass spectrometry analyses of CCL1 complexes identified AMFR as a CCL1 receptor, and deletion of Amfr impaired fibroblast activation. Mechanistically, CCL1 binding triggered ubiquitination of the ERK inhibitor Spry1 by AMFR, thus activating Ras-mediated profibrotic protein synthesis. Antibody blockade of CCL1 ameliorated PF pathology, supporting the therapeutic potential of targeting this pathway for treating fibroproliferative lung diseases.

NBS1 lactylation is required for efficient DNA repair and chemotherapy resistance.

The Warburg effect is a hallmark of cancer that refers to the preference of cancer cells to metabolize glucose anaerobically rather than aerobically1,2. This results in substantial accumulation of lacate, the end product of anaerobic glycolysis, in cancer cells3. However, how cancer metabolism affects chemotherapy response and DNA repair in general remains incompletely understood. Here we report that lactate-driven lactylation of NBS1 promotes homologous recombination (HR)-mediated DNA repair. Lactylation of NBS1 at lysine 388 (K388) is essential for MRE11-RAD50-NBS1 (MRN) complex formation and the accumulation of HR repair proteins at the sites of DNA double-strand breaks. Furthermore, we identify TIP60 as the NBS1 lysine lactyltransferase and the 'writer' of NBS1 K388 lactylation, and HDAC3 as the NBS1 de-lactylase. High levels of NBS1 K388 lactylation predict poor patient outcome of neoadjuvant chemotherapy, and lactate reduction using either genetic depletion of lactate dehydrogenase A (LDHA) or stiripentol, a lactate dehydrogenase A inhibitor used clinically for anti-epileptic treatment, inhibited NBS1 K388 lactylation, decreased DNA repair efficacy and overcame resistance to chemotherapy. In summary, our work identifies NBS1 lactylation as a critical mechanism for genome stability that contributes to chemotherapy resistance and identifies inhibition of lactate production as a promising therapeutic cancer strategy.

Coastal phytoplankton blooms expand and intensify in the 21st century.

Phytoplankton blooms in coastal oceans can be beneficial to coastal fisheries production and ecosystem function, but can also cause major environmental problems1,2-yet detailed characterizations of bloom incidence and distribution are not available worldwide. Here we map daily marine coastal algal blooms between 2003 and 2020 using global satellite observations at 1-km spatial resolution. We found that algal blooms occurred in 126 out of the 153 coastal countries examined. Globally, the spatial extent (+13.2%) and frequency (+59.2%) of blooms increased significantly (P < 0.05) over the study period, whereas blooms weakened in tropical and subtropical areas of the Northern Hemisphere. We documented the relationship between the bloom trends and ocean circulation, and identified the stimulatory effects of recent increases in sea surface temperature. Our compilation of daily mapped coastal phytoplankton blooms provides the basis for global assessments of bloom risks and benefits, and for the formulation or evaluation of management or policy actions.

Genome editing of a rice CDP-DAG synthase confers multipathogen resistance.

The discovery and application of genome editing introduced a new era of plant breeding by giving researchers efficient tools for the precise engineering of crop genomes1. Here we demonstrate the power of genome editing for engineering broad-spectrum disease resistance in rice (Oryza sativa). We first isolated a lesion mimic mutant (LMM) from a mutagenized rice population. We then demonstrated that a 29-base-pair deletion in a gene we named RESISTANCE TO BLAST1 (RBL1) caused broad-spectrum disease resistance and showed that this mutation caused an approximately 20-fold reduction in yield. RBL1 encodes a cytidine diphosphate diacylglycerol synthase that is required for phospholipid biosynthesis2. Mutation of RBL1 results in reduced levels of phosphatidylinositol and its derivative phosphatidylinositol 4,5-bisphosphate (PtdIns(4,5)P2). In rice, PtdIns(4,5)P2 is enriched in cellular structures that are specifically associated with effector secretion and fungal infection, suggesting that it has a role as a disease-susceptibility factor3. By using targeted genome editing, we obtained an allele of RBL1, named RBL1Δ12, which confers broad-spectrum disease resistance but does not decrease yield in a model rice variety, as assessed in small-scale field trials. Our study has demonstrated the benefits of editing an LMM gene, a strategy relevant to diverse LMM genes and crops.

Targeting Degradation of the Transcription Factor C/EBPß Reduces Lung Fibrosis by Restoring Activity of the Ubiquitin-Editing Enzyme A20 in Macrophages.

Although recent progress provides mechanistic insights into the pathogenesis of pulmonary fibrosis (PF), rare anti-PF therapeutics show definitive promise for treating this disease. Repeated lung epithelial injury results in injury-repairing response and inflammation, which drive the development of PF. Here, we report that chronic lung injury inactivated the ubiquitin-editing enzyme A20, causing progressive accumulation of the transcription factor C/EBPß in alveolar macrophages (AMs) from PF patients and mice, which upregulated a number of immunosuppressive and profibrotic factors promoting PF development. In response to chronic lung injury, elevated glycogen synthase kinase-3ß (GSK-3ß) interacted with and phosphorylated A20 to suppress C/EBPß degradation. Ectopic expression of A20 or pharmacological restoration of A20 activity by disturbing the A20-GSK-3ß interaction accelerated C/EBPß degradation and showed potent therapeutic efficacy against experimental PF. Our study indicates that a regulatory mechanism of the GSK-3ß-A20-C/EBPß axis in AMs may be a potential target for treating PF and fibroproliferative lung diseases.

Inhibition of ASGR1 decreases lipid levels by promoting cholesterol excretion.

High cholesterol is a major risk factor for cardiovascular disease1. Currently, no drug lowers cholesterol through directly promoting cholesterol excretion. Human genetic studies have identified that the loss-of-function Asialoglycoprotein receptor 1 (ASGR1) variants associate with low cholesterol and a reduced risk of cardiovascular disease2. ASGR1 is exclusively expressed in liver and mediates internalization and lysosomal degradation of blood asialoglycoproteins3. The mechanism by which ASGR1 affects cholesterol metabolism is unknown. Here, we find that Asgr1 deficiency decreases lipid levels in serum and liver by stabilizing LXRα. LXRα upregulates ABCA1 and ABCG5/G8, which promotes cholesterol transport to high-density lipoprotein and excretion to bile and faeces4, respectively. ASGR1 deficiency blocks endocytosis and lysosomal degradation of glycoproteins, reduces amino-acid levels in lysosomes, and thereby inhibits mTORC1 and activates AMPK. On one hand, AMPK increases LXRα by decreasing its ubiquitin ligases BRCA1/BARD1. On the other hand, AMPK suppresses SREBP1 that controls lipogenesis. Anti-ASGR1 neutralizing antibody lowers lipid levels by increasing cholesterol excretion, and shows synergistic beneficial effects with atorvastatin or ezetimibe, two widely used hypocholesterolaemic drugs. In summary, this study demonstrates that targeting ASGR1 upregulates LXRα, ABCA1 and ABCG5/G8, inhibits SREBP1 and lipogenesis, and therefore promotes cholesterol excretion and decreases lipid levels.

Robust Hi-C Maps of Enhancer-Promoter Interactions Reveal the Function of Non-coding Genome in Neural Development and Diseases.

Genome-wide mapping of chromatin interactions at high resolution remains experimentally and computationally challenging. Here we used a low-input "easy Hi-C" protocol to map the 3D genome architecture in human neurogenesis and brain tissues and also demonstrated that a rigorous Hi-C bias-correction pipeline (HiCorr) can significantly improve the sensitivity and robustness of Hi-C loop identification at sub-TAD level, especially the enhancer-promoter (E-P) interactions. We used HiCorr to compare the high-resolution maps of chromatin interactions from 10 tissue or cell types with a focus on neurogenesis and brain tissues. We found that dynamic chromatin loops are better hallmarks for cellular differentiation than compartment switching. HiCorr allowed direct observation of cell-type- and differentiation-specific E-P aggregates spanning large neighborhoods, suggesting a mechanism that stabilizes enhancer contacts during development. Interestingly, we concluded that Hi-C loop outperforms eQTL in explaining neurological GWAS results, revealing a unique value of high-resolution 3D genome maps in elucidating the disease etiology.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Enhanced mapping of small-molecule binding sites in cells.

Photoaffinity probes are routinely utilized to identify proteins that interact with small molecules. However, despite this common usage, resolving the specific sites of these interactions remains a challenge. Here we developed a chemoproteomic workflow to determine precise protein binding sites of photoaffinity probes in cells. Deconvolution of features unique to probe-modified peptides, such as their tendency to produce chimeric spectra, facilitated the development of predictive models to confidently determine labeled sites. This yielded an expansive map of small-molecule binding sites on endogenous proteins and enabled the integration with multiplexed quantitation, increasing the throughput and dimensionality of experiments. Finally, using structural information, we characterized diverse binding sites across the proteome, providing direct evidence of their tractability to small molecules. Together, our findings reveal new knowledge for the analysis of photoaffinity probes and provide a robust method for high-resolution mapping of reversible small-molecule interactions en masse in native systems.

Cell-type-specific 3D epigenomes in the developing human cortex.

Lineage-specific epigenomic changes during human corticogenesis have been difficult to study owing to challenges with sample availability and tissue heterogeneity. For example, previous studies using single-cell RNA sequencing identified at least 9 major cell types and up to 26 distinct subtypes in the dorsal cortex alone1,2. Here we characterize cell-type-specific cis-regulatory chromatin interactions, open chromatin peaks, and transcriptomes for radial glia, intermediate progenitor cells, excitatory neurons, and interneurons isolated from mid-gestational samples of the human cortex. We show that chromatin interactions underlie several aspects of gene regulation, with transposable elements and disease-associated variants enriched at distal interacting regions in a cell-type-specific manner. In addition, promoters with increased levels of chromatin interactivity-termed super-interactive promoters-are enriched for lineage-specific genes, suggesting that interactions at these loci contribute to the fine-tuning of transcription. Finally, we develop CRISPRview, a technique that integrates immunostaining, CRISPR interference, RNAscope, and image analysis to validate cell-type-specific cis-regulatory elements in heterogeneous populations of primary cells. Our findings provide insights into cell-type-specific gene expression patterns in the developing human cortex and advance our understanding of gene regulation and lineage specification during this crucial developmental window.

Genetic regulators of sputum mucin concentration and their associations with COPD phenotypes.

Hyper-secretion and/or hyper-concentration of mucus is a defining feature of multiple obstructive lung diseases, including chronic obstructive pulmonary disease (COPD). Mucus itself is composed of a mixture of water, ions, salt and proteins, of which the gel-forming mucins, MUC5AC and MUC5B, are the most abundant. Recent studies have linked the concentrations of these proteins in sputum to COPD phenotypes, including chronic bronchitis (CB) and acute exacerbations (AE). We sought to determine whether common genetic variants influence sputum mucin concentrations and whether these variants are also associated with COPD phenotypes, specifically CB and AE. We performed a GWAS to identify quantitative trait loci for sputum mucin protein concentration (pQTL) in the Sub-Populations and InteRmediate Outcome Measures in COPD Study (SPIROMICS, n = 708 for total mucin, n = 215 for MUC5AC, MUC5B). Subsequently, we tested for associations of mucin pQTL with CB and AE using regression modeling (n = 822-1300). Replication analysis was conducted using data from COPDGene (n = 5740) and by examining results from the UK Biobank. We identified one genome-wide significant pQTL for MUC5AC (rs75401036) and two for MUC5B (rs140324259, rs10001928). The strongest association for MUC5B, with rs140324259 on chromosome 11, explained 14% of variation in sputum MUC5B. Despite being associated with lower MUC5B, the C allele of rs140324259 conferred increased risk of CB (odds ratio (OR) = 1.42; 95% confidence interval (CI): 1.10-1.80) as well as AE ascertained over three years of follow up (OR = 1.41; 95% CI: 1.02-1.94). Associations between rs140324259 and CB or AE did not replicate in COPDGene. However, in the UK Biobank, rs140324259 was associated with phenotypes that define CB, namely chronic mucus production and cough, again with the C allele conferring increased risk. We conclude that sputum MUC5AC and MUC5B concentrations are associated with common genetic variants, and the top locus for MUC5B may influence COPD phenotypes, in particular CB.

Cargo-selective and adaptive delivery of nucleic acid therapeutics by bola-amphiphilic dendrimers.

Nucleic acid therapeutics are becoming an important drug modality, offering the unique opportunity to address "undruggable" targets, respond rapidly to evolving pathogens, and treat diseases at the gene level for precision medicine. However, nucleic acid therapeutics have poor bioavailability and are chemolabile and enzymolabile, imposing the need for delivery vectors. Dendrimers, by virtue of their well-defined structure and cooperative multivalence, represent precision delivery systems. We synthesized and studied bola-amphiphilic dendrimers for cargo-selective and on-demand delivery of DNA and small interfering RNA (siRNA), both important nucleic acid therapeutics. Remarkably, superior performances were achieved for siRNA delivery with the second-generation dendrimer, yet for DNA delivery with the third generation. We systematically studied these dendrimers with regard to cargo binding, cellular uptake, endosomal release, and in vivo delivery. Differences in size both of the dendrimers and their nucleic acid cargos impacted the cooperative multivalent interactions for cargo binding and release, leading to cargo-adaptive and selective delivery. Moreover, both dendrimers harnessed the advantages of lipid and polymer vectors, while offering nanotechnology-based tumor targeting and redox-responsive cargo release. Notably, they allowed tumor- and cancer cell-specific delivery of siRNA and DNA therapeutics for effective treatment in different cancer models, including aggressive and metastatic malignancies, outperforming the currently available vectors. This study provides avenues to engineer tailor-made vectors for nucleic acid delivery and precision medicine.

A spatiotemporal barrier formed by Follistatin is required for left-right patterning.

How left-right (LR) asymmetry emerges in a patterning field along the anterior-posterior axis remains an unresolved problem in developmental biology. Left-biased Nodal emanating from the LR organizer propagates from posterior to anterior (PA) and establishes the LR pattern of the whole embryo. However, little is known about the regulatory mechanism of the PA spread of Nodal and its asymmetric activation in the forebrain. Here, we identify bilaterally expressed Follistatin (Fst) as a regulator blocking the propagation of the zebrafish Nodal ortholog Southpaw (Spaw) in the right lateral plate mesoderm (LPM), and restricting Spaw transmission in the left LPM to facilitate the establishment of a robust LR asymmetric Nodal patterning. In addition, Fst inhibits the Activin-Nodal signaling pathway in the forebrain thus preventing Nodal activation prior to the arrival, at a later time, of Spaw emanating from the left LPM. This contributes to the orderly propagation of asymmetric Nodal activation along the PA axis. The LR regulation function of Fst is further confirmed in chick and frog embryos. Overall, our results suggest that a robust LR patterning emerges by counteracting a Fst barrier formed along the PA axis.

Canonical correlation analysis for multi-omics: Application to cross-cohort analysis.

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Proteome-wide association study using cis and trans variants and applied to blood cell and lipid-related traits in the Women's Health Initiative study.

In most Proteome-Wide Association Studies (PWAS), variants near the protein-coding gene (±1 Mb), also known as cis single nucleotide polymorphisms (SNPs), are used to predict protein levels, which are then tested for association with phenotypes. However, proteins can be regulated through variants outside of the cis region. An intermediate GWAS step to identify protein quantitative trait loci (pQTL) allows for the inclusion of trans SNPs outside the cis region in protein-level prediction models. Here, we assess the prediction of 540 proteins in 1002 individuals from the Women's Health Initiative (WHI), split equally into a GWAS set, an elastic net training set, and a testing set. We compared the testing r2 between measured and predicted protein levels using this proposed approach, to the testing r2 using only cis SNPs. The two methods usually resulted in similar testing r2, but some proteins showed a significant increase in testing r2 with our method. For example, for cartilage acidic protein 1, the testing r2 increased from 0.101 to 0.351. We also demonstrate reproducible findings for predicted protein association with lipid and blood cell traits in WHI participants without proteomics data and in UK Biobank utilizing our PWAS weights.

STAMENLESS1 activates SUPERWOMAN 1 and FLORAL ORGAN NUMBER 1 to control floral organ identities and meristem fate in rice.

Floral patterns are unique to rice and contribute significantly to its reproductive success. SL1 encodes a C2H2 transcription factor that plays a critical role in flower development in rice, but the molecular mechanism regulated by it remains poorly understood. Here, we describe interactions of the SL1 with floral homeotic genes, SPW1, and DL in specifying floral organ identities and floral meristem fate. First, the sl1 spw1 double mutant exhibited a stamen-to-pistil transition similar to that of sl1, spw1, suggesting that SL1 and SPW1 may located in the same pathway regulating stamen development. Expression analysis revealed that SL1 is located upstream of SPW1 to maintain its high level of expression and that SPW1, in turn, activates the B-class genes OsMADS2 and OsMADS4 to suppress DL expression indirectly. Secondly, sl1 dl displayed a severe loss of floral meristem determinacy and produced amorphous tissues in the third/fourth whorl. Expression analysis revealed that the meristem identity gene OSH1 was ectopically expressed in sl1 dl in the fourth whorl, suggesting that SL1 and DL synergistically terminate the floral meristem fate. Another meristem identity gene, FON1, was significantly decreased in expression in sl1 background mutants, suggesting that SL1 may directly activate its expression to regulate floral meristem fate. Finally, molecular evidence supported the direct genomic binding of SL1 to SPW1 and FON1 and the subsequent activation of their expression. In conclusion, we present a model to illustrate the roles of SL1, SPW1, and DL in floral organ specification and regulation of floral meristem fate in rice.