ScRNA-seq reveals dark- and light-induced differentially expressed gene atlases of seedling leaves in Arachis hypogaea L.

Although the regulatory mechanisms of dark and light-induced plant morphogenesis have been broadly investigated, the biological process in peanuts has not been systematically explored on single-cell resolution. Herein, 10 cell clusters were characterized using scRNA-seq-identified marker genes, based on 13 409 and 11 296 single cells from 1-week-old peanut seedling leaves grown under dark and light conditions. 6104 genes and 50 transcription factors (TFs) displayed significant expression patterns in distinct cell clusters, which provided gene resources for profiling dark/light-induced candidate genes. Further pseudo-time trajectory and cell cycle evidence supported that dark repressed the cell division and perturbed normal cell cycle, especially the PORA abundances correlated with 11 TFs highly enriched in mesophyll to restrict the chlorophyllide synthesis. Additionally, light repressed the epidermis cell developmental trajectory extending by inhibiting the growth hormone pathway, and 21 TFs probably contributed to the different genes transcriptional dynamic. Eventually, peanut AHL17 was identified from the profile of differentially expressed TFs, which encoded protein located in the nucleus promoted leaf epidermal cell enlargement when ectopically overexpressed in Arabidopsis through the regulatory phytohormone pathway. Overall, our study presents the different gene atlases in peanut etiolated and green seedlings, providing novel biological insights to elucidate light-induced leaf cell development at the single-cell level.

Phosphatidyl Ethanolamine Binding Protein FLOWERING LOCUS T-like 12 (OsFTL12) Regulates the Rice Heading Date under Different Day-Length Conditions.

Plant FLOWERING LOCUS T-Like (FTL) genes often redundantly duplicate on chromosomes and functionally diverge to modulate reproductive traits. Rice harbors thirteen FTL genes, the functions of which are still not clear, except for the Hd3a and RFT genes. Here, we identified the molecular detail of OsFTL12 in rice reproductive stage. OsFTL12 encoding protein contained PEBP domain and localized into the nucleus, which transcripts specifically expressed in the shoot and leaf blade with high abundance. Further GUS-staining results show the OsFTL12 promoter activity highly expressed in the leaf and stem. OsFTL12 knock-out concurrently exhibited early flowering phenotype under the short- and long-day conditions as compared with wild-type and over-expression plants, which independently regulates flowering without an involved Hd1/Hd3a and Ehd1/RFT pathway. Further, an AT-hook protein OsATH1 was identified to act as upstream regulator of OsFTL12, as the knock-out OsATH1 elevated the OsFTL12 expression by modifying Histone H3 acetylation abundance. According to the dissection of OsFTL12 molecular functions, our study expanded the roles intellectual function of OsFTL12 in the mediating of a rice heading date.

Silicon Application for the Modulation of Rhizosphere Soil Bacterial Community Structures and Metabolite Profiles in Peanut under Ralstonia solanacearum Inoculation.

Silicon (Si) has been shown to promote peanut growth and yield, but whether Si can enhance the resistance against peanut bacterial wilt (PBW) caused by Ralstonia solanacearum, identified as a soil-borne pathogen, is still unclear. A question regarding whether Si enhances the resistance of PBW is still unclear. Here, an in vitro R. solanacearum inoculation experiment was conducted to study the effects of Si application on the disease severity and phenotype of peanuts, as well as the microbial ecology of the rhizosphere. Results revealed that Si treatment significantly reduced the disease rate, with a decrement PBW severity of 37.50% as compared to non-Si treatment. The soil available Si (ASi) significantly increased by 13.62-44.87%, and catalase activity improved by 3.01-3.10%, which displayed obvious discrimination between non-Si and Si treatments. Furthermore, the rhizosphere soil bacterial community structures and metabolite profiles dramatically changed under Si treatment. Three significantly changed bacterial taxa were observed, which showed significant abundance under Si treatment, whereas the genus Ralstonia genus was significantly suppressed by Si. Similarly, nine differential metabolites were identified to involve into unsaturated fatty acids via a biosynthesis pathway. Significant correlations were also displayed between soil physiochemical properties and enzymes, the bacterial community, and the differential metabolites by pairwise comparisons. Overall, this study reports that Si application mediated the evolution of soil physicochemical properties, the bacterial community, and metabolite profiles in the soil rhizosphere, which significantly affects the colonization of the Ralstonia genus and provides a new theoretical basis for Si application in PBW prevention.

Single-cell RNA-seq describes the transcriptome landscape and identifies critical transcription factors in the leaf blade of the allotetraploid peanut (Arachis hypogaea L.).

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.

Comparative Transcriptome Analysis Identified Candidate Genes for Late Leaf Spot Resistance and Cause of Defoliation in Groundnut.

Late leaf spot (LLS) caused by fungus Nothopassalora personata in groundnut is responsible for up to 50% yield loss. To dissect the complex nature of LLS resistance, comparative transcriptome analysis was performed using resistant (GPBD 4), susceptible (TAG 24) and a resistant introgression line (ICGV 13208) and identified a total of 12,164 and 9954 DEGs (differentially expressed genes) respectively in A- and B-subgenomes of tetraploid groundnut. There were 135 and 136 unique pathways triggered in A- and B-subgenomes, respectively, upon N. personata infection. Highly upregulated putative disease resistance genes, an RPP-13 like (Aradu.P20JR) and a NBS-LRR (Aradu.Z87JB) were identified on chromosome A02 and A03, respectively, for LLS resistance. Mildew resistance Locus (MLOs)-like proteins, heavy metal transport proteins, and ubiquitin protein ligase showed trend of upregulation in susceptible genotypes, while tetratricopeptide repeats (TPR), pentatricopeptide repeat (PPR), chitinases, glutathione S-transferases, purple acid phosphatases showed upregulation in resistant genotypes. However, the highly expressed ethylene responsive factor (ERF) and ethylene responsive nuclear protein (ERF2), and early responsive dehydration gene (ERD) might be related to the possible causes of defoliation in susceptible genotypes. The identified disease resistance genes can be deployed in genomics-assisted breeding for development of LLS resistant cultivars to reduce the yield loss in groundnut.

Improving Gene Annotation of the Peanut Genome by Integrated Proteogenomics Workflow.

Peanut (Arachis hypogaea L.) is a staple crop in semiarid tropical and subtropical regions. Although the genome of peanut has been fully sequenced, the current gene annotations are still incomplete. New technologies in genomics and proteomics have resulted in the emergence of proteogenomics, which can integrate genomic, transcriptomic, and proteomic data for improving gene annotation. In the present study, we collected RNA-seq and proteomic data from multiple tissues such as seed, shell, and gynophore of peanut and utilized a proteogenomic approach to improve the gene annotation of peanut based on these data. A total of 1 935 655 904 RNA-seq reads and 7 490 280 MS/MS spectra were collected. Ultimately, 13 767 annotated genes were found with evidence at the protein level, and seven novel protein-coding genes were found with both RNA-seq and proteomics evidence. In addition, 35 gene models were updated based on proteomics data. Proteogenomic approaches improved the gene annotation in certain aspects by integrating both RNA-seq and proteomic data. We expect that these approaches could help improve existing genome annotations of other species.

Genome-wide identification of microsatellite markers from cultivated peanut (Arachis hypogaea L.).

BACKGROUND: Microsatellites, or simple sequence repeats (SSRs), represent important DNA variations that are widely distributed across the entire plant genome and can be used to develop SSR markers, which can then be used to conduct genetic analyses and molecular breeding. Cultivated peanut (A. hypogaea L.), an important oil crop worldwide, is an allotetraploid (AABB, 2n = 4× = 40) plant species. Because of its complex genome, genomic marker development has been very challenging. However, sequencing of cultivated peanut genome allowed us to develop genomic markers and construct a high-density physical map. RESULTS: A total of 8,329,496 SSRs were identified, including 3,772,653, 4,414,961, and 141,882 SSRs that were distributed in subgenome A, B, and nine scaffolds, respectively. Based on the flanking sequences of the identified SSRs, a total of 973,984 newly developed SSR markers were developed in subgenome A (462,267), B (489,394), and nine scaffolds (22,323), with an average density of 392.45 markers per Mb. In silico PCR evaluation showed that an average of 88.32% of the SSR markers generated only one in silico-specific product in two tetraploid A. hypogaea varieties, Tifrunner and Shitouqi. A total of 39,599 common SSR markers were identified among the two A. hypogaea varieties and two progenitors, A. duranensis and A. ipaensis. Additionally, an amplification effectiveness of 44.15% was observed by real PCR validation. Moreover, a total of 1276 public SSR loci were integrated with the newly developed SSR markers. Finally, a previously known leaf spot quantitative trait locus (QTL), qLLS_T13_A05_7, was determined to be in a 1.448-Mb region on chromosome A05. In this region, a total of 819 newly developed SSR markers were located and 108 candidate genes were detected. CONCLUSIONS: The availability of these newly developed and public SSR markers both provide a large number of molecular markers that could potentially be used to enhance the process of trait genetic analyses and improve molecular breeding strategies for cultivated peanut.

Draft genome of the peanut A-genome progenitor (Arachis duranensis) provides insights into geocarpy, oil biosynthesis, and allergens.

Peanut or groundnut (Arachis hypogaea L.), a legume of South American origin, has high seed oil content (45-56%) and is a staple crop in semiarid tropical and subtropical regions, partially because of drought tolerance conferred by its geocarpic reproductive strategy. We present a draft genome of the peanut A-genome progenitor, Arachis duranensis, and 50,324 protein-coding gene models. Patterns of gene duplication suggest the peanut lineage has been affected by at least three polyploidizations since the origin of eudicots. Resequencing of synthetic Arachis tetraploids reveals extensive gene conversion in only three seed-to-seed generations since their formation by human hands, indicating that this process begins virtually immediately following polyploid formation. Expansion of some specific gene families suggests roles in the unusual subterranean fructification of Arachis For example, the S1Fa-like transcription factor family has 126 Arachis members, in contrast to no more than five members in other examined plant species, and is more highly expressed in roots and etiolated seedlings than green leaves. The A. duranensis genome provides a major source of candidate genes for fructification, oil biosynthesis, and allergens, expanding knowledge of understudied areas of plant biology and human health impacts of plants, informing peanut genetic improvement and aiding deeper sequencing of Arachis diversity.

Transcriptomic Analysis Reveals the High-Oleic Acid Feedback Regulating the Homologous Gene Expression of Stearoyl-ACP Desaturase 2 (SAD2) in Peanuts.

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 H2O2 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.

Consensus map integration and QTL meta-analysis narrowed a locus for yield traits to 0.7 cM and refined a region for late leaf spot resistance traits to 0.38 cM on linkage group A05 in peanut (Arachis hypogaea L.).

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

TALEN-mediated targeted mutagenesis of fatty acid desaturase 2 (FAD2) in peanut (Arachis hypogaea L.) promotes the accumulation of oleic acid.

KEY MESSAGE: A first creation of high oleic acid peanut varieties by using transcription activator-like effecter nucleases (TALENs) mediated targeted mutagenesis of Fatty Acid Desaturase 2 (FAD2). Transcription activator like effector nucleases (TALENs), which allow the precise editing of DNA, have already been developed and applied for genome engineering in diverse organisms. However, they are scarcely used in higher plant study and crop improvement, especially in allopolyploid plants. In the present study, we aimed to create targeted mutagenesis by TALENs in peanut. Targeted mutations in the conserved coding sequence of Arachis hypogaea fatty acid desaturase 2 (AhFAD2) were created by TALENs. Genetic stability of AhFAD2 mutations was identified by DNA sequencing in up to 9.52 and 4.11% of the regeneration plants at two different targeted sites, respectively. Mutation frequencies among AhFAD2 mutant lines were significantly correlated to oleic acid accumulation. Genetically, stable individuals of positive mutant lines displayed a 0.5-2 fold increase in the oleic acid content compared with non-transgenic controls. This finding suggested that TALEN-mediated targeted mutagenesis could increase the oleic acid content in edible peanut oil. Furthermore, this was the first report on peanut genome editing event, and the obtained high oleic mutants could serve for peanut breeding project.

Identification of the Candidate Proteins Related to Oleic Acid Accumulation during Peanut (Arachis hypogaea L.) Seed Development through Comparative Proteome Analysis.

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.

Transcriptome-wide sequencing provides insights into geocarpy in peanut (Arachis hypogaea L.).

A characteristic feature of peanut is the subterranean fructification, geocarpy, in which the gynophore ('peg'), a specialized organ that transitions from upward growth habit to downward outgrowth upon fertilization, drives the developing pod into the soil for subsequent development underground. As a step towards understanding this phenomenon, we explore the developmental dynamics of the peanut pod transcriptome at 11 successive stages. We identified 110 217 transcripts across developmental stages and quantified their abundance along a pod developmental gradient in pod wall. We found that the majority of transcripts were differentially expressed along the developmental gradient as well as identified temporal programs of gene expression, including hundreds of transcription factors. Thought to be an adaptation to particularly harsh subterranean environments, both up- and down-regulated gene sets in pod wall were enriched for response to a broad array of stimuli, like gravity, light and subterranean environmental factors. We also identified hundreds of transcripts associated with gravitropism and photomorphogenesis, which may be involved in the geocarpy. Collectively, this study forms a transcriptional baseline for geocarpy in peanut as well as provides a considerable body of evidence that transcriptional regulation in peanut aerial and subterranean fruits is complex.

Comparative transcriptome analysis of aerial and subterranean pods development provides insights into seed abortion in peanut.

The peanut is a special plant for its aerial flowering but subterranean fructification. The failure of peg penetration into the soil leads to form aerial pod and finally seed abortion. However, the mechanism of seed abortion during aerial pod development remains obscure. Here, a comparative transcriptome analysis between aerial and subterranean pods at different developmental stages was produced using a customized NimbleGen microarray representing 36,158 unigenes. By comparing 4 consecutive time-points, totally 6,203 differentially expressed genes, 4,732 stage-specific expressed genes and 2,401 specific expressed genes only in aerial or subterranean pods were identified in this study. Functional annotation showed their mainly involvement in biosynthesis, metabolism, transcription regulation, transporting, stress response, photosynthesis, signal transduction, cell division, apoptosis, embryonic development, hormone response and light signaling, etc. Emphasis was focused on hormone response, cell apoptosis, embryonic development and light signaling relative genes. These genes might function as potential candidates to provide insights into seed abortion during aerial pod development. Ten candidate genes were validated by Real-time RT-PCR. Additionally, consistent with up-regulation of auxin response relative genes in aerial pods, endogenous IAA content was also significantly increased by HPLC analysis. This study will further provide new molecular insight that auxin and auxin response genes potentially contribute to peanut seed and pod development.

A Single-Nucleus Resolution Atlas of Transcriptome and Chromatin Accessibility for Peanut (Arachis Hypogaea L.) Leaves.

The peanut is an important worldwide cash-crop for edible oil and protein. However, the kinetic mechanisms that determine gene expression and chromatin accessibility during leaf development in peanut represented allotetraploid leguminous crops are poorly understood at single-cell resolution. Here, a single-nucleus atlas of peanut leaves is developed by simultaneously profiling the transcriptome and chromatin accessibility in the same individual-cell using fluorescence-activated sorted single-nuclei. In total, 5930 cells with 50 890 expressed genes are classified into 18 cell-clusters, and 5315 chromatin fragments are enriched with 26 083 target genes in the chromatin accessible landscape. The developmental trajectory analysis reveals the involvement of the ethylene-AP2 module in leaf cell differentiation, and cell-cycle analysis demonstrated that genome replication featured in distinct cell-types with circadian rhythms transcription factors (TFs). Furthermore, dual-omics illustrates that the fatty acid pathway modulates epidermal-guard cells differentiation and providescritical TFs interaction networks for understanding mesophyll development, and the cytokinin module (LHY/LOG) that regulates vascular growth. Additionally, an AT-hook protein AhAHL11 is identified that promotes leaf area expansion by modulating the auxin content increase. In summary, the simultaneous profiling of transcription and chromatin accessibility landscapes using snRNA/ATAC-seq provides novel biological insights into the dynamic processes of peanut leaf cell development at the cellular level.

Silica nanoparticles conferring resistance to bacterial wilt in peanut (Arachis hypogaea L.).

Peanut bacterial wilt (PBW) caused by the pathogen Ralstonia solanacearum severely affects the growth and yield potential of peanut crop. In this study, we synthesized silica nanoparticles (SiO2 NPs), a prospective efficient management approach to control PBW, and conducted a hydroponic experiment to investigate the effects of different SiO2 NPs treatments (i.e., 0, 100, and 500 mg L-1 as NP0, NP100, and NP500, respectively) on promoting plant growth and resistance to R. solanacearum. Results indicated that the disease indices of NP100 and NP500 decreased by 51.5 % and 55.4 % as compared with NP0 under R. solanacearum inoculation, respectively, while the fresh and dry weights and shoot length of NP100 and NP500 increased by 7.62-42.05 %, 9.45-32.06 %, and 2.37-17.83 %, respectively. Furthermore, SiO2 NPs induced an improvement in physio-biochemical enzymes (superoxide dismutase, peroxidase, catalase, ascorbate peroxidase, and lipoxygenase) which eliminated the excess production of hydrogen peroxide, superoxide anions, and malondialdehyde to alleviate PBW stress. Notably, the targeted metabolomic analysis indicated that SiO2 NPs enhanced salicylic acid (SA) contents, which involved the induction of systemic acquired resistance (SAR). Moreover, the transcriptomic analysis revealed that SiO2 NPs modulated the expression of multiple transcription factors (TFs) involved in the hormone pathway, such as AHLs, and the identification of hormone pathways related to plant defense responses, such as the SA pathway, which activated SA-dependent defense mechanisms. Meanwhile, the up-regulated expression of the SA-metabolism gene, salicylate carboxymethyltransferase (SAMT), initiated SAR to promote PBW resistance. Overall, our findings revealed that SiO2 NPs, functioning as a plant elicitor, could effectively modulate physiological enzyme activities and enhance SA contents through the regulation of SA-metabolism genes to confer the PBW resistance in peanuts, which highlighted the potential of SiO2 NPs for sustainable agricultural practices.

A genomic variation map provides insights into peanut diversity in China and associations with 28 agronomic traits.

Peanut (Arachis hypogaea L.) is an important allotetraploid oil and food legume crop. China is one of the world's largest peanut producers and consumers. However, genomic variations underlying the migration and divergence of peanuts in China remain unclear. Here we reported a genome-wide variation map based on the resequencing of 390 peanut accessions, suggesting that peanuts might have been introduced into southern and northern China separately, forming two cultivation centers. Selective sweep analysis highlights asymmetric selection between the two subgenomes during peanut improvement. A classical pedigree from South China offers a context for the examination of the impact of artificial selection on peanut genome. Genome-wide association studies identified 22,309 significant associations with 28 agronomic traits, including candidate genes for plant architecture and oil biosynthesis. Our findings shed light on peanut migration and diversity in China and provide valuable genomic resources for peanut improvement.

Deep sequencing analysis of the transcriptomes of peanut aerial and subterranean young pods identifies candidate genes related to early embryo abortion.

The failure of peg penetration into the soil leads to seed abortion in peanut. Knowledge of genes involved in these processes is comparatively deficient. Here, we used RNA-seq to gain insights into transcriptomes of aerial and subterranean pods. More than 2 million transcript reads with an average length of 396 bp were generated from one aerial (AP) and two subterranean (SP1 and SP2) pod libraries using pyrosequencing technology. After assembly, sets of 49 632, 49 952 and 50 494 from a total of 74 974 transcript assembly contigs (TACs) were identified in AP, SP1 and SP2, respectively. A clear linear relationship in the gene expression level was observed between these data sets. In brief, 2194 differentially expressed TACs with a 99.0% true-positive rate were identified, among which 859 and 1068 TACs were up-regulated in aerial and subterranean pods, respectively. Functional analysis showed that putative function based on similarity with proteins catalogued in UniProt and gene ontology term classification could be determined for 59 342 (79.2%) and 42 955 (57.3%) TACs, respectively. A total of 2968 TACs were mapped to 174 KEGG pathways, of which 168 were shared by aerial and subterranean transcriptomes. TACs involved in photosynthesis were significantly up-regulated and enriched in the aerial pod. In addition, two senescence-associated genes were identified as significantly up-regulated in the aerial pod, which potentially contribute to embryo abortion in aerial pods, and in turn, to cessation of swelling. The data set generated in this study provides evidence for some functional genes as robust candidates underlying aerial and subterranean pod development and contributes to an elucidation of the evolutionary implications resulting from fruit development under light and dark conditions.

Overexpression of peanut (Arachis hypogaea L.) AhGRFi gene enhanced root growth inhibition under exogenous NAA treatment in Arabidopsis thaliana.

The 14-3-3 protein is a kind of evolutionary ubiquitous protein family highly conserved in eukaryotes. Initially, 14-3-3 proteins were reported in mammalian nervous tissues, but in the last decade, their role in various metabolic pathways in plants established the importance of 14-3-3 proteins. In the present study, a total of 22 14-3-3 genes, also called general regulatory factors (GRF), were identified in the peanut (Arachis hypogaea) genome, out of which 12 belonged to the ε group, whereas 10 of them belonged to the non- ε-group. Tissue-specific expression of identified 14-3-3 genes were studied using transcriptome analysis. The peanut AhGRFi gene was cloned and transformed into Arabidopsis thaliana. The investigation of subcellular localization indicated that AhGRFi is localized in the cytoplasm. Overexpression of the AhGRFi gene in transgenic Arabidopsis showed that under exogenous 1-naphthaleneacetic acid (NAA) treatment, root growth inhibition in transgenic plants was enhanced. Further analysis indicated that the expression of auxin-responsive genes IAA3, IAA7, IAA17, and SAUR-AC1 was upregulated and GH3.2 and GH3.3 were downregulated in transgenic plants, but the expression of GH3.2, GH3.3, and SAUR-AC1 showed opposite trends of change under NAA treatment. These results suggest that AhGRFi may be involved in auxin signaling during seedling root development. An in-depth study of the molecular mechanism of this process remains to be further explored.

scRNA-seq Reveals the Mechanism of Fatty Acid Desaturase 2 Mutation to Repress Leaf Growth in Peanut (Arachis hypogaea L.).

Fatty Acid Desaturase 2 (FAD2) controls the conversion of oleic acids into linoleic acids. Mutations in FAD2 not only increase the high-oleic content, but also repress the leaf growth. However, the mechanism by which FAD2 regulates the growth pathway has not been elucidated in peanut leaves with single-cell resolution. In this study, we isolated fad2 mutant leaf protoplast cells to perform single-cell RNA sequencing. Approximately 24,988 individual cells with 10,249 expressed genes were classified into five major cell types. A comparative analysis of 3495 differentially expressed genes (DEGs) in distinct cell types demonstrated that fad2 inhibited the expression of the cytokinin synthesis gene LOG in vascular cells, thereby repressing leaf growth. Further, pseudo-time trajectory analysis indicated that fad2 repressed leaf cell differentiation, and cell-cycle evidence displayed that fad2 perturbed the normal cell cycle to induce the majority of cells to drop into the S phase. Additionally, important transcription factors were filtered from the DEG profiles that connected the network involved in high-oleic acid accumulation (WRKY6), activated the hormone pathway (WRKY23, ERF109), and potentially regulated leaf growth (ERF6, MYB102, WRKY30). Collectively, our study describes different gene atlases in high-oleic and normal peanut seedling leaves, providing novel biological insights to elucidate the molecular mechanism of the high-oleic peanut-associated agronomic trait at the single-cell level.