Summary: Single‐cell RNA‐seq (scRNA‐seq) has been highlighted as a powerful tool for the description of human cell transcriptome, but the technology has not been broadly applied in plant cells. Herein, we describe the successful development of a robust protoplast cell isolation system in the peanut leaf. A total of 6,815 single cells were divided into eight cell clusters based on reported marker genes by applying scRNA‐seq. Further, a pseudo‐time analysis was used to describe the developmental trajectory and interaction network of transcription factors (TFs) of distinct cell types during leaf growth. The trajectory enabled re‐investigation of the primordium‐driven development processes of the mesophyll and epidermis. These results suggest that palisade cells likely differentiate into spongy cells, while the epidermal cells originated earlier than the primordium. Subsequently, the developed method integrated multiple technologies to efficiently validate the scRNA‐seq result in a homogenous cell population. The expression levels of several TFs were strongly correlated with epidermal ontogeny in accordance with obtained scRNA‐seq values. Additionally, peanut AHL23 (AT‐HOOK MOTIF NUCLEAR LOCALIZED PROTEIN 23), which is localized in nucleus, promoted leaf growth when ectopically expressed in Arabidopsis by modulating the phytohormone pathway. Together, our study displays that application of scRNA‐seq can provide new hypotheses regarding cell differentiation in the leaf blade of Arachis hypogaea. We believe that this approach will enable significant advances in the functional study of leaf blade cells in the allotetraploid peanut and other plant species.

Many large-effect quantitative trait loci (QTLs) for yield and disease resistance related traits have been identified in different mapping populations of peanut (Arachis hypogaea L.) under multiple environments. However, only a limited number of QTLs have been used in marker-assisted selection (MAS) because of unfavorable epistatic interactions between QTLs in different genetic backgrounds. Thus, it is essential to identify consensus QTLs across different environments and genetic backgrounds for use in MAS. Here, we used QTL meta-analysis to identify a set of consensus QTLs for yield and disease resistance related traits in peanut. A new integrated consensus genetic map with 5874 loci was constructed. The map comprised 20 linkage groups (LGs) and was up to a total length of 2918.62 cM with average marker density of 2.01 loci per centimorgan (cM). A total of 292 initial QTLs were projected on the new consensus map, and 40 meta-QTLs (MQTLs) for yield and disease resistance related traits were detected on four LGs. The genetic intervals of these consensus MQTLs varied from 0.20 cM to 7.4 cM, which is narrower than the genetic intervals of the initial QTLs, meaning they may be suitable for use in MAS. Importantly, a region of the map that previously co-localized multiple major QTLs for pod traits was narrowed from 3.7 cM to 0.7 cM using an overlap region of four MQTLs for yield related traits on LG A05, which corresponds to a physical region of about 630.3 kb on the A05 pseudomolecule of peanut, including 38 annotated candidate genes (54 transcripts) related to catalytic activity and metabolic process. Additionally, one major MQTL for late leaf spot (LLS) was identified in a region of about 0.38 cM. BLAST searches identified 26 candidate genes (30 different transcripts) in this region, some of which were annotated as related to regulation of disease resistance in different plant species. Combined with the high-density marker consensus map, all the detected MQTLs could be useful in MAS. The biological functions of the 64 candidate genes should be validated to unravel the molecular mechanisms of yield and disease resistance in peanut.

Peanuts with high oleic acid content are usually considered to be beneficial for human health and edible oil storage. In breeding practice, peanut lines with high monounsaturated fatty acids are selected using fatty acid desaturase 2 (FAD2), which is responsible for the conversion of oleic acid (C18:1) to linoleic acid (C18:2). Here, comparative transcriptomics were used to analyze the global gene expression profile of high- and normal-oleic peanut cultivars at six time points during seed development. First, the mutant type of FAD2 was determined in the high-oleic peanut (H176). The result suggested that early translation termination occurred simultaneously in the coding sequence of FAD2-A and FAD2-B, and the cultivar H176 is capable of utilizing a potential germplasm resource for future high-oleic peanut breeding. Furthermore, transcriptomic analysis identified 74 differentially expressed genes (DEGs) involved in lipid metabolism in high-oleic peanut seed, of which five DEGs encoded the fatty acid desaturase. Aradu.XM2MR belonged to the homologous gene of stearoyl-ACP (acyl carrier protein) desaturase 2 (SAD2) that converted the C18:0 into C18:1. Further subcellular localization studies indicated that FAD2 was located at the endoplasmic reticulum (ER), and Aradu.XM2MR was targeted to the plastid in Arabidopsis protoplast cells. To examine the dynamic mechanism of this finding, we focused on the peroxidase (POD)-mediated fatty acid (FA) degradation pathway. The fad2 mutant significantly increased the POD activity and H₂O₂ concentration at the early stage of seed development, implying that redox signaling likely acted as a messenger to connect the signaling transduction between the high-oleic content and Aradu.XM2MR transcription level. Taken together, transcriptome analysis revealed the feedback mechanism of SAD2 (Aradu.XM2MR) associated with FAD2 mutation during the seed developmental stage, which could provide a potential peanut breeding strategy based on identified candidate genes to improve the content of oleic acid.

Peanuts (Arachis hypogaea L.) are an important oilseed crop, containing high contents of protein and fatty acids (FA). The major components of FA found in peanut oil are unsaturated FAs, including oleic acid (OA, C18:1) and linoleic acid (LOA, C18:2). Moreover, the high content of OA in peanut oil is beneficial for human health and long-term storage due to its antioxidant activity. However, the dynamic changes in proteomics related to OA accumulation during seed development still remain largely unexplored. In the present study, a comparative proteome analysis based on iTRAQ (isobaric Tags for Relative and Absolute Quantification) was performed to identify the critical candidate factors involved in OA formation. A total of 389 differentially expressed proteins (DEPs) were identified between high-oleate cultivar Kainong176 and low-oleate cultivar Kainong70. Among these DEPs, 201 and 188 proteins were upregulated and downregulated, respectively. In addition, these DEPs were categorized into biosynthesis pathways of unsaturated FAs at the early stage during the high-oleic peanut seed development, and several DEPs involved in lipid oxidation pathway were found at the stage of seed maturation. Meanwhile, 28 DEPs were sporadically distributed in distinct stages of seed formation, and their molecular functions were directly correlated to FA biosynthesis and degradation. Fortunately, the expression of FAB2 (stearoyl-acyl carrier protein desaturase), the rate-limiting enzyme in the upstream biosynthesis process of OA, was significantly increased in the early stage and then decreased in the late stage of seed development in the high-oleate cultivar Kainong176. Furthermore, real-time PCR verified the expression pattern of FAB2 at the mRNA level, which was consistent with its protein abundance. However, opposite results were found for the low-oleate cultivar Kainong70. Overall, the comparative proteome analysis provided valuable insight into the molecular dynamics of OA accumulation during peanut seed development.