Jurassic shuotheriids show earliest dental diversification of mammaliaforms.

Shuotheriids are Jurassic mammaliaforms that possess pseudotribosphenic teeth in which a pseudotalonid is anterior to the trigonid in the lower molar, contrasting with the tribosphenic pattern of therian mammals (placentals, marsupials and kin) in which the talonid is posterior to the trigonid1-4. The origin of the pseudotribosphenic teeth remains unclear, obscuring our perception of shuotheriid affinities and the early evolution of mammaliaforms1,5-9. Here we report a new Jurassic shuotheriid represented by two skeletal specimens. Their complete pseudotribosphenic dentitions allow reidentification of dental structures using serial homology and the tooth occlusal relationship. Contrary to the conventional view1,2,6,10,11, our findings show that dental structures of shuotheriids can be homologized to those of docodontans and partly support homologous statements for some dental structures between docodontans and other mammaliaforms6,12. The phylogenetic analysis based on new evidence removes shuotheriids from the tribosphenic ausktribosphenids (including monotremes) and clusters them with docodontans to form a new clade, Docodontiformes, that is characterized by pseudotribosphenic features. In the phylogeny, docodontiforms and 'holotherians' (Kuehneotherium, monotremes and therians)13 evolve independently from a Morganucodon-like ancestor with triconodont molars by labio-lingual widening their posterior teeth for more efficient food processing. The pseudotribosphenic pattern passed a cusp semitriangulation stage9, whereas the tribosphenic pattern and its precursor went through a stage of cusp triangulation. The two different processes resulted in complex tooth structures and occlusal patterns that elucidate the earliest diversification of mammaliaforms.

Fossils document evolutionary changes of jaw joint to mammalian middle ear.

The dual jaw joint of Morganucodon1,2 consists of the dentary-squamosal joint laterally and the articular-quadrate one medially. The articular-quadrate joint and its associated post-dentary bones constitute the precursor of the mammalian middle ear. Fossils documenting the transition from such a precursor to the mammalian middle ear are poor, resulting in inconsistent interpretations of this hallmark apparatus in the earliest stage of mammaliaform evolution1-5. Here we report mandibular middle ears from two Jurassic mammaliaforms: a new morganucodontan-like species and a pseudotribosphenic shuotheriid species6. The morganucodontan-like species shows many previously unknown post-dentary bone morphologies1,2 and exhibits features that suggest a loss of load-bearing function in its articular-quadrate joint. The middle ear of the shuotheriid approaches the mammalian condition in that it has features that are suitable for an exclusively auditory function, although the post-dentary bones are still attached to the dentary. With size reduction of the jaw-joint bones, the quadrate shifts medially at different degrees in relation to the articular in the two mammaliaforms. These changes provide evidence of a gradual loss of load-bearing function in the articular-quadrate jaw joint-a prerequisite for the detachment of the post-dentary bones from the dentary7-12 and the eventual breakdown of the Meckel's cartilage13-15 during the evolution of mammaliaforms.

Antiferromagnetic phase transition in a 3D fermionic Hubbard model.

The fermionic Hubbard model (FHM)1 describes a wide range of physical phenomena resulting from strong electron-electron correlations, including conjectured mechanisms for unconventional superconductivity. Resolving its low-temperature physics is, however, challenging theoretically or numerically. Ultracold fermions in optical lattices2,3 provide a clean and well-controlled platform offering a path to simulate the FHM. Doping the antiferromagnetic ground state of a FHM simulator at half-filling is expected to yield various exotic phases, including stripe order4, pseudogap5, and d-wave superfluid6, offering valuable insights into high-temperature superconductivity7-9. Although the observation of antiferromagnetic correlations over short10 and extended distances11 has been obtained, the antiferromagnetic phase has yet to be realized as it requires sufficiently low temperatures in a large and uniform quantum simulator. Here we report the observation of the antiferromagnetic phase transition in a three-dimensional fermionic Hubbard system comprising lithium-6 atoms in a uniform optical lattice with approximately 800,000 sites. When the interaction strength, temperature and doping concentration are finely tuned to approach their respective critical values, a sharp increase in the spin structure factor is observed. These observations can be well described by a power-law divergence, with a critical exponent of 1.396 from the Heisenberg universality class12. At half-filling and with optimal interaction strength, the measured spin structure factor reaches 123(8), signifying the establishment of an antiferromagnetic phase. Our results provide opportunities for exploring the low-temperature phase diagram of the FHM.

Targeting 7-Dehydrocholesterol Reductase Integrates Cholesterol Metabolism and IRF3 Activation to Eliminate Infection.

Recent work suggests that cholesterol metabolism impacts innate immune responses against infection. However, the key enzymes or the natural products and mechanisms involved are not well elucidated. Here, we have shown that upon DNA and RNA viral infection, macrophages reduced 7-dehydrocholesterol reductase (DHCR7) expression. DHCR7 deficiency or treatment with the natural product 7-dehydrocholesterol (7-DHC) could specifically promote phosphorylation of IRF3 (not TBK1) and enhance type I interferon (IFN-I) production in macrophages. We further elucidated that viral infection or 7-DHC treatment enhanced AKT3 expression and activation. AKT3 directly bound and phosphorylated IRF3 at Ser385, together with TBK1-induced phosphorylation of IRF3 Ser386, to achieve IRF3 dimerization. Deletion of DHCR7 and the DHCR7 inhibitors including AY9944 and the chemotherapy drug tamoxifen promoted clearance of Zika virus and multiple viruses in vitro or in vivo. Taken together, we propose that the DHCR7 inhibitors and 7-DHC are potential therapeutics against emerging or highly pathogenic viruses.

A new avialan theropod from an emerging Jurassic terrestrial fauna.

Birds are descended from non-avialan theropod dinosaurs of the Late Jurassic period, but the earliest phase of this evolutionary process remains unclear owing to the exceedingly sparse and spatio-temporally restricted fossil record1-5. Information about the early-diverging species along the avialan line is crucial to understand the evolution of the characteristic bird bauplan, and to reconcile phylogenetic controversies over the origin of birds3,4. Here we describe one of the stratigraphically youngest and geographically southernmost Jurassic avialans, Fujianvenator prodigiosus gen. et sp. nov., from the Tithonian age of China. This specimen exhibits an unusual set of morphological features that are shared with other stem avialans, troodontids and dromaeosaurids, showing the effects of evolutionary mosaicism in deep avialan phylogeny. F. prodigiosus is distinct from all other Mesozoic avialan and non-avialan theropods in having a particularly elongated hindlimb, suggestive of a terrestrial or wading lifestyle-in contrast with other early avialans, which exhibit morphological adaptations to arboreal or aerial environments. During our fieldwork in Zhenghe where F. prodigiosus was found, we discovered a diverse assemblage of vertebrates dominated by aquatic and semi-aquatic species, including teleosts, testudines and choristoderes. Using in situ radioisotopic dating and stratigraphic surveys, we were able to date the fossil-containing horizons in this locality-which we name the Zhenghe Fauna-to 148-150 million years ago. The diversity of the Zhenghe Fauna and its precise chronological framework will provide key insights into terrestrial ecosystems of the Late Jurassic.

Fossoriality and evolutionary development in two Cretaceous mammaliamorphs.

Mammaliamorpha comprises the last common ancestor of Tritylodontidae and Mammalia plus all its descendants1. Tritylodontids are nonmammaliaform herbivorous cynodonts that originated in the Late Triassic epoch, diversified in the Jurassic period2-5 and survived into the Early Cretaceous epoch6,7. Eutriconodontans have generally been considered to be an extinct mammalian group, although different views exist8. Here we report a newly discovered tritylodontid and eutriconodontan from the Early Cretaceous Jehol Biota of China. Eutriconodontans are common in this biota9, but it was not previously known to contain tritylodontids. The two distantly related species show convergent features that are adapted for fossorial life, and are the first 'scratch-diggers' known from this biota. Both species also show an increased number of presacral vertebrae, relative to the ancestral state in synapsids or mammals10,11, that display meristic and homeotic changes. These fossils shed light on the evolutionary development of the axial skeleton in mammaliamorphs, which has been the focus of numerous studies in vertebrate evolution12-17 and developmental biology18-28. The phenotypes recorded by these fossils indicate that developmental plasticity in somitogenesis and HOX gene expression in the axial skeleton-similar to that observed in extant mammals-was already in place in stem mammaliamorphs. The interaction of these developmental mechanisms with natural selection may have underpinned the diverse phenotypes of body plan that evolved independently in various clades of mammaliamorph.

Non-KREEP origin for Chang'e-5 basalts in the Procellarum KREEP Terrane.

Mare volcanics on the Moon are the key record of thermo-chemical evolution throughout most of lunar history1-3. Young mare basalts-mainly distributed in a region rich in potassium, rare-earth elements and phosphorus (KREEP) in Oceanus Procellarum, called the Procellarum KREEP Terrane (PKT)4-were thought to be formed from KREEP-rich sources at depth5-7. However, this hypothesis has not been tested with young basalts from the PKT. Here we present a petrological and geochemical study of the basalt clasts from the PKT returned by the Chang'e-5 mission8. These two-billion-year-old basalts are the youngest lunar samples reported so far9. Bulk rock compositions have moderate titanium and high iron contents  with KREEP-like rare-earth-element and high thorium concentrations. However, strontium-neodymium isotopes indicate that these basalts were derived from a non-KREEP mantle source. To produce the high abundances of rare-earth elements and thorium, low-degree partial melting and extensive fractional crystallization are required. Our results indicate that the KREEP association may not be a prerequisite for young mare volcanism. Absolving the need to invoke heat-producing elements in their source implies a more sustained cooling history of the lunar interior to generate the Moon's youngest melts.

A highly reproducible and efficient method for retinal organoid differentiation from human pluripotent stem cells.

Human pluripotent stem cell (hPSC)-derived retinal organoids are three-dimensional cellular aggregates that differentiate and self-organize to closely mimic the spatial and temporal patterning of the developing human retina. Retinal organoid models serve as reliable tools for studying human retinogenesis, yet limitations in the efficiency and reproducibility of current retinal organoid differentiation protocols have reduced the use of these models for more high-throughput applications such as disease modeling and drug screening. To address these shortcomings, the current study aimed to standardize prior differentiation protocols to yield a highly reproducible and efficient method for generating retinal organoids. Results demonstrated that through regulation of organoid size and shape using quick reaggregation methods, retinal organoids were highly reproducible compared to more traditional methods. Additionally, the timed activation of BMP signaling within developing cells generated pure populations of retinal organoids at 100% efficiency from multiple widely used cell lines, with the default forebrain fate resulting from the inhibition of BMP signaling. Furthermore, given the ability to direct retinal or forebrain fates at complete purity, mRNA-seq analyses were then utilized to identify some of the earliest transcriptional changes that occur during the specification of these two lineages from a common progenitor. These improved methods also yielded retinal organoids with expedited differentiation timelines when compared to traditional methods. Taken together, the results of this study demonstrate the development of a highly reproducible and minimally variable method for generating retinal organoids suitable for analyzing the earliest stages of human retinal cell fate specification.

AKT2 reduces IFNß1 production to modulate antiviral responses and systemic lupus erythematosus.

Interferon regulatory factor 3 (IRF3)-induced type I interferon (I-IFN) production plays key roles in both antiviral and autoimmune responses. IRF3 phosphorylation, dimerization, and nuclear localization are needed for its activation and function, but the precise regulatory mechanisms remain to be explored. Here, we show that the serine/threonine kinase AKT2 interacts with IRF3 and phosphorylates it on Thr207, thereby attenuating IRF3 nuclear translocation in a 14-3-3ε-dependent manner and reducing I-IFN production. We further find that AKT2 expression is downregulated in viral-infected macrophages or in monocytes and tissue samples from systemic lupus erythematosus (SLE) patients and mouse models. Akt2-deficient mice exhibit increased I-IFN induction and reduced mortality in response to viral infection, but aggravated severity of SLE. Overexpression of AKT2 kinase-inactive or IRF3-T207A mutants in zebrafish supports that AKT2 negatively regulates I-IFN production and antiviral response in a kinase-dependent manner. This negative role of AKT2 in IRF3-induced I-IFN production suggests that AKT2 may be therapeutically targeted to differentially regulate antiviral infection and SLE.

Whole-genome sequencing of diverse wheat accessions uncovers genetic changes during modern breeding in China and the United States.

Breeding has dramatically changed the plant architecture of wheat (Triticum aestivum), resulting in the development of high-yielding varieties adapted to modern farming systems. However, how wheat breeding shaped the genomic architecture of this crop remains poorly understood. Here, we performed a comprehensive comparative analysis of a whole-genome resequencing panel of 355 common wheat accessions (representing diverse landraces and modern cultivars from China and the United States) at the phenotypic and genomic levels. The genetic diversity of modern wheat cultivars was clearly reduced compared to landraces. Consistent with these genetic changes, most phenotypes of cultivars from China and the United States were significantly altered. Of the 21 agronomic traits investigated, 8 showed convergent changes between the 2 countries. Moreover, of the 207 loci associated with these 21 traits, more than half overlapped with genomic regions that showed evidence of selection. The distribution of selected loci between the Chinese and American cultivars suggests that breeding for increased productivity in these 2 regions was accomplished by pyramiding both shared and region-specific variants. This work provides a framework to understand the genetic architecture of the adaptation of wheat to diverse agricultural production environments, as well as guidelines for optimizing breeding strategies to design better wheat varieties.

Integrated optical multi-ion quantum logic.

Practical and useful quantum information processing requires substantial improvements with respect to current systems, both in the error rates of basic operations and in scale. The fundamental qualities of individual trapped-ion1 qubits are promising for long-term systems2, but the optics involved in their precise control are a barrier to scaling3. Planar-fabricated optics integrated within ion-trap devices can make such systems simultaneously more robust and parallelizable, as suggested by previous work with single ions4. Here we use scalable optics co-fabricated with a surface-electrode ion trap to achieve high-fidelity multi-ion quantum logic gates, which are often the limiting elements in building up the precise, large-scale entanglement that is essential to quantum computation. Light is efficiently delivered to a trap chip in a cryogenic environment via direct fibre coupling on multiple channels, eliminating the need for beam alignment into vacuum systems and cryostats and lending robustness to vibrations and beam-pointing drifts. This allows us to perform ground-state laser cooling of ion motion and to implement gates generating two-ion entangled states with fidelities greater than 99.3(2) per cent. This work demonstrates hardware that reduces noise and drifts in sensitive quantum logic, and simultaneously offers a route to practical parallelization for high-fidelity quantum processors5. Similar devices may also find applications in atom- and ion-based quantum sensing and timekeeping6.

Submicrosecond entangling gate between trapped ions via Rydberg interaction.

Generating quantum entanglement in large systems on timescales much shorter than the coherence time is key to powerful quantum simulation and computation. Trapped ions are among the most accurately controlled and best isolated quantum systems1 with low-error entanglement gates operated within tens of microseconds using the vibrational motion of few-ion crystals2,3. To exceed the level of complexity tractable by classical computers the main challenge is to realize fast entanglement operations in crystals made up of many ions (large ion crystals)4. The strong dipole-dipole interactions in polar molecule5 and Rydberg atom6,7 systems allow much faster entangling gates, yet stable state-independent confinement comparable with trapped ions needs to be demonstrated in these systems8. Here we combine the benefits of these approaches: we report a two-ion entangling gate with 700-nanosecond gate time that uses the strong dipolar interaction between trapped Rydberg ions, which we use to produce a Bell state with 78 per cent fidelity. The sources of gate error are identified and a total error of less than 0.2 per cent is predicted for experimentally achievable parameters. Furthermore, we predict that residual coupling to motional modes contributes an approximate gate error of 10-4 in a large ion crystal of 100 ions. This provides a way to speed up and scale up trapped-ion quantum computers and simulators substantially.

Effect of non-pharmaceutical interventions to contain COVID-19 in China.

On 11 March 2020, the World Health Organization (WHO) declared coronavirus disease 2019 (COVID-19) a pandemic1. The strategies based on non-pharmaceutical interventions that were used to contain the outbreak in China appear to be effective2, but quantitative research is still needed to assess the efficacy of non-pharmaceutical interventions and their timings3. Here, using epidemiological data on COVID-19 and anonymized data on human movement4,5, we develop a modelling framework that uses daily travel networks to simulate different outbreak and intervention scenarios across China. We estimate that there were a total of 114,325 cases of COVID-19 (interquartile range 76,776-164,576) in mainland China as of 29 February 2020. Without non-pharmaceutical interventions, we predict that the number of cases would have been 67-fold higher (interquartile range 44-94-fold) by 29 February 2020, and we find that the effectiveness of different interventions varied. We estimate that early detection and isolation of cases prevented more infections than did travel restrictions and contact reductions, but that a combination of non-pharmaceutical interventions achieved the strongest and most rapid effect. According to our model, the lifting of travel restrictions from 17 February 2020 does not lead to an increase in cases across China if social distancing interventions can be maintained, even at a limited level of an on average 25% reduction in contact between individuals that continues until late April. These findings improve our understanding of the effects of non-pharmaceutical interventions on COVID-19, and will inform response efforts across the world.

The water lily genome and the early evolution of flowering plants.

Water lilies belong to the angiosperm order Nymphaeales. Amborellales, Nymphaeales and Austrobaileyales together form the so-called ANA-grade of angiosperms, which are extant representatives of lineages that diverged the earliest from the lineage leading to the extant mesangiosperms1-3. Here we report the 409-megabase genome sequence of the blue-petal water lily (Nymphaea colorata). Our phylogenomic analyses support Amborellales and Nymphaeales as successive sister lineages to all other extant angiosperms. The N. colorata genome and 19 other water lily transcriptomes reveal a Nymphaealean whole-genome duplication event, which is shared by Nymphaeaceae and possibly Cabombaceae. Among the genes retained from this whole-genome duplication are homologues of genes that regulate flowering transition and flower development. The broad expression of homologues of floral ABCE genes in N. colorata might support a similarly broadly active ancestral ABCE model of floral organ determination in early angiosperms. Water lilies have evolved attractive floral scents and colours, which are features shared with mesangiosperms, and we identified their putative biosynthetic genes in N. colorata. The chemical compounds and biosynthetic genes behind floral scents suggest that they have evolved in parallel to those in mesangiosperms. Because of its unique phylogenetic position, the N. colorata genome sheds light on the early evolution of angiosperms.

JANUS, a spliceosome-associated protein, promotes miRNA biogenesis in Arabidopsis.

MicroRNAs (miRNAs) are important regulators of genes expression. Their levels are precisely controlled through modulating the activity of the microprocesser complex (MC). Here, we report that JANUS, a homology of the conserved U2 snRNP assembly factor in yeast and human, is required for miRNA accumulation. JANUS associates with MC components Dicer-like 1 (DCL1) and SERRATE (SE) and directly binds the stem-loop of pri-miRNAs. In a hypomorphic janus mutant, the activity of DCL1, the numbers of MC, and the interaction of primary miRNA transcript (pri-miRNAs) with MC are reduced. These data suggest that JANUS promotes the assembly and activity of MC through its interaction with MC and/or pri-miRNAs. In addition, JANUS modulates the transcription of some pri-miRNAs as it binds the promoter of pri-miRNAs and facilitates Pol II occupancy of at their promoters. Moreover, global splicing defects are detected in janus. Taken together, our study reveals a novel role of a conserved splicing factor in miRNA biogenesis.

GPS-SUMO 2.0: an updated online service for the prediction of SUMOylation sites and SUMO-interacting motifs.

Small ubiquitin-like modifiers (SUMOs) are tiny but important protein regulators involved in orchestrating a broad spectrum of biological processes, either by covalently modifying protein substrates or by noncovalently interacting with other proteins. Here, we report an updated server, GPS-SUMO 2.0, for the prediction of SUMOylation sites and SUMO-interacting motifs (SIMs). For predictor training, we adopted three machine learning algorithms, penalized logistic regression (PLR), a deep neural network (DNN), and a transformer, and used 52 404 nonredundant SUMOylation sites in 8262 proteins and 163 SIMs in 102 proteins. To further increase the accuracy of predicting SUMOylation sites, a pretraining model was first constructed using 145 545 protein lysine modification sites, followed by transfer learning to fine-tune the model. GPS-SUMO 2.0 exhibited greater accuracy in predicting SUMOylation sites than did other existing tools. For users, one or multiple protein sequences or identifiers can be input, and the prediction results are shown in a tabular list. In addition to the basic statistics, we integrated knowledge from 35 public resources to annotate SUMOylation sites or SIMs. The GPS-SUMO 2.0 server is freely available at https://sumo.biocuckoo.cn/. We believe that GPS-SUMO 2.0 can serve as a useful tool for further analysis of SUMOylation and SUMO interactions.

Establishing Wiring Specificity in Visual System Circuits: From the Retina to the Brain.

The retina is a tremendously complex image processor, containing numerous cell types that form microcircuits encoding different aspects of the visual scene. Each microcircuit exhibits a distinct pattern of synaptic connectivity. The developmental mechanisms responsible for this patterning are just beginning to be revealed. Furthermore, signals processed by different retinal circuits are relayed to specific, often distinct, brain regions. Thus, much work has focused on understanding the mechanisms that wire retinal axonal projections to their appropriate central targets. Here, we highlight recently discovered cellular and molecular mechanisms that together shape stereotypic wiring patterns along the visual pathway, from within the retina to the brain. Although some mechanisms are common across circuits, others play unconventional and circuit-specific roles. Indeed, the highly organized connectivity of the visual system has greatly facilitated the discovery of novel mechanisms that establish precise synaptic connections within the nervous system.

Spatiotemporal genotype replacement of H5N8 avian influenza viruses contributed to H5N1 emergence in 2021/2022 panzootic.

Since 2020, clade 2.3.4.4b highly pathogenic avian influenza H5N8 and H5N1 viruses have swept through continents, posing serious threats to the world. Through comprehensive analyses of epidemiological, genetic, and bird migration data, we found that the dominant genotype replacement of the H5N8 viruses in 2020 contributed to the H5N1 outbreak in the 2021/2022 wave. The 2020 outbreak of the H5N8 G1 genotype instead of the G0 genotype produced reassortment opportunities and led to the emergence of a new H5N1 virus with G1's HA and MP genes. Despite extensive reassortments in the 2021/2022 wave, the H5N1 virus retained the HA and MP genes, causing a significant outbreak in Europe and North America. Furtherly, through the wild bird migration flyways investigation, we found that the temporal-spatial coincidence between the outbreak of the H5N8 G1 virus and the bird autumn migration may have expanded the H5 viral spread, which may be one of the main drivers of the emergence of the 2020-2022 H5 panzootic.IMPORTANCESince 2020, highly pathogenic avian influenza (HPAI) H5 subtype variants of clade 2.3.4.4b have spread across continents, posing unprecedented threats globally. However, the factors promoting the genesis and spread of H5 HPAI viruses remain unclear. Here, we found that the spatiotemporal genotype replacement of H5N8 HPAI viruses contributed to the emergence of the H5N1 variant that caused the 2021/2022 panzootic, and the viral evolution in poultry of Egypt and surrounding area and autumn bird migration from the Russia-Kazakhstan region to Europe are important drivers of the emergence of the 2020-2022 H5 panzootic. These findings provide important targets for early warning and could help control the current and future HPAI epidemics.

Benchmarking of local genetic correlation estimation methods using summary statistics from genome-wide association studies.

Local genetic correlation evaluates the correlation of additive genetic effects between different traits across the same genetic variants at a genomic locus. It has been proven informative for understanding the genetic similarities of complex traits beyond that captured by global genetic correlation calculated across the whole genome. Several summary-statistics-based approaches have been developed for estimating local genetic correlation, including $\rho$-hess, SUPERGNOVA and LAVA. However, there has not been a comprehensive evaluation of these methods to offer practical guidelines on the choices of these methods. In this study, we conduct benchmark comparisons of the performance of these three methods through extensive simulation and real data analyses. We focus on two technical difficulties in estimating local genetic correlation: sample overlaps across traits and local linkage disequilibrium (LD) estimates when only the external reference panels are available. Our simulations suggest the likelihood of incorrectly identifying correlated regions and local correlation estimation accuracy are highly dependent on the estimation of the local LD matrix. These observations are corroborated by real data analyses of 31 complex traits. Overall, our findings illuminate the distinct results yielded by different methods applied in post-genome-wide association studies (post-GWAS) local correlation studies. We underscore the sensitivity of local genetic correlation estimates and inferences to the precision of local LD estimation. These observations accentuate the vital need for ongoing refinement in methodologies.

Transcription factor PpNAC1 and DNA demethylase PpDML1 synergistically regulate peach fruit ripening.

Fruit ripening is accompanied by dramatic changes in color, texture, and flavor and is regulated by transcription factors (TFs) and epigenetic factors. However, the detailed regulatory mechanism remains unclear. Gene expression patterns suggest that PpNAC1 (NAM/ATAF1/2/CUC) TF plays a major role in peach (Prunus persica) fruit ripening. DNA affinity purification (DAP)-seq combined with transactivation tests demonstrated that PpNAC1 can directly activate the expression of multiple ripening-related genes, including ACC synthase1 (PpACS1) and ACC oxidase1 (PpACO1) involved in ethylene biosynthesis, pectinesterase1 (PpPME1), pectate lyase1 (PpPL1), and polygalacturonase1 (PpPG1) related to cell wall modification, and lipase1 (PpLIP1), fatty acid desaturase (PpFAD3-1), and alcohol acyltransferase1 (PpAAT1) involved in volatiles synthesis. Overexpression of PpNAC1 in the tomato (Solanum lycopersicum) nor (nonripening) mutant restored fruit ripening, and its transient overexpression in peach fruit induced target gene expression, supporting a positive role of PpNAC1 in fruit ripening. The enhanced transcript levels of PpNAC1 and its target genes were associated with decreases in their promoter mCG methylation during ripening. Declining DNA methylation was negatively associated with increased transcripts of DNA demethylase1 (PpDML1), whose promoter is recognized and activated by PpNAC1. We propose that decreased methylation of the promoter region of PpNAC1 leads to a subsequent decrease in DNA methylation levels and enhanced transcription of ripening-related genes. These results indicate that positive feedback between PpNAC1 and PpDML1 plays an important role in directly regulating expression of multiple genes required for peach ripening and quality formation.