A plant immune protein enables broad antitumor response by rescuing microRNA deficiency.

Cancer cells are featured with uncontrollable activation of cell cycle, and microRNA deficiency drives tumorigenesis. The RNA-dependent RNA polymerase (RDR) is essential for small-RNA-mediated immune response in plants but is absent in vertebrates. Here, we show that ectopic expression of plant RDR1 can generally inhibit cancer cell proliferation. In many human primary tumors, abnormal microRNA isoforms with 1-nt-shorter 3' ends are widely accumulated. RDR1 with nucleotidyltransferase activity can recognize and modify the problematic AGO2-free microRNA duplexes with mononucleotides to restore their 2 nt overhang structure, which eventually rescues AGO2-loading efficiency and elevates global miRNA expression to inhibit cancer cell-cycle specifically. The broad antitumor effects of RDR1, which can be delivered by an adeno-associated virus, are visualized in multiple xenograft tumor models in vivo. Altogether, we reveal the widespread accumulation of aberrant microRNA isoforms in tumors and develop a plant RDR1-mediated antitumor stratagem by editing and repairing defective microRNAs.

Paradoxical Role for Wild-Type p53 in Driving Therapy Resistance in Melanoma.

Metastatic melanoma is an aggressive disease, despite recent improvements in therapy. Eradicating all melanoma cells even in drug-sensitive tumors is unsuccessful in patients because a subset of cells can transition to a slow-cycling state, rendering them resistant to most targeted therapy. It is still unclear what pathways define these subpopulations and promote this resistant phenotype. In the current study, we show that Wnt5A, a non-canonical Wnt ligand that drives a metastatic, therapy-resistant phenotype, stabilizes the half-life of p53 and uses p53 to initiate a slow-cycling state following stress (DNA damage, targeted therapy, and aging). Inhibiting p53 blocks the slow-cycling phenotype and sensitizes melanoma cells to BRAF/MEK inhibition. In vivo, this can be accomplished with a single dose of p53 inhibitor at the commencement of BRAF/MEK inhibitor therapy. These data suggest that taking the paradoxical approach of inhibiting rather than activating wild-type p53 may sensitize previously resistant metastatic melanoma cells to therapy.

UB2/UB3/TSH4-anchored transcriptional networks regulate early maize inflorescence development in response to simulated shade.

Increasing planting density has been adopted as an effective means to increase maize (Zea mays) yield. Competition for light from neighbors can trigger plant shade avoidance syndrome, which includes accelerated flowering. However, the regulatory networks of maize inflorescence development in response to high-density planting remain poorly understood. In this study, we showed that shade-mimicking treatments cause precocious development of the tassels and ears. Comparative transcriptome profiling analyses revealed the enrichment of phytohormone-related genes and transcriptional regulators among the genes co-regulated by developmental progression and simulated shade. Network analysis showed that three homologous Squamosa promoter binding protein (SBP)-like (SPL) transcription factors, Unbranched2 (UB2), Unbranched3 (UB3), and Tasselsheath4 (TSH4), individually exhibited connectivity to over 2,400 genes across the V3-to-V9 stages of tassel development. In addition, we showed that the ub2 ub3 double mutant and tsh4 single mutant were almost insensitive to simulated shade treatments. Moreover, we demonstrated that UB2/UB3/TSH4 could directly regulate the expression of Barren inflorescence2 (BIF2) and Zea mays teosinte branched1/cycloidea/proliferating cell factor30 (ZmTCP30). Furthermore, we functionally verified a role of ZmTCP30 in regulating tassel branching and ear development. Our results reveal a UB2/UB3/TSH4-anchored transcriptional regulatory network of maize inflorescence development and provide valuable targets for breeding shade-tolerant maize cultivars.

Double-strand break repair-independent role for BRCA2 in blocking stalled replication fork degradation by MRE11.

Breast cancer suppressor BRCA2 is critical for maintenance of genomic integrity and resistance to agents that damage DNA or collapse replication forks, presumably through homology-directed repair of double-strand breaks (HDR). Using single-molecule DNA fiber analysis, we show here that nascent replication tracts created before fork stalling with hydroxyurea are degraded in the absence of BRCA2 but are stable in wild-type cells. BRCA2 mutational analysis reveals that a conserved C-terminal site involved in stabilizing RAD51 filaments, but not in loading RAD51 onto DNA, is essential for this fork protection but dispensable for HDR. RAD51 filament disruption in wild-type cells phenocopies BRCA2 deficiency. BRCA2 prevents chromosomal aberrations on replication stalling, which are alleviated by inhibition of MRE11, the nuclease responsible for this form of fork instability. Thus, BRCA2 prevents rather than repairs nucleolytic lesions at stalled replication forks to maintain genomic integrity and hence likely suppresses tumorigenesis through this replication-specific function.

Genomic analysis reveals a cryptic pangolin species.

Eight extant species of pangolins are currently recognized. Recent studies found that two mitochondrial haplotypes identified in confiscations in Hong Kong could not be assigned to any known pangolin species, implying the existence of a species. Here, we report that two additional mitochondrial haplotypes identified in independent confiscations from Yunnan align with the putative species haplotypes supporting the existence of this mysterious species/population. To verify the new species scenario we performed a comprehensive analysis of scale characteristics and 138 whole genomes representing all recognized pangolin species and the cryptic new species, 98 of which were generated here. Our morphometric results clearly attributed this cryptic species to Asian pangolins (Manis sp.) and the genomic data provide robust and compelling evidence that it is a pangolin species distinct from those recognized previously, which separated from the Philippine pangolin and Malayan pangolin over 5 Mya. Our study provides a solid genomic basis for its formal recognition as the ninth pangolin species or the fifth Asian one, supporting a new taxonomic classification of pangolins. The effects of glacial climate changes and recent anthropogenic activities driven by illegal trade are inferred to have caused its population decline with the genomic signatures showing low genetic diversity, a high level of inbreeding, and high genetic load. Our finding greatly expands current knowledge of pangolin diversity and evolution and has vital implications for conservation efforts to prevent the extinction of this enigmatic and endangered species from the wild.

Tiam1 methylation by NSD2 promotes Rac1 signaling activation and colon cancer metastasis.

Metastasis is a major cause of cancer therapy failure and mortality. However, targeting metastatic seeding and colonization remains a significant challenge. In this study, we identified NSD2, a histone methyltransferase responsible for dimethylating histone 3 at lysine 36, as being overexpressed in metastatic tumors. Our findings suggest that NSD2 overexpression enhances tumor metastasis both in vitro and in vivo. Further analysis revealed that NSD2 promotes tumor metastasis by activating Rac1 signaling. Mechanistically, NSD2 combines with and activates Tiam1 (T lymphoma invasion and metastasis 1) and promotes Rac1 signaling by methylating Tiam1 at K724. In vivo and in vitro studies revealed that Tiam1 K724 methylation could be a predictive factor for cancer prognosis and a potential target for metastasis inhibition. Furthermore, we have developed inhibitory peptide which was proved to inhibit tumor metastasis through blocking the interaction between NSD2 and Tiam1. Our results demonstrate that NSD2-methylated Tiam1 promotes Rac1 signaling and cancer metastasis. These results provide insights into the inhibition of tumor metastasis.

A prominent gene activation role for C-terminal binding protein in mediating PcG/trxG proteins through Hox gene regulation.

The evolutionarily conserved C-terminal binding protein (CtBP) has been well characterized as a transcriptional co-repressor. Herein, we report a previously unreported function for CtBP, showing that lowering CtBP dosage genetically suppresses Polycomb group (PcG) loss-of-function phenotypes while enhancing that of trithorax group (trxG) in Drosophila, suggesting that the role of CtBP in gene activation is more pronounced in fly development than previously thought. In fly cells, we show that CtBP is required for the derepression of the most direct PcG target genes, which are highly enriched by homeobox transcription factors, including Hox genes. Using ChIP and co-IP assays, we demonstrate that CtBP is directly required for the molecular switch between H3K27me3 and H3K27ac in the derepressed Hox loci. In addition, CtBP physically interacts with many proteins, such as UTX, CBP, Fs(1)h and RNA Pol II, that have activation roles, potentially assisting in their recruitment to promoters and Polycomb response elements that control Hox gene expression. Therefore, we reveal a prominent activation function for CtBP that confers a major role for the epigenetic program of fly segmentation and development.

BLISTER promotes seed maturation and fatty acid biosynthesis by interacting with WRINKLED1 to regulate chromatin dynamics in Arabidopsis.

WRINKLED1 (WRI1) is an important transcription factor that regulates seed oil biosynthesis. However, how WRI1 regulates gene expression during this process remains poorly understood. Here, we found that BLISTER (BLI) is expressed in maturing Arabidopsis thaliana seeds and acts as an interacting partner of WRI1. bli mutant seeds showed delayed maturation, a wrinkled seed phenotype, and reduced oil content, similar to the phenotypes of wri1. In contrast, BLI overexpression resulted in enlarged seeds and increased oil content. Gene expression and genetic analyses revealed that BLI plays a role in promoting the expression of WRI1 targets involved in fatty acid biosynthesis and regulates seed maturation together with WRI1. BLI is recruited by WRI1 to the AW boxes in the promoters of fatty acid biosynthesis genes. BLI shows a mutually exclusive interaction with the Polycomb-group protein CURLY LEAF (CLF) or the chromatin remodeling factor SWITCH/SUCROSE NONFERMENTING 3B (SWI3B), which facilitates gene expression by modifying nucleosomal occupancy and histone modifications. Together, these data suggest that BLI promotes the expression of fatty acid biosynthesis genes by interacting with WRI1 to regulate chromatin dynamics, leading to increased fatty acid production. These findings provide insights into the roles of the WRI1-BLI-CLF-SWI3B module in mediating seed maturation and gene expression.

Uncovering the enigmatic evolution of bears in greater depth: The hybrid origin of the Asiatic black bear.

Bears are fascinating mammals because of their complex pattern of speciation and rapid evolution of distinct phenotypes. Interspecific hybridization has been common and has shaped the complex evolutionary history of bears. In this study, based on the largest population-level genomic dataset to date involving all Ursinae species and recently developed methods for detecting hybrid speciation, we provide explicit evidence for the hybrid origin of Asiatic black bears, which arose through historical hybridization between the ancestor of polar bear/brown bear/American black bears and the ancestor of sun bear/sloth bears. This was inferred to have occurred soon after the divergence of the two parental lineages in Eurasia due to climate-driven population expansion and dispersal. In addition, we found that the intermediate body size of this hybrid species arose from its combination of relevant genes derived from two parental lineages of contrasting sizes. This and alternate fixation of numerous other loci that had diverged between parental lineages may have initiated the reproductive isolation of the Asiatic black bear from its two parents. Our study sheds further light on the evolutionary history of bears and documents the importance of hybridization in new species formation and phenotypic evolution in mammals.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales.

BACKGROUND: The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants. RESULTS: We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection. CONCLUSIONS: The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Fluorescent chemosensors facilitate the visualization of plant health and their living environment in sustainable agriculture.

Globally, 91% of plant production encounters diverse environmental stresses that adversely affect their growth, leading to severe yield losses of 50-60%. In this case, monitoring the connection between the environment and plant health can balance population demands with environmental protection and resource distribution. Fluorescent chemosensors have shown great progress in monitoring the health and environment of plants due to their high sensitivity and biocompatibility. However, to date, no comprehensive analysis and systematic summary of fluorescent chemosensors used in monitoring the correlation between plant health and their environment have been reported. Thus, herein, we summarize the current fluorescent chemosensors ranging from their design strategies to applications in monitoring plant-environment interaction processes. First, we highlight the types of fluorescent chemosensors with design strategies to resolve the bottlenecks encountered in monitoring the health and living environment of plants. In addition, the applications of fluorescent small-molecule, nano and supramolecular chemosensors in the visualization of the health and living environment of plants are discussed. Finally, the major challenges and perspectives in this field are presented. This work will provide guidance for the design of efficient fluorescent chemosensors to monitor plant health, and then promote sustainable agricultural development.

A LcDOF5.6-LcRbohD regulatory module controls the reactive oxygen species-mediated fruitlet abscission in litchi.

Reactive oxygen species (ROS) have been emerging as a key regulator in plant organ abscission. However, the mechanism underlying the regulation of ROS homeostasis in the abscission zone (AZ) is not completely established. Here, we report that a DOF (DNA binding with one finger) transcription factor LcDOF5.6 can suppress the litchi fruitlet abscission through repressing the ROS accumulation in fruitlet AZ (FAZ). The expression of LcRbohD, a homolog of the Arabidopsis RBOHs that are critical for ROS production, was significantly increased during the litchi fruitlet abscission, in parallel with an increased accumulation of ROS in FAZ. In contrast, silencing of LcRbohD reduced the ROS accumulation in FAZ and decreased the fruitlet abscission in litchi. Using in vitro and in vivo assays, we revealed that LcDOF5.6 was shown to inhibit the expression of LcRbohD via direct binding to its promoter. Consistently, silencing of LcDOF5.6 increased the expression of LcRbohD, concurrently with higher ROS accumulation in FAZ and increased fruitlet abscission. Furthermore, the expression of key genes (LcIDL1, LcHSL2, LcACO2, LcACS1, and LcEIL3) in INFLORESCENCE DEFICIENT IN ABSCISSION signaling and ethylene pathways were altered in LcRbohD-silenced and LcDOF5.6-silenced FAZ cells. Taken together, our results demonstrate an important role of the LcDOF5.6-LcRbohD module during litchi fruitlet abscission. Our findings provide new insights into the molecular regulatory network of organ abscission.

Therapeutic vaccination with lentiviral vector in HBV-persistent mice and two inactive HBsAg carriers.

BACKGROUND & AIMS: Immunotherapy for chronic hepatitis B virus (HBV) infection has not yet demonstrated sufficient efficacy. We developed a non-integrative lentiviral-vectored therapeutic vaccine for chronic hepatitis B and tested its antiviral effects in HBV-persistent mice and two inactive HBsAg carriers. METHODS: Lentiviral vectors (LVs) encoding the core, preS1, or large HBsAg (LHBs) proteins of HBV were evaluated for immunogenicity in HBV-naïve mice and therapeutic efficacy in a murine model of chronic HBV infection. In addition, two inactive HBsAg carriers each received two doses of 5×107 transduction units (TU) or 1×108 TU of lentiviral-vectored LHBs (LV-LHBs), respectively. The endpoints were safety, LHBs-specific T-cell responses, and serum HBsAg levels during a 24-week follow-up. RESULTS: In the mouse models, LV-LHBs was the most promising in eliciting robust antigen-specific T cells and in reducing the levels of serum HBsAg and viral load. By the end of the 34-week observation period, six out of ten (60%) HBV-persistent mice vaccinated with LV-LHBs achieved serum HBsAg loss and significant depletion of HBV-positive hepatocytes in the liver. In the two inactive HBsAg carriers, vaccination with LV-LHBs induced a considerable increase in the number of peripheral LHBs-specific T cells in one patient, and a weak but detectable response in the other, accompanied by a sustained reduction of HBsAg (-0.31 log10 IU/ml and -0.46 log10 IU/ml, respectively) from baseline to nadir. CONCLUSIONS: A lentiviral-vectored therapeutic vaccine for chronic HBV infection demonstrated the potential to improve HBV-specific T-cell responses and deplete HBV-positive hepatocytes, leading to a sustained loss or reduction of serum HBsAg. IMPACT AND IMPLICATIONS: Chronic HBV infection is characterized by an extremely low number and profound hypo-responsiveness of HBV-specific T cells. Therapeutic vaccines are designed to improve HBV-specific T-cell responses. We show that immunization with a lentiviral-vectored therapeutic HBV vaccine was able to expand HBV-specific T cells in vivo, leading to reductions of HBV-positive hepatocytes and serum HBsAg.

SIRPG expression positively associates with an inflamed tumor microenvironment and response to PD-1 blockade.

BACKGROUND: This study aimed to investigate the relationship between signal regulatory protein gamma (SIRPG) and tumor immune microenvironment phenotypes or T cell mediated-adaptive antitumor immunity, and its predictive value for response to PD-1 blockade in cancers. METHODS: Pan-cancer analysis of SIRPG expression and immune deconvolution was performed using transcriptomic data across 33 tumor types. Transcriptomic and clinical data from 157 patients with non-small-cell lung cancer (NSCLC) and melanoma received PD-1 blockade were analyzed. Expression characteristics of SIRPG were investigated using single-cell RNA sequencing (scRNA-seq) data of 103,599 cells. The effect of SIRPG expression was evaluated via SIRPG knockdown or overexpression in Jurkat T cells. RESULTS: The results showed that most cancers with high SIRPG expression had significantly higher abundance of T cells, B cells, NK cells, M1 macrophages and cytotoxic lymphocytes and increased expression level of immunomodulatory factors regulating immune cell recruitment, antigen presentation, T cell activation and cytotoxicity, but markedly lower abundance of neutrophils, M2 macrophages, and myeloid-derived suppressor cells. High SIRPG expression was associated with favorable response to PD-1 blockade in both NSCLC and melanoma. scRNA-seq data suggested SIRPG was mainly expressed in CD8+ exhausted T and CD4+ regulatory T cells, and positively associated with immune checkpoint expression including PDCD1 and CTLA4. In vitro test showed SIRPG expression in T cells could facilitate expression of PDCD1 and CTLA4. CONCLUSION: High SIRPG expression is associated with an inflamed immune phenotype in cancers and favorable response to PD-1 blockade, suggesting it would be a promising predictive biomarker for PD-1 blockade and novel immunotherapeutic target.

Modifying electronic and optical properties of violet phosphorus through variable fluorine coverage.

We present a comprehensive investigation of the electronic properties of fluorinated monolayer violet phosphorus using first-principles calculations. Our results reveal a strong dependence of the electronic properties on the different fluorine coverages of fluorination. As the fluorine coverage increases, monolayer violet phosphorus undergoes a significant transition from a wide direct bandgap semiconductor to a narrow indirect bandgap semiconductor. Moreover, both semi-fluorinated and fully fluorinated monolayer violet phosphorus exhibit advantageous semiconducting characteristics, with a tunable bandgap of 0.50 ~ 1.04 eV under biaxial strain ranging from -6% to 6%. Notably, the fully fluorinated monolayer violet phosphorus demonstrates a higher coefficient of light absorption within the visible range. Therefore, our findings highlight the tunability of monolayer violet phosphorus properties through the absorption of various fluorine coverages, providing valuable insights for the design and development of novel semiconductor devices based on this material.

Biomineralization-Inspired Copper Sulfide Decorated Aramid Textiles via In Situ Anchoring toward Versatile Wearable Thermal Management.

Designing smart textiles for personal thermal management (PTM) is an effective strategy for thermoregulation and energy saving. However, the manufacture of versatile high-performance thermal management textiles for complex real-world environments remains a challenge due to the limitations of functional integration, material properties, and preparation procedures. In this study, an aramid fabric based on in situ anchored copper sulfide nanostructure is developed. The textile with excellent solar and Joule heating properties can effectively keep the body warm even at low energy inputs. Meanwhile, the reduced infrared emissivity of the textile decreases the thermal radiation losses and helps to maintain a constant body temperature. Impressively, the textile integrates superb electromagnetic shielding, near-complete UV protection properties, and ideal resistance to fire and bacteria. This work provides a simple strategy for fabricating multi-functional integrated wearable devices with flexibility and breathability, which is highly promising in versatile PTM applications.

Identification of a cinnamoyl-CoA reductase from Cinnamomum cassia involved in trans-cinnamaldehyde biosynthesis.

MAIN CONCLUSION: The identification of a functional cinnamoyl-CoA reductase enzyme from Cinnamomum cassia involved in trans-cinnamaldehyde biosynthesis offers the potential for enhancing trans-cinnamaldehyde production through genetic engineering. A significant accumulation of trans-cinnamaldehyde has been found in the bark tissues of C. cassia, used in traditional Chinese medicine. trans-Cinnamaldehyde exhibits various pharmacological properties such as anti-inflammatory, analgesic, and protection of the stomach and the digestive tract. However, further elucidation and characterization of the biosynthetic pathway for trans-cinnamaldehyde is required. In this study, we conducted an integrated analysis of trans-cinnamaldehyde accumulation profiles and transcriptomic data from five different C. cassia tissues to identify the genes involved in its biosynthesis. The transcriptome data we obtained included nearly all genes associated with the trans-cinnamaldehyde pathway, with the majority demonstrating high abundance in branch barks and trunk barks. We successfully cloned four C. cassia cinnamoyl-CoA reductases (CcCCRs), a key gene in trans-cinnamaldehyde biosynthesis. We found that the recombinant CcCCR1 protein was the only one that more efficiently converted cinnamoyl-CoA into trans-cinnamaldehyde. CcCCR1 exhibited approximately 14.7-fold higher catalytic efficiency (kcat/Km) compared to the Arabidopsis thaliana cinnamoyl-CoA reductase 1 (AtCCR1); therefore, it can be utilized for engineering higher trans-cinnamaldehyde production as previously reported. Molecular docking studies and mutagenesis experiments also validated the superior catalytic activity of CcCCR1 compared to AtCCR1. These findings provide valuable insights for the functional characterization of enzyme-coding genes and hold potential for future engineering of trans-cinnamaldehyde biosynthetic pathways.

Distinct motifs in the E protein are required for SARS-CoV-2 virus particle formation and lysosomal deacidification in host cells.

IMPORTANCE: Severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2), the virus responsible for coronavirus disease 2019 (COVID-19), has caused a global public health crisis. The E protein, a structural protein found in this virus particle, is also known to be a viroporin. As such, it forms oligomeric ion channels or pores in the host cell membrane. However, the relationship between these two functions is poorly understood. In this study, we showed that the roles of E protein in virus particle and viroporin formation are distinct. This study contributes to the development of drugs that inhibit SARS-CoV-2 virus particle formation. Additionally, we designed a highly sensitive and high-throughput virus-like particle detection system using the HiBiT tag, which is a useful tool for studying the release of SARS-CoV-2.

Predictive, preventive, and personalized medicine in breast cancer: targeting the PI3K pathway.

Breast cancer (BC) is a multifaceted disease characterized by distinct molecular subtypes and varying responses to treatment. In BC, the phosphatidylinositol 3-kinase (PI3K) pathway has emerged as a crucial contributor to the development, advancement, and resistance to treatment. This review article explores the implications of the PI3K pathway in predictive, preventive, and personalized medicine for BC. It emphasizes the identification of predictive biomarkers, such as PIK3CA mutations, and the utility of molecular profiling in guiding treatment decisions. The review also discusses the potential of targeting the PI3K pathway for preventive strategies and the customization of therapy based on tumor stage, molecular subtypes, and genetic alterations. Overcoming resistance to PI3K inhibitors and exploring combination therapies are addressed as important considerations. While this field holds promise in improving patient outcomes, further research and clinical trials are needed to validate these approaches and translate them into clinical practice.

ZmSPL10, ZmSPL14 and ZmSPL26 act together to promote stigmatic papilla formation in maize through regulating auxin signaling and ZmWOX3A expression.

Maize silk is a specialized type of stigma, covered with numerous papillae for pollen grain capture. However, the developmental process of stigmatic papillae and the underlying regulatory mechanisms have remained largely unknown. Here, we combined the cytological, genetic and molecular studies to demonstrate that three homologous genes ZmSPL10, ZmSPL14 and ZmSPL26 play a central role in promoting stigmatic papilla formation in maize. We show that their triple knockout mutants are nearly complete lack of stigmatic papilla, resulting in a severe reduction in kernel setting. Cellular examination reveals that stigmatic papilla is developed from a precursor cell, which is the smaller daughter cell resulting from asymmetric cell division of a silk epidermal cell. In situ hybridization shows that ZmSPL10, ZmSPL14 and their target genes SPI1, ZmPIN1b, ZmARF28 and ZmWOX3A are preferentially expressed in the precursor cells of stigmatic papillae. Moreover, ZmSPL10, ZmSPL14 and ZmSPL26 directly bind to the promoters of SPI1, ZmPIN1b, ZmARF28 and ZmWOX3A and promote their expression. Further, Zmwox3a knockout mutants display severe defects in stigmatic papilla formation and reduced seed setting. Collectively, our results demonstrate that ZmSPL10, ZmSPL14 and ZmSPL26 act together to promote stigmatic papilla development through regulating auxin signaling and ZmWOX3A expression.