NAC transcription factors (TFs) participate in multiple biological processes, including biotic and abiotic stress responses, signal transduction and development. Cold stress can adversely impact plant growth and development, thereby limiting agricultural productivity. Prunus mume, an excellent horticultural crop, is widely cultivated in Asian countries. Its flower can tolerate freezing-stress in the early spring. To investigate the putative NAC genes responsible for cold-stress, we identified and analyzed 113 high-confidence PmNAC genes and characterized them by bioinformatics tools and expression profiles. These PmNACs were clustered into 14 sub-families and distributed on eight chromosomes and scaffolds, with the highest number located on chromosome 3. Duplicated events resulted in a large gene family; 15 and 8 pairs of PmNACs were the result of tandem and segmental duplicates, respectively. Moreover, three membrane-bound proteins (PmNAC59/66/73) and three miRNA-targeted genes (PmNAC40/41/83) were identified. Most PmNAC genes presented tissue-specific and time-specific expression patterns. Sixteen PmNACs (PmNAC11/19/20/23/41/48/58/74/75/76/78/79/85/86/103/111) exhibited down-regulation during flower bud opening and are, therefore, putative candidates for dormancy and cold-tolerance. Seventeen genes (PmNAC11/12/17/21/29/42/30/48/59/66/73/75/85/86/93/99/111) were highly expressed in stem during winter and are putative candidates for freezing resistance. The cold-stress response pattern of 15 putative PmNACs was observed under 4 °C at different treatment times. The expression of 10 genes (PmNAC11/20/23/40/42/48/57/60/66/86) was upregulated, while 5 genes (PmNAC59/61/82/85/107) were significantly inhibited. The putative candidates, thus identified, have the potential for breeding the cold-tolerant horticultural plants. This study increases our understanding of functions of the NAC gene family in cold tolerance, thereby potentially intensifying the molecular breeding programs of woody plants.

Summary: Plant with naturally twisted branches is referred to as a tortuous‐branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome‐level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi‐C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi‐omics in deciphering the properties of P. mume plant architecture.

Background: Mei (Prunus mume Sieb. et Zucc.) is a famous ornamental plant and fruit crop grown in East Asian countries. Limited genetic resources, especially molecular markers, have hindered the progress of mei breeding projects. Here, we performed low-depth whole-genome sequencing of Prunus mume ‘Fenban’ and Prunus mume ‘Kouzi Yudie’ to identify high-quality polymorphic markers between the two cultivars on a large scale. Results: A total of 1464.1 Mb and 1422.1 Mb of ‘Fenban’ and ‘Kouzi Yudie’ sequencing data were uniquely mapped to the mei reference genome with about 6-fold coverage, respectively. We detected a large number of putative polymorphic markers from the 196.9 Mb of sequencing data shared by the two cultivars, which together contained 200,627 SNPs, 4,900 InDels, and 7,063 SSRs. Among these markers, 38,773 SNPs, 174 InDels, and 418 SSRs were distributed in the 22.4 Mb CDS region, and 63.0% of these marker-containing CDS sequences were assigned to GO terms. Subsequently, 670 selected SNPs were validated using an Agilent’s SureSelect solution phase hybridization assay. A subset of 599 SNPs was used to assess the genetic similarity of a panel of mei germplasm samples and a plum (P. salicina) cultivar, producing a set of informative diversity data. We also analyzed the frequency and distribution of detected InDels and SSRs in mei genome and validated their usefulness as DNA markers. These markers were successfully amplified in the cultivars and in their segregating progeny. Conclusions: A large set of high-quality polymorphic SNPs, InDels, and SSRs were identified in parallel between ‘Fenban’ and ‘Kouzi Yudie’ using low-depth whole-genome sequencing. The study presents extensive data on these polymorphic markers, which can be useful for constructing high-resolution genetic maps, performing genome-wide association studies, and designing genomic selection strategies in mei.

Mei, Prunus mume Sieb. et Zucc., is an ornamental plant popular in East Asia and, as an important member of genus Prunus, has played a pivotal role in systematic studies of the Rosaceae. However, the genetic architecture of botanical traits in this species remains elusive. This paper represents the first genome-wide mapping study of quantitative trait loci (QTLs) that affect stem growth and form, leaf morphology and leaf anatomy in an intraspecific cross derived from two different mei cultivars. Genetic mapping based on a high-density linkage map constricted from 120 SSRs and 1,484 SNPs led to the detection of multiple QTLs for each trait, some of which exert pleiotropic effects on correlative traits. Each QTL explains 3-12% of the phenotypic variance. Several leaf size traits were found to share common QTLs, whereas growth-related traits and plant form traits might be controlled by a different set of QTLs. Our findings provide unique insights into the genetic control of tree growth and architecture in mei and help to develop an efficient breeding program for selecting superior mei cultivars.

Protein kinases of the MAPK cascade family (MAPKKK–MAPKK–MAPK) play an essential role in plant stress response and hormone signal transduction. However, their role in the cold hardiness of Prunus mume (Mei), a class of ornamental woody plant, remains unclear. In this study, we use bioinformatic approaches to assess and analyze two related protein kinase families, namely, MAP kinases (MPKs) and MAPK kinases (MKKs), in wild P. mume and its variety P. mume var. tortuosa. We identify 11 PmMPK and 7 PmMKK genes in the former species and 12 PmvMPK and 7 PmvMKK genes in the latter species, and we investigate whether and how these gene families contribute to cold stress responses. Members of the MPK and MKK gene families located on seven and four chromosomes of both species are free of tandem duplication. Four, three, and one segment duplication events are exhibited in PmMPK, PmvMPK, and PmMKK, respectively, suggesting that segment duplications play an essential role in the expansion and evolution of P. mume and its gene variety. Moreover, synteny analysis suggests that most MPK and MKK genes have similar origins and involved similar evolutionary processes in P. mume and its variety. A cis-acting regulatory element analysis shows that MPK and MKK genes may function in P. mume and its variety’s development, modulating processes such as light response, anaerobic induction, and abscisic acid response as well as responses to a variety of stresses, such as low temperature and drought. Most PmMPKs and PmMKKs exhibited tissue-specifific expression patterns, as well as time-specific expression patterns that protect them through cold. In a low-temperature treatment experiment with the cold-tolerant cultivar P. mume ‘Songchun’ and the cold-sensitive cultivar ‘Lve’, we find that almost all PmMPK and PmMKK genes, especially PmMPK3/5/6/20 and PmMKK2/3/6, dramatically respond to cold stress as treatment duration increases. This study introduces the possibility that these family members contribute to P. mume’s cold stress response. Further investigation is warranted to understand the mechanistic functions of MAPK and MAPKK proteins in P. mume development and response to cold stress.

High risk neuroblastoma (HR-NB) remains difficult to treat, and its overall survival (OS) is still below 50%. Although HR-NB is a heterogeneous disease, HR-NB patients are currently treated in a similar fashion. Through unsupervised biclustering, we further stratified HR-NB patients into two reproducible and clinically distinct subtypes, including an ultra-high risk neuroblastoma (UHR-NB) and high risk neuroblastoma (HR-NB). The UHR-NB subtype consistently had the worst OS in multiple independent cohorts ( $p<0.008$ ). Out of 283 neuroblastoma-specific immune genes that were used for stratification, 39 of them were differentiated in UHR-NB, including four upregulated and 35 downregulated, as compared to HR-NB. The four UHR-NB upregulated genes (ADAM22, GAL, KLHL13 and TWIST1) were all upregulated in MYCN amplified neuroblastoma in 5 additional cohorts. TWIST1 and ADAM22 were also positively correlated with cancer stage, while GAL was an independent OS predictor in addition to MYCN and age. Furthermore, we identified 26 commonly upregulated and 311 downregulated genes in UHR-NB from all 4723 immune-related genes. While 43 KEGG pathways with molecular functions were enriched in the downregulated immune-related genes, only the P53 signaling pathway was enriched in the upregulated ones, which suggested that UHR-NB was a TP53 related subtype with reduced immune activities.