Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops.

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

Maize GOLDEN2-LIKE proteins enhance drought tolerance in rice by promoting stomatal closure.

Drought has become one of the most severe abiotic stresses experienced in agricultural production across the world. Plants respond to water deficit via stomatal movements in the leaves, which are mainly regulated by abscisic acid (ABA). A previous study from our lab showed that constitutive expression of maize (Zea mays L.) GOLDEN2-LIKE (GLK) transcription factors in rice (Oryza sativa L.) can improve stomatal conductance and plant photosynthetic capacity under field conditions. In the present study, we uncovered a function of ZmGLK regulation of stomatal movement in rice during drought stress. We found that elevated drought tolerance in rice plants overexpressing ZmGLK1 or GOLDEN2 (ZmG2) was conferred by rapid ABA-mediated stomatal closure. Comparative analysis of RNA-sequencing (RNA-seq) data from the rice leaves and DNA affinity purification sequencing (DAP-seq) results obtained in vitro revealed that ZmGLKs played roles in regulating ABA-related and stress-responsive pathways. Four upregulated genes closely functioning in abiotic stress tolerance with strong binding peaks in the DAP-seq data were identified as putative target genes of ZmGLK1 and ZmG2 in rice. These results demonstrated that maize GLKs play an important role in regulating stomatal movements to coordinate photosynthesis and stress tolerance. This trait is a valuable target for breeding drought-tolerant crop plants without compromising photosynthetic capacity.

Photoperiod-responsive changes in chromatin accessibility in phloem companion and epidermis cells of Arabidopsis leaves.

Photoperiod plays a key role in controlling the phase transition from vegetative to reproductive growth in flowering plants. Leaves are the major organs perceiving day-length signals, but how specific leaf cell types respond to photoperiod remains unknown. We integrated photoperiod-responsive chromatin accessibility and transcriptome data in leaf epidermis and vascular companion cells of Arabidopsis thaliana by combining isolation of nuclei tagged in specific cell/tissue types with assay for transposase-accessible chromatin using sequencing and RNA-sequencing. Despite a large overlap, vasculature and epidermis cells responded differently. Long-day predominantly induced accessible chromatin regions (ACRs); in the vasculature, more ACRs were induced and these were located at more distal gene regions, compared with the epidermis. Vascular ACRs induced by long days were highly enriched in binding sites for flowering-related transcription factors. Among the highly ranked genes (based on chromatin and expression signatures in the vasculature), we identified TREHALOSE-PHOSPHATASE/SYNTHASE 9 (TPS9) as a flowering activator, as shown by the late flowering phenotypes of T-DNA insertion mutants and transgenic lines with phloem-specific knockdown of TPS9. Our cell-type-specific analysis sheds light on how the long-day photoperiod stimulus impacts chromatin accessibility in a tissue-specific manner to regulate plant development.

Stable unmethylated DNA demarcates expressed genes and their cis-regulatory space in plant genomes.

The genomic sequences of crops continue to be produced at a frenetic pace. It remains challenging to develop complete annotations of functional genes and regulatory elements in these genomes. Chromatin accessibility assays enable discovery of functional elements; however, to uncover the full portfolio of cis-elements would require profiling of many combinations of cell types, tissues, developmental stages, and environments. Here, we explore the potential to use DNA methylation profiles to develop more complete annotations. Using leaf tissue in maize, we define ∼100,000 unmethylated regions (UMRs) that account for 5.8% of the genome; 33,375 UMRs are found greater than 2 kb from genes. UMRs are highly stable in multiple vegetative tissues, and they capture the vast majority of accessible chromatin regions from leaf tissue. However, many UMRs are not accessible in leaf, and these represent regions with potential to become accessible in specific cell types or developmental stages. These UMRs often occur near genes that are expressed in other tissues and are enriched for binding sites of transcription factors. The leaf-inaccessible UMRs exhibit unique chromatin modification patterns and are enriched for chromatin interactions with nearby genes. The total UMR space in four additional monocots ranges from 80 to 120 megabases, which is remarkably similar considering the range in genome size of 271 megabases to 4.8 gigabases. In summary, based on the profile from a single tissue, DNA methylation signatures provide powerful filters to distill large genomes down to the small fraction of putative functional genes and regulatory elements.

HY5-HDA9 orchestrates the transcription of HsfA2 to modulate salt stress response in Arabidopsis.

Integration of light signaling and diverse abiotic stress responses contribute to plant survival in a changing environment. Some reports have indicated that light signals contribute a plant's ability to deal with heat, cold, and stress. However, the molecular link between light signaling and the salt-response pathways remains unclear. We demonstrate here that increasing light intensity elevates the salt stress tolerance of plants. Depletion of HY5, a key component of light signaling, causes Arabidopsis thaliana to become salinity sensitive. Interestingly, the small heat shock protein (sHsp) family genes are upregulated in hy5-215 mutant plants, and HsfA2 is commonly involved in the regulation of these sHsps. We found that HY5 directly binds to the G-box motifs in the HsfA2 promoter, with the cooperation of HISTONE DEACETYLASE 9 (HDA9), to repress its expression. Furthermore, the accumulation of HDA9 and the interaction between HY5 and HDA9 are significantly enhanced by salt stress. On the contrary, high temperature triggers HY5 and HDA9 degradation, which leads to dissociation of HY5-HDA9 from the HsfA2 promoter, thereby reducing salt tolerance. Under salt and heat stress conditions, fine tuning of protein accumulation and an interaction between HY5 and HDA9 regulate HsfA2 expression. This implies that HY5, HDA9, and HsfA2 play important roles in the integration of light signaling with salt stress and heat shock response.

Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana.

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

Monitoring the interplay between transposable element families and DNA methylation in maize.

DNA methylation and epigenetic silencing play important roles in the regulation of transposable elements (TEs) in many eukaryotic genomes. A majority of the maize genome is derived from TEs that can be classified into different orders and families based on their mechanism of transposition and sequence similarity, respectively. TEs themselves are highly methylated and it can be tempting to view them as a single uniform group. However, the analysis of DNA methylation profiles in flanking regions provides evidence for distinct groups of chromatin properties at different TE families. These differences among TE families are reproducible in different tissues and different inbred lines. TE families with varying levels of DNA methylation in flanking regions also show distinct patterns of chromatin accessibility and modifications within the TEs. The differences in the patterns of DNA methylation flanking TE families arise from a combination of non-random insertion preferences of TE families, changes in DNA methylation triggered by the insertion of the TE and subsequent selection pressure. A set of nearly 70,000 TE polymorphisms among four assembled maize genomes were used to monitor the level of DNA methylation at haplotypes with and without the TE insertions. In many cases, TE families with high levels of DNA methylation in flanking sequence are enriched for insertions into highly methylated regions. The majority of the >2,500 TE insertions into unmethylated regions result in changes in DNA methylation in haplotypes with the TE, suggesting the widespread potential for TE insertions to condition altered methylation in conserved regions of the genome. This study highlights the interplay between TEs and the methylome of a major crop species.

Chromatin accessibility profiling in Neurospora crassa reveals molecular features associated with accessible and inaccessible chromatin.

BACKGROUND: Regulation of chromatin accessibility and transcription are tightly coordinated processes. Studies in yeast and higher eukaryotes have described accessible chromatin regions, but little work has been done in filamentous fungi. RESULTS: Here we present a genome-scale characterization of accessible chromatin regions in Neurospora crassa, which revealed characteristic molecular features of accessible and inaccessible chromatin. We present experimental evidence of inaccessibility within heterochromatin regions in Neurospora, and we examine features of both accessible and inaccessible chromatin, including the presence of histone modifications, types of transcription, transcription factor binding, and relative nucleosome turnover rates. Chromatin accessibility is not strictly correlated with expression level. Accessible chromatin regions in the model filamentous fungus Neurospora are characterized the presence of H3K27 acetylation and commonly associated with pervasive non-coding transcription. Conversely, methylation of H3 lysine-36 catalyzed by ASH1 is correlated with inaccessible chromatin within promoter regions. CONCLUSIONS: In N. crassa, H3K27 acetylation is the most predictive histone modification for open chromatin. Conversely, our data show that H3K36 methylation is a key marker of inaccessible chromatin in gene-rich regions of the genome. Our data are consistent with an expanded role for H3K36 methylation in intergenic regions of filamentous fungi compared to the model yeasts, S. cerevisiae and S. pombe, which lack homologs of the ASH1 methyltransferase.

Tissue-Specific Ubiquitination by IPA1 INTERACTING PROTEIN1 Modulates IPA1 Protein Levels to Regulate Plant Architecture in Rice.

Plant architecture, a collection of genetically controlled agronomic traits, is one of the decisive factors that determine grain production. IDEAL PLANT ARCHITECTURE1 (IPA1) encodes a key transcription factor with pleiotropic effects on regulating plant architecture in rice (Oryza sativa), and IPA1 expression is controlled at the posttranscriptional level by microRNA156 and microRNA529. Here, we report the identification and characterization of IPA1 INTERACTING PROTEIN1 (IPI1), a RING-finger E3 ligase that can interact with IPA1 in the nucleus. IPI1 promotes the degradation of IPA1 in panicles, while it stabilizes IPA1 in shoot apexes. Consistent with these findings, the ipi1 loss-of-function mutants showed markedly altered plant architecture, including more tillers, enlarged panicles, and increased yield per plant. Moreover, IPI1 could ubiquitinate the IPA1-mediated complex with different polyubiquitin chains, adding K48-linked polyubiquitin chains in panicles and K63-linked polyubiquitin chains in the shoot apex. These results demonstrate that IPI1 affects plant architecture through precisely tuning IPA1 protein levels in different tissues in rice and provide new insight into the tissue-specific regulation of plant architecture and important genetic resources for molecular breeding.

Combining ATAC-seq with nuclei sorting for discovery of cis-regulatory regions in plant genomes.

Chromatin structure plays a pivotal role in facilitating proper control of gene expression. Transcription factor (TF) binding of cis-elements is often associated with accessible chromatin regions. Therefore, the ability to identify these accessible regions throughout plant genomes will advance understanding of the relationship between TF binding, chromatin status and the regulation of gene expression. Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is a recently developed technique used to map open chromatin zones in animal genomes. However, in plants, the existence of cell walls, subcellular organelles and the lack of stable cell lines have prevented routine application of this technique. Here, we describe an assay combining ATAC-seq with fluorescence-activated nuclei sorting (FANS) to identify and map open chromatin and TF-binding sites in plant genomes. FANS-ATAC-seq compares favorably with published DNaseI sequencing (DNase-seq) results and it requires less than 50 000 nuclei for accurate identification of accessible genomic regions. SUMMARY: Application of ATAC-seq to sorted nuclei identifies accessible regions genome-wide.

On the origin and evolutionary consequences of gene body DNA methylation.

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Strigolactone Signaling in Arabidopsis Regulates Shoot Development by Targeting D53-Like SMXL Repressor Proteins for Ubiquitination and Degradation.

Strigolactones (SLs) are carotenoid-derived phytohormones that control many aspects of plant development, including shoot branching, leaf shape, stem secondary thickening, and lateral root growth. In rice (Oryza sativa), SL signaling requires the degradation of DWARF53 (D53), mediated by a complex including D14 and D3, but in Arabidopsis thaliana, the components and mechanism of SL signaling involving the D3 ortholog MORE AXILLARY GROWTH2 (MAX2) are unknown. Here, we show that SL-dependent regulation of shoot branching in Arabidopsis requires three D53-like proteins, SUPPRESSOR OF MORE AXILLARY GROWTH2-LIKE6 (SMXL6), SMXL7, and SMXL8. The smxl6 smxl7 smxl8 triple mutant suppresses the highly branched phenotypes of max2 and the SL-deficient mutant max3. Overexpression of a mutant form of SMXL6 that is resistant to SL-induced ubiquitination and degradation enhances shoot branching. Exogenous application of the SL analog rac-GR24 causes ubiquitination and degradation of SMXL6, 7, and 8; this requires D14 and MAX2. D53-like SMXLs form complexes with MAX2 and TOPLESS-RELATED PROTEIN2 (TPR2) and interact with D14 in a GR24-responsive manner. Furthermore, D53-like SMXLs exhibit TPR2-dependent transcriptional repression activity and repress the expression of BRANCHED1. Our findings reveal that in Arabidopsis, D53-like SMXLs act with TPR2 to repress transcription and so allow lateral bud outgrowth but that SL-induced degradation of D53-like proteins activates transcription to inhibit outgrowth.

Genome-wide binding analysis of the transcription activator ideal plant architecture1 reveals a complex network regulating rice plant architecture.

Ideal plant architecture1 (IPA1) is critical in regulating rice (Oryza sativa) plant architecture and substantially enhances grain yield. To elucidate its molecular basis, we first confirmed IPA1 as a functional transcription activator and then identified 1067 and 2185 genes associated with IPA1 binding sites in shoot apices and young panicles, respectively, through chromatin immunoprecipitation sequencing assays. The Squamosa promoter binding protein-box direct binding core motif GTAC was highly enriched in IPA1 binding peaks; interestingly, a previously uncharacterized indirect binding motif TGGGCC/T was found to be significantly enriched through the interaction of IPA1 with proliferating cell nuclear antigen promoter binding factor1 or promoter binding factor2. Genome-wide expression profiling by RNA sequencing revealed IPA1 roles in diverse pathways. Moreover, our results demonstrated that IPA1 could directly bind to the promoter of rice teosinte branched1, a negative regulator of tiller bud outgrowth, to suppress rice tillering, and directly and positively regulate dense and erect panicle1, an important gene regulating panicle architecture, to influence plant height and panicle length. The elucidation of target genes of IPA1 genome-wide will contribute to understanding the molecular mechanisms underlying plant architecture and to facilitating the breeding of elite varieties with ideal plant architecture.

The unusual predominance of maintenance DNA methylation in Spirodela polyrhiza.

Duckweeds are among the fastest reproducing plants, able to clonally divide at exponential rates. However, the genetic and epigenetic impact of clonality on plant genomes is poorly understood. 5-methylcytosine (5mC) is a modified base often described as necessary for the proper regulation of certain genes and transposons and for the maintenance of genome integrity in plants. However, the extent of this dogma is limited by the current phylogenetic sampling of land plant species diversity. Here we analyzed DNA methylomes, small RNAs, mRNA-seq, and H3K9me2 histone modification for Spirodela polyrhiza. S. polyrhiza has lost highly conserved genes involved in de novo methylation of DNA at sites often associated with repetitive DNA, and within genes, however, symmetrical DNA methylation and heterochromatin are maintained during cell division at certain transposons and repeats. Consequently, small RNAs that normally guide methylation to silence repetitive DNA like retrotransposons are diminished. Despite the loss of a highly conserved methylation pathway, and the reduction of small RNAs that normally target repetitive DNA, transposons have not proliferated in the genome, perhaps due in part to the rapid, clonal growth lifestyle of duckweeds.

Enhancing rice panicle branching and grain yield through tissue-specific brassinosteroid inhibition.

Crop yield potential is constrained by the inherent trade-offs among traits such as between grain size and number. Brassinosteroids (BRs) promote grain size, yet their role in regulating grain number is unclear. By deciphering the clustered-spikelet rice germplasm, we show that activation of the BR catabolic gene BRASSINOSTEROID-DEFICIENT DWARF3 (BRD3) markedly increases grain number. We establish a molecular pathway in which the BR signaling inhibitor GSK3/SHAGGY-LIKE KINASE2 phosphorylates and stabilizes OsMADS1 transcriptional factor, which targets TERMINAL FLOWER1-like gene RICE CENTRORADIALIS2. The tissue-specific activation of BRD3 in the secondary branch meristems enhances panicle branching, minimizing negative effects on grain size, and improves grain yield. Our study showcases the power of tissue-specific hormonal manipulation in dismantling the trade-offs among various traits and thus unleashing crop yield potential in rice.

Chromatin accessibility landscapes revealed the subgenome-divergent regulation networks during wheat grain development.

Development of wheat (Triticum aestivum L.) grain mainly depends on the processes of starch synthesis and storage protein accumulation, which are critical for grain yield and quality. However, the regulatory network underlying the transcriptional and physiological changes of grain development is still not clear. Here, we combined ATAC-seq and RNA-seq to discover the chromatin accessibility and gene expression dynamics during these processes. We found that the chromatin accessibility changes are tightly associated with differential transcriptomic expressions, and the proportion of distal ACRs was increased gradually during grain development. Specific transcription factor (TF) binding sites were enriched at different stages and were diversified among the 3 subgenomes. We further predicted the potential interactions between key TFs and genes related with starch and storage protein biosynthesis and found different copies of some key TFs played diversified roles. Overall, our findings have provided numerous resources and illustrated the regulatory network during wheat grain development, which would shed light on the improvement of wheat yields and qualities. Supplementary Information: The online version contains supplementary material available at 10.1007/s42994-023-00095-8.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat.

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat.

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

Centromere repositioning and shifts in wheat evolution.

The centromere is the region of a chromosome that directs its separation and plays an important role in cell division and reproduction of organisms. Elucidating the dynamics of centromeres is an alternative strategy for exploring the evolution of wheat. Here, we comprehensively analyzed centromeres from the de novo-assembled common wheat cultivar Aikang58 (AK58), Chinese Spring (CS), and all sequenced diploid and tetraploid ancestors by chromatin immunoprecipitation sequencing, whole-genome bisulfite sequencing, RNA sequencing, assay for transposase-accessible chromatin using sequencing, and comparative genomics. We found that centromere-associated sequences were concentrated during tetraploidization and hexaploidization. Centromeric repeats of wheat (CRWs) have undergone expansion during wheat evolution, with strong interweaving between the A and B subgenomes post tetraploidization. We found that CENH3 prefers to bind with younger CRWs, as directly supported by immunocolocalization on two chromosomes (1A and 2A) of wild emmer wheat with dicentromeric regions, only one of which bound with CENH3. In a comparison of AK58 with CS, obvious centromere repositioning was detected on chromosomes 1B, 3D, and 4D. The active centromeres showed a unique combination of lower CG but higher CHH and CHG methylation levels. We also found that centromeric chromatin was more open than pericentromeric chromatin, with higher levels of gene expression but lower gene density. Frequent introgression between tetraploid and hexaploid wheat also had a strong influence on centromere position on the same chromosome. This study also showed that active wheat centromeres were genetically and epigenetically determined.

Low phosphorus promotes NSP1-NSP2 heterodimerization to enhance strigolactone biosynthesis and regulate shoot and root architecture in rice.

Phosphorus is an essential macronutrient for plant development and metabolism, and plants have evolved ingenious mechanisms to overcome phosphate (Pi) starvation. However, the molecular mechanisms underlying the regulation of shoot and root architecture by low phosphorus conditions and the coordinated utilization of Pi and nitrogen remain largely unclear. Here, we show that Nodulation Signaling Pathway 1 (NSP1) and NSP2 regulate rice tiller number by promoting the biosynthesis of strigolactones (SLs), a class of phytohormones with fundamental effects on plant architecture and environmental responses. We found that NSP1 and NSP2 are induced by Oryza sativa PHOSPHATE STARVATION RESPONSE2 (OsPHR2) in response to low-Pi stress and form a complex to directly bind the promoters of SL biosynthesis genes, thus markedly increasing SL biosynthesis in rice. Interestingly, the NSP1/2-SL signaling module represses the expression of CROWN ROOTLESS 1 (CRL1), a newly identified early SL-responsive gene in roots, to restrain lateral root density under Pi deficiency. We also demonstrated that GR244DO treatment under normal conditions inhibits the expression of OsNRTs and OsAMTs to suppress nitrogen absorption but enhances the expression of OsPTs to promote Pi absorption, thus facilitating the balance between nitrogen and phosphorus uptake in rice. Importantly, we found that NSP1p:NSP1 and NSP2p:NSP2 transgenic plants show improved agronomic traits and grain yield under low- and medium-phosphorus conditions. Taken together, these results revealed a novel regulatory mechanism of SL biosynthesis and signaling in response to Pi starvation, providing genetic resources for improving plant architecture and nutrient-use efficiency in low-Pi environments.