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

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.

Single-nucleus RNA and ATAC sequencing analyses provide molecular insights into early pod development of peanut fruit

Peanut (Arachis hypogaea L.) is an important leguminous oil and economic crop that produces flowers aboveground and fruits belowground. Subterranean fruit-pod development, which significantly affects peanut production, involves complex molecular mechanisms that likely require the coordinated regulation of multiple genes in different tissues. To investigate the molecular mechanisms that underlie peanut fruit-pod development, we characterized the anatomical features of early fruit-pod development and integrated single-nucleus RNA-sequencing (snRNA-seq) and single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) data at the single-cell level. We identified distinct cell types, such as meristem, embryo, vascular tissue, cuticular layer, and stele cells within the shell wall. These specific cell types were used to examine potential molecular changes unique to each cell type during pivotal stages of fruit-pod development. snRNA-seq analyses of differentially expressed genes revealed cell-type-specific insights that were not previously obtainable from transcriptome analyses of bulk RNA. For instance, we identified MADS-box genes that contributes to the formation of parenchyma cells and gravity-related genes that are present in the vascular cells, indicating an essential role for the vascular cells in peg gravitropism. Overall, our single-nucleus analysis provides comprehensive and novel information on specific cell types, gene expression, and chromatin accessibility during the early stages of fruit-pod development. This information will enhance our understanding of the mechanisms that underlie fruit-pod development in peanut and contribute to efforts aimed at improving peanut production.

Omics-driven advances in the understanding of regulatory landscape of peanut seed development

Peanuts (Arachis hypogaea) are an essential oilseed crop known for their unique developmental process, characterized by aerial flowering followed by subterranean fruit development. This crop is polyploid, consisting of A and B subgenomes, which complicates its genetic analysis. The advent and progression of omics technologies-encompassing genomics, transcriptomics, proteomics, epigenomics, and metabolomics-have significantly advanced our understanding of peanut biology, particularly in the context of seed development and the regulation of seed-associated traits. Following the completion of the peanut reference genome, research has utilized omics data to elucidate the quantitative trait loci (QTL) associated with seed weight, oil content, protein content, fatty acid composition, sucrose content, and seed coat color as well as the regulatory mechanisms governing seed development. This review aims to summarize the advancements in peanut seed development regulation and trait analysis based on reference genome-guided omics studies. It provides an overview of the significant progress made in understanding the molecular basis of peanut seed development, offering insights into the complex genetic and epigenetic mechanisms that influence key agronomic traits. These studies highlight the significance of omics data in profoundly elucidating the regulatory mechanisms of peanut seed development. Furthermore, they lay a foundational basis for future research on trait-related functional genes, highlighting the pivotal role of comprehensive genomic analysis in advancing our understanding of plant biology.

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.

Depletion of protective microbiota promotes the incidence of fruit disease

Plant-associated microbiomes play important roles in plant health and productivity. However, despite fruits being directly linked to plant productivity, little is known about the microbiomes of fruits and their potential association with fruit health. Here, by integrating 16S rRNA gene, ITS high-throughput sequencing data, and microbiological culturable approaches, we reported that roots and fruits (pods) of peanut, a typical plant that bears fruits underground, recruit different bacterial and fungal communities independently of cropping conditions and that the incidence of pod disease under monocropping conditions is attributed to the depletion of Bacillus genus and enrichment of Aspergillus genus in geocarposphere. On this basis, we constructed a synthetic community (SynCom) consisting of three Bacillus strains from geocarposphere soil under rotation conditions with high culturable abundance. Comparative transcriptome, microbiome profiling, and plant phytohormone signaling analysis reveal that the SynCom exhibited more effective Aspergillus growth inhibition and pod disease control than individual strain, which was underpinned by a combination of molecular mechanisms related to fungal cell proliferation interference, mycotoxins biosynthesis impairment, and jasmonic acid-mediated plant immunity activation. Overall, our results reveal the filter effect of plant organs on the microbiome and that depletion of key protective microbial community promotes the fruit disease incidence.

Genome-wide analysis of PIN genes in cultivated peanuts (Arachis hypogaea L.): identification, subcellular localization, evolution, and expression patterns

BACKGROUND

Auxin is an important hormone in plants and the PIN-FORMED (PIN) genes are essential to auxin distribution in growth and developmental processes of plants. Peanut is an influential cash crop, but research into PIN genes in peanuts remains limited.

RESULTS

In this study, 16 PIN genes were identified in the genome of cultivated peanut, resolving into four subfamilies. All PIN genes were predicted to be located in the plasma membrane and a subcellular location experiment confirmed this prediction for eight of them. The gene structure, cis-elements in the promoter, and evolutionary relationships were elucidated, facilitating our understanding of peanut PINs and their evolution. In addition, the expression patterns of these PINs in various tissues were analyzed according to a previously published transcriptome dataset and qRT-PCR, which gave us a clear understanding of the temporal and spatial expression of PIN genes in different growth stages and different tissues. The expression trend of homologous genes was similar. AhPIN2A and AhPIN2B exhibited predominant expression in roots. AhPIN1A-1 and AhPIN1B-1 displayed significant upregulation following peg penetration, suggesting a potential close association with peanut pod development. Furthermore, we presented the gene network and gene ontology enrichment of these PINs. Notably, AhABCB19 exhibited a co-expression relationship with AhPIN1A and AhPIN1B-1, with all three genes displaying higher expression levels in peanut pegs and pods. These findings reinforce their potential role in peanut pod development.

CONCLUSIONS

This study details a comprehensive analysis of PIN genes in cultivated peanuts and lays the foundation for subsequent studies of peanut gene function and phenotype.

Diversity and Pathobiology of an Ilarvirus Unexpectedly Detected in Diverse Plants and Global Sequencing Data

High-throughput sequencing (HTS) and sequence mining tools revolutionized virus detection and discovery in recent years, and implementing them with classical plant virology techniques results in a powerful approach to characterize viruses. An example of a virus discovered through HTS is Solanum nigrum ilarvirus 1 (SnIV1) (Bromoviridae), which was recently reported in various solanaceous plants from France, Slovenia, Greece, and South Africa. It was likewise detected in grapevines (Vitaceae) and several Fabaceae and Rosaceae plant species. Such a diverse set of source organisms is atypical for ilarviruses, thus warranting further investigation. In this study, modern and classical virological tools were combined to accelerate the characterization of SnIV1. Through HTS-based virome surveys, mining of sequence read archive datasets, and a literature search, SnIV1 was further identified from diverse plant and non-plant sources globally. SnIV1 isolates showed relatively low variability compared with other phylogenetically related ilarviruses. Phylogenetic analyses showed a distinct basal clade of isolates from Europe, whereas the rest formed clades of mixed geographic origin. Furthermore, systemic infection of SnIV1 in Solanum villosum and its mechanical and graft transmissibility to solanaceous species were demonstrated. Near-identical SnIV1 genomes from the inoculum (S. villosum) and inoculated Nicotiana benthamiana were sequenced, thus partially fulfilling Koch's postulates. SnIV1 was shown to be seed-transmitted and potentially pollen-borne, has spherical virions, and possibly induces histopathological changes in infected N. benthamiana leaf tissues. Overall, this study provides information to better understand the diversity, global presence, and pathobiology of SnIV1; however, its possible emergence as a destructive pathogen remains uncertain. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Fine mapping of qAHPS07 and functional studies of AhRUVBL2 controlling pod size in peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.) is an important oil and cash crop. Pod size is one of the major traits determining yield and commodity characteristic of peanut. Fine mapping of quantitative trait locus (QTL) and identification of candidate genes associated with pod size are essential for genetic improvement and molecular breeding of peanut varieties. In this study, a major QTL related to pod size, qAHPS07, was fine mapped to a 36.46 kb interval on chromosome A07 using F2 , recombinant inbred line (RIL) and secondary F2 populations. qAHPS07 explained 38.6%, 23.35%, 37.48%, 25.94% of the phenotypic variation for single pod weight (SPW), pod length (PL), pod width (PW) and pod shell thickness (PST), respectively. Whole genome resequencing and gene expression analysis revealed that a RuvB-like 2 protein coding gene AhRUVBL2 was the most likely candidate for qAHPS07. Overexpression of AhRUVBL2 in Arabidopsis led to larger seeds and plants than the wild type. AhRUVBL2-silenced peanut seedlings represented small leaves and shorter main stems. Three haplotypes were identified according to three SNPs in the promoter of AhRUVBL2 among 119 peanut accessions. Among them, SPW, PW and PST of accessions carrying Hap_ATT represent 17.6%, 11.2% and 26.3% higher than those carrying Hap_GAC，respectively. In addition, a functional marker of AhRUVBL2 was developed. Taken together, our study identified a key functional gene of peanut pod size, which provides new insights into peanut pod size regulation mechanism and offers practicable markers for the genetic improvement of pod size-related traits in peanut breeding.

Global Transcriptome and Co-Expression Network Analyses Revealed Hub Genes Controlling Seed Size/Weight and/or Oil Content in Peanut

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin-proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut.

Multi-Omics Profiling Identifies Candidate Genes Controlling Seed Size in Peanut

Seed size is the major yield component and a key target trait that is selected during peanut breeding. However, the mechanisms that regulate peanut seed size are unknown. Two peanut mutants with bigger seed size were isolated in this study by ⁶⁰Co treatment of a common peanut landrace, Huayu 22, and were designated as the "big seed" mutant lines (hybs). The length and weight of the seed in hybs were about 118% and 170% of those in wild-type (WT), respectively. We adopted a multi-omics approach to identify the genomic locus underlying the hybs mutants. We performed whole genome sequencing (WGS) of WT and hybs mutants and identified thousands of large-effect variants (SNPs and indels) that occurred in about four hundred genes in hybs mutants. Seeds from both WT and hybs lines were sampled 20 days after flowering (DAF) and were used for RNA-Seq analysis; the results revealed about one thousand highly differentially expressed genes (DEGs) in hybs compared to WT. Using a method that combined large-effect variants with DEGs, we identified 45 potential candidate genes that shared gene product mutations and expression level changes in hybs compared to WT. Among the genes, two candidate genes encoding cytochrome P450 superfamily protein and NAC transcription factors may be associated with the increased seed size in hybs. The present findings provide new information on the identification and functional research into candidate genes responsible for the seed size phenotype in peanut.

Analysis of the Transcriptional Dynamics of Regulatory Genes During Peanut Pod Development Caused by Darkness and Mechanical Stress

Peanut is an oil crop with important economic value that is widely cultivated around the world. It blooms on the ground but bears fruit underground. When the peg penetrates the ground, it enters a dark environment, is subjected to mechanical stress from the soil, and develops into a normal pod. When a newly developed pod emerges from the soil, it turns green and stops growing. It has been reported that both darkness and mechanical stress are necessary for normal pod development. In this study, we investigated changes in gene expression during the reverse process of peg penetration: developmental arrest caused by pod (Pattee 3 pods) excavation. Bagging the aerial pods was used to simulate loss of mechanical pressure, while direct exposure of the aerial pods was used to simulate loss of both mechanical pressure and darkness. After the loss of mechanical stress and darkness, the DEGs were significantly enriched in photosynthesis, photosynthesis-antenna proteins, plant-pathogen interaction, DNA replication, and circadian rhythm pathways. The DNA replication pathway was enriched by down-regulated genes, and the other four pathways were enriched by upregulated genes. Upregulated genes were also significantly enriched in protein ubiquitination and calmodulin-related genes, highlighting the important role of ubiquitination and calcium signaling in pod development. Further analysis of DEGs showed that phytochrome A (Phy A), auxin response factor 9 (IAA9), and mechanosensitive ion channel protein played important roles in geocarpy. The expression of these two genes increased in subterranean pods but decreased in aerial pods. Based on a large number of chloroplast-related genes, calmodulin, kinases, and ubiquitin-related proteins identified in this study, we propose two possible signal transduction pathways involved in peanut geocarpy, namely, one begins in chloroplasts and signals down through phosphorylation, and the other begins during abiotic stress and signals down through calcium signaling, phosphorylation, and ubiquitination. Our study provides valuable information about putative regulatory genes for peanut pod development and contributes to a better understanding of the biological phenomenon of geocarpy.

Genome-Wide Identification and Expression of FAR1 Gene Family Provide Insight Into Pod Development in Peanut (Arachis hypogaea)

The far-red-impaired response 1 (FAR1) transcription family were initially identified as important factors for phytochrome A (phyA)-mediated far-red light signaling in Arabidopsis; they play crucial roles in controlling the growth and development of plants. The reported reference genome sequences of Arachis, including A. duranensis, A. ipaensis, A. monticola, and A. hypogaea, and its related species Glycine max provide an opportunity to systematically perform a genome-wide identification of FAR1 homologous genes and investigate expression patterns of these members in peanut species. Here, a total of 650 FAR1 genes were identified from four Aarchis and its closely related species G. max. Of the studied species, A. hypogaea contained the most (246) AhFAR1 genes, which can be classified into three subgroups based on phylogenic relationships. The synonymous (Ks) and non-synonymous (Ka) substitution rates, phylogenetic relationship and synteny analysis of the FAR1 family provided deep insight into polyploidization, evolution and domestication of peanut AhFAR1 genes. The transcriptome data showed that the AhFAR1 genes exhibited distinct tissue- and stage-specific expression patterns in peanut. Three candidate genes including Ahy_A10g049543, Ahy_A06g026579, and Ahy_A10g048401, specifically expressed in peg and pod, might participate in pod development in the peanut. The quantitative real-time PCR (qRT-PCR) analyses confirmed that the three selected genes were highly and specifically expressed in the peg and pod. This study systematically analyzed gene structure, evolutionary characteristics and expression patterns of FAR1 gene family, which will provide a foundation for the study of genetic and biological function in the future.

Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut

Combining QTL-seq, QTL-mapping and RNA-seq identified a major QTL and candidate genes, which contributed to the development of KASP markers and understanding of molecular mechanisms associated with seed weight in peanut. Seed weight, as an important component of seed yield, is a significant target of peanut breeding. However, relatively little is known about the quantitative trait loci (QTLs) and candidate genes associated with seed weight in peanut. In this study, three major QTLs on chromosomes A05, B02, and B06 were determined by applying the QTL-seq approach in a recombinant inbred line (RIL) population. Based on conventional QTL-mapping, these three QTL regions were successfully narrowed down through newly developed single nucleotide polymorphism (SNP) and simple sequence repeat markers. Among these three QTL regions, qSWB06.3 exhibited stable expression, contributing mainly to phenotypic variance across environments. Furthermore, differentially expressed genes (DEGs) were identified at the three seed developmental stages between the two parents of the RIL population. It was found that the DEGs were widely distributed in the ubiquitin-proteasome pathway, the serine/threonine-protein pathway, signal transduction of hormones and transcription factors. Notably, DEGs at the early stage were mostly involved in regulating cell division, whereas DEGs at the middle and late stages were primarily involved in cell expansion during seed development. The expression patterns of candidate genes related to seed weight in qSWB06.3 were investigated using quantitative real-time PCR. In addition, the allelic diversity of qSWB06.3 was investigated in peanut germplasm accessions. The marker Ah011475 has higher efficiency for discriminating accessions with different seed weights, and it would be useful as a diagnostic marker in marker-assisted breeding. This study provided insights into the genetic and molecular mechanisms of seed weight in peanut.

Deep transcriptomic study reveals the role of cell wall biosynthesis and organization networks in the developing shell of peanut pod

BACKGROUND

Peanut (Arachis hypogaea L.) belongs to an exceptional group of legume plants, wherein the flowers are produced aerially, but the pods develop under the ground. In such a unique environment, the pod's outer shell plays a vital role as a barrier against mechanical damage and soilborne pathogens. Recent studies have reported the uniqueness and importance of gene expression patterns that accompany peanut pods' biogenesis. These studies focused on biogenesis and pod development during the early stages, but the late developmental stages and disease resistance aspects still have gaps. To extend this information, we analyzed the transcriptome generated from four pod developmental stages of two genotypes, Hanoch (Virginia-type) and IGC53 (Peruvian-type), which differs significantly in their pod shell characteristics and pathogen resistance.

RESULTS

The transcriptome study revealed a significant reprogramming of the number and nature of differentially expressed (DE) genes during shell development. Generally, the numbers of DE genes were higher in IGC53 than in Hanoch, and the R5-R6 transition was the most dynamic in terms of transcriptomic changes. Genes related to cell wall biosynthesis, modification and transcription factors (TFs) dominated these changes therefore, we focused on their differential, temporal and spatial expression patterns. Analysis of the cellulose synthase superfamily identified specific Cellulose synthase (CesAs) and Cellulose synthase-like (Csl) genes and their coordinated interplay with other cell wall-related genes during the peanut shell development was demonstrated. TFs were also identified as being involved in the shell development process, and their pattern of expression differed in the two peanut genotypes. The shell component analysis showed that overall crude fiber, cellulose, lignin, hemicelluloses and dry matter increased with shell development, whereas K, N, protein, and ash content decreased. Genotype IGC53 contained a higher level of crude fiber, cellulose, NDF, ADF, K, ash, and dry matter percentage, while Hanoch had higher protein and nitrogen content.

CONCLUSIONS

The comparative transcriptome analysis identified differentially expressed genes, enriched processes, and molecular processes like cell wall biosynthesis/modifications, carbohydrate metabolic process, signaling, transcription factors, transport, stress, and lignin biosynthesis during the peanut shell development between two contrasting genotypes. TFs and other genes like chitinases were also enriched in peanut shells known for pathogen resistance against soilborne major pathogens causing pod wart disease and pod damages. This study will shed new light on the biological processes involved with underground pod development in an important legume crop.

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.

Genome-Wide Identification and Characterization of HSP90-RAR1-SGT1-Complex Members From Arachis Genomes and Their Responses to Biotic and Abiotic Stresses

The molecular chaperone complex HSP90-RAR1-SGT1 (HRS) plays important roles in both biotic and abiotic stress responses in plants. A previous study showed that wild peanut Arachis diogoi SGT1 (AdSGT1) could enhance disease resistance in transgenic tobacco and peanut. However, no systematic analysis of the HRS complex in Arachis has been conducted to date. In this study, a comprehensive analysis of the HRS complex were performed in Arachis. Nineteen HSP90, two RAR1 and six SGT1 genes were identified from the allotetraploid peanut Arachis hypogaea, a number close to the sum of those from the two wild diploid peanut species Arachis duranensis and Arachis ipaensis. According to phylogenetic and chromosomal location analyses, thirteen orthologous gene pairs from Arachis were identified, all of which except AhHSP90-A8, AhHSP90-B9, AdHSP90-9, and AiHSP90-9 were localized on the syntenic locus, and they shared similar exon-intron structures, conserved motifs and expression patterns. Phylogenetic analysis showed that HSP90 and RAR1 from dicot and monocot plants diverged into different clusters throughout their evolution. Chromosomal location analysis indicated that AdSGT1 (the orthologous gene of AhSGT1-B3 in this study) might provide resistance to leaf late spot disease dependent on the orthologous genes of AhHSP90-B10 and AhRAR1-B in the wild peanut A. diogoi. Several HRS genes exhibited tissue-specific expression patterns, which may reflect the sites where they perform functions. By exploring published RNA-seq data, we found that several HSP90 genes play major roles in both biotic and abiotic stress responses, especially salt and drought responses. Autoactivation assays showed that AhSGT1-B1 could not be used as bait for yeast two-hybrid (Y2H) library screening. AhRAR1 and AhSGT1 could strongly interact with each other and interact with AhHSP90-B8. The present study represents the first systematic analysis of HRS complex genes in Arachis and provides valuable information for functional analyses of HRS complex genes. This study also offers potential stress-resistant genes for peanut improvement.

Comprehensive Transcriptome Analyses Reveal Candidate Genes for Variation in Seed Size/Weight During Peanut (Arachis hypogaea L.) Domestication

Seed size/weight, a key domestication trait, is also an important selection target during peanut breeding. However, the mechanisms that regulate peanut seed development are unknown. We re-sequenced 12 RNA samples from developing seeds of two cultivated peanut accessions (Lines 8106 and 8107) and wild Arachis monticola at 15, 30, 45, and 60 days past flowering (DPF). Transcriptome analyses showed that ∼36,000 gene loci were expressed in each of the 12 RNA samples, with nearly half exhibiting moderate (2 ≤ FPKM < 10) expression levels. Of these genes, 12.2% (4,523) were specifically expressed during seed development, mainly at 15 DPF. Also, ∼12,000 genes showed significant differential expression at 30, 45, and/or 60 DPF within each of the three peanut accessions, accounting for 31.8-34.1% of the total expressed genes. Using a method that combined comprehensive transcriptome analysis and previously mapped QTLs, we identified several candidate genes that encode transcription factor TGA7, topless-related protein 2, IAA-amino acid hydrolase ILR1-like 5, and putative pentatricopeptide repeat-containing (PPR) protein. Based on sequence variations identified in these genes, SNP markers were developed and used to genotype both 30 peanut landraces and a genetic segregated population, implying that EVM0025654 encoding a PPR protein may be associated with the increased seed size/weight of the cultivated accessions in comparison with the allotetraploid wild peanut. Our results provide additional knowledge for the identification and functional research into candidate genes responsible for the seed size/weight phenotype in peanut.

A proteomic analysis of peanut seed at different stages of underground development to understand the changes of seed proteins

In order to obtain more valuable insights into the protein dynamics and accumulation of allergens in seeds during underground development, we performed a proteomic study on developing peanut seeds at seven different stages. A total of 264 proteins with altered abundance and contained at least one unique peptide was detected by matrix-assisted laser desorption ionization time-of-flight/time-of-flight mass spectrometry (MALDI-TOF/TOF MS). All identified proteins were classified into five functional categories as level 1 and 20 secondary functional categories as level 2. Among them, 88 identified proteins (IPs) were related to carbohydrate/ amino acid/ lipid transport and metabolism, indicating that carbohydrate/amino acid/ lipid metabolism played a key role in the underground development of peanut seeds. Hierarchical cluster analysis showed that all IPs could be classified into eight cluster groups according to the abundance profiles, suggesting that the modulatory patterns of these identified proteins were complicated during seed development. The largest group contained 41 IPs, the expression of which decreased at R 2 and reached a maximum at R3 but gradually decreased from R4. A total of 14 IPs were identified as allergen-like proteins by BLAST with A genome (Arachis duranensis) or B genome (Arachis ipaensis) translated allergen sequences. Abundance profile analysis of 14 identified allergens showed that the expression of all allergen proteins was low or undetectable by 2-DE at the early stages (R1 to R4), and began to accumulate from the R5 stage and gradually increased. Network analysis showed that most of the significant proteins were involved in active metabolic pathways in early development. Real time RT-PCR analysis revealed that transcriptional regulation was approximately consistent with expression at the protein level for 8 selected identified proteins. In addition, some amino acid sequences that may be associated with new allergens were also discussed.

Arachis hypogaea gene expression atlas for fastigiata subspecies of cultivated groundnut to accelerate functional and translational genomics applications

Spatio-temporal and developmental stage-specific transcriptome analysis plays a crucial role in systems biology-based improvement of any species. In this context, we report here the Arachis hypogaea gene expression atlas (AhGEA) for the world's widest cultivated subsp. fastigiata based on RNA-seq data using 20 diverse tissues across five key developmental stages. Approximately 480 million paired-end filtered reads were generated followed by identification of 81 901 transcripts from an early-maturing, high-yielding, drought-tolerant groundnut variety, ICGV 91114. Further, 57 344 genome-wide transcripts were identified with ≥1 FPKM across different tissues and stages. Our in-depth analysis of the global transcriptome sheds light into complex regulatory networks namely gravitropism and photomorphogenesis, seed development, allergens and oil biosynthesis in groundnut. Importantly, interesting insights into molecular basis of seed development and nodulation have immense potential for translational genomics research. We have also identified a set of stable expressing transcripts across the selected tissues, which could be utilized as internal controls in groundnut functional genomics studies. The AhGEA revealed potential transcripts associated with allergens, which upon appropriate validation could be deployed in the coming years to develop consumer-friendly groundnut varieties. Taken together, the AhGEA touches upon various important and key features of cultivated groundnut and provides a reference for further functional, comparative and translational genomics research for various economically important traits.

Alternative splicing profiling provides insights into the molecular mechanisms of peanut peg development

BACKGROUND

The cultivated peanut (Arachis hypogaea) is one of the most important oilseed crops worldwide, and the generation of pegs and formation of subterranean pods are essential processes in peanut reproductive development. However, little information has been reported about alternative splicing (AS) in peanut peg formation and development.

RESULTS

Herein, we presented a comprehensive full-length (FL) transcriptome profiling of AS isoforms during peanut peg and early pod development. We identified 1448, 1102, 832, and 902 specific spliced transcripts in aerial pegs, subterranean pegs, subterranean unswollen pegs, and early swelling pods, respectively. A total of 184 spliced transcripts related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins were identified as possibly involved in peanut peg development. For aerial pegs, spliced transcripts we got were mainly involved in gravity stimulation and cell wall morphogenetic processes. The genes undergoing AS in subterranean peg were possibly involved in gravity stimulation, cell wall morphogenetic processes, and abiotic response. For subterranean unswollen pegs, spliced transcripts were predominantly related to the embryo development and root formation. The genes undergoing splice in early swelling pods were mainly related to ovule development, root hair cells enlargement, root apex division, and seed germination.

CONCLUSION

This study provides evidence that multiple genes are related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins undergoing AS express development-specific spliced isoforms or exhibit an obvious isoform switch during the peanut peg development. AS isoforms in subterranean pegs and pods provides valuable sources to further understand post-transcriptional regulatory mechanisms of AS in the generation of pegs and formation of subterranean pods.

Global transcriptome analysis of subterranean pod and seed in peanut (Arachis hypogaea L.) unravels the complexity of fruit development under dark condition

Peanut pods develop underground, which is the most salient characteristic in peanut. However, its developmental transcriptome remains largely unknown. In the present study, we sequenced over one billion transcripts to explore the developmental transcriptome of peanut pod using Illumina sequencing. Moreover, we identified and quantified the abundances of 165,689 transcripts in seed and shell tissues along with a pod developmental gradient. The dynamic changes of differentially expressed transcripts (DETs) were described in seed and shell. Additionally, we found that photosynthetic genes were not only pronouncedly enriched in aerial pod, but also played roles in developing pod under dark condition. Genes functioning in photomorphogenesis showed distinct expression profiles along subterranean pod development. Clustering analysis unraveled a dynamic transcriptome, in which transcripts for DNA synthesis and cell division during pod expansion were transitioning to transcripts for cell expansion and storage activity during seed filling. Collectively, our study formed a transcriptional baseline for peanut fruit development under dark condition.

Defining the function of SUMO system in pod development and abiotic stresses in Peanut

BACKGROUND

Posttranslational modification of proteins by small ubiquitin like modifier (SUMO) proteins play an important role during the developmental process and in response to abiotic stresses in plants. However, little is known about SUMOylation in peanut (Arachis hypogaea L.), one of the world's major food legume crops. In this study, we characterized the SUMOylation system from the diploid progenitor genomes of peanut, Arachis duranensis (AA) and Arachis ipaensis (BB).

RESULTS

Genome-wide analysis revealed the presence of 40 SUMO system genes in A. duranensis and A. ipaensis. Our results showed that peanut also encodes a novel class II isotype of the SCE1, which was previously reported to be uniquely present in cereals. RNA-seq data showed that the core components of the SUMOylation cascade SUMO1/2 and SCE1 genes exhibited pod-specific expression patterns, implying coordinated regulation during pod development. Furthermore, both transcripts and conjugate profiles revealed that SUMOylation has significant roles during the pod development. Moreover, dynamic changes in the SUMO conjugates were observed in response to abiotic stresses.

CONCLUSIONS

The identification and organization of peanut SUMO system revealed SUMOylation has important roles during stress defense and pod development. The present study will serve as a resource for providing new strategies to enhance agronomic yield and reveal the mechanism of peanut pod development.

Peg Biology: Deciphering the Molecular Regulations Involved During Peanut Peg Development

Peanut or groundnut is one of the most important legume crops with high protein and oil content. The high nutritional qualities of peanut and its multiple usage have made it an indispensable component of our daily life, in both confectionary and therapeutic food industries. Given the socio-economic significance of peanut, understanding its developmental biology is important in providing a molecular framework to support breeding activities. In peanut, the formation and directional growth of a specialized reproductive organ called a peg, or gynophore, is especially relevant in genetic improvement. Several studies have indicated that peanut yield can be improved by improving reproductive traits including peg development. Therefore, we aim to identify unifying principles for the genetic control, underpinning molecular and physiological basis of peg development for devising appropriate strategy for peg improvement. This review discusses the current understanding of the molecular aspects of peanut peg development citing several studies explaining the key mechanisms. Deciphering and integrating recent transcriptomic, proteomic, and miRNA-regulomic studies provide a new perspective for understanding the regulatory events of peg development that participate in pod formation and thus control yield.

The characteristic of Arachis duranensis-specific genes and their potential function

Arachis species produce flowers aerially, and then grow into the ground, where they develop into fruits; a feature that is unique to Arachis species. We hypothesized that Arachis species evolved genes specifically involved in the control of aerial flowers and the formation of underground fruits. Arachis duranensis is more resistant to biotic and abiotic stressors. Here, we compared different legume species and identified Arachis duranensis-specific genes. We analyzed gene expression patterns, base substitution patterns and sequence features between genes that are conserved across legume plants and A. duranensis-specific genes. Furthermore, we tested the role of A. duranensis-specific genes during seed development, response to nematode Meloidogyne arenaria infection and drought stress. We found that A. duranensis-specific genes had characteristics of young genes. The gene expression level and breadth were lower in the A. duranensis-specific genes compared to conserved genes. The A. duranensis-specific genes had higher codon usage bias than conserved genes, and the polypeptide length and GC content at the three codon sites were lower compared to conserved genes. Of the A. duranensis-specific genes, single-copy and duplicated genes had different features. The RNA-seq result showed A. duranensis-specific genes were involved in seed development, as well as response to nematode infection and drought stress. In addition, we detected asymmetric functions in A. duranensis-specific duplicated genes in response to nematode infection and drought stress.

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

Development of chloroplast genome resources for peanut (Arachis hypogaea L.) and other species of Arachis

ABSRACT

Peanut (Arachis hypogaea L.) is an important oilseed and cash crop worldwide. Wild Arachis spp. are potental sources of novel genes for the genetic improvement of cultivated peanut. Understanding the genetic relationships with cultivated peanut is important for the efficient use of wild species in breeding programmes. However, for this genus, only a few genetic resources have been explored so far. In this study, new chloroplast genomic resources have been developed for the genus Arachis based on whole chloroplast genomes from seven species that were sequenced using next-generation sequencing technologies. The chloroplast genomes ranged in length from 156,275 to 156,395 bp, and their gene contents, gene orders, and GC contents were similar to those for other Fabaceae species. Comparative analyses among the seven chloroplast genomes revealed 643 variable sites that included 212 singletons and 431 parsimony-informative sites. We also identified 101 SSR loci and 85 indel mutation events. Thirty-seven SSR loci were found to be polymorphic by in silico comparative analyses. Eleven highly divergent DNA regions, suitable for phylogenetic and species identification, were detected in the seven chloroplast genomes. A molecular phylogeny based on the complete chloroplast genome sequences provided the best resolution of the seven Arachis species.

Mutant Transcriptome Sequencing Provides Insights into Pod Development in Peanut (Arachis hypogaea L.)

Pod size is the major yield component and a key target trait that is selected for in peanut breeding. However, although numerous quantitative trait loci (QTLs) for peanut pod size have been described, the molecular mechanisms underlying the development of this characteristic remain elusive. A peanut mutant with a narrower pod was developed in this study using ethyl methanesulfonate (EMS) mutagenesis and designated as the "pod width" mutant line (pw). The fresh pod weight of pw was only about 40% of that seen in the wild-type (WT) Zhonghua16, while the hull and seed filling of the mutant both also developed at earlier stages. Pods from both pw and WT lines were sampled 20, 40, and 60 days after flowering (DAF) and used for RNA-Seq analysis; the results revealed highly differentially expressed lignin metabolic pathway genes at all three stages, but especially at DAF 20 and DAF 40. At the same time, expression of genes related to auxin signal transduction was found to be significantly repressed during the pw early pod developmental stage. A genome-wide comparative analysis of expression profiles revealed 260 differentially expressed genes (DEGs) across all three stages, and two candidate genes, c26901_g1 (CAD) and c37339_g1 (ACS), responsible for pod width were identified by integrating expression patterns and function annotation of the common DEGs within the three stages. Taken together, the information provided in this study illuminates the processes underlying peanut pod development, and will facilitate further identification of causal genes and the development of improved peanut varieties with higher yields.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

A Developmental Transcriptome Map for Allotetraploid Arachis hypogaea

The advent of the genome sequences of Arachis duranensis and Arachis ipaensis has ushered in a new era for peanut genomics. With the goal of producing a gene atlas for cultivated peanut (Arachis hypogaea), 22 different tissue types and ontogenies that represent the full development of peanut were sequenced, including a complete reproductive series from flower to peg elongation and peg tip immersion in the soil to fully mature seed. Using a genome-guided assembly pipeline, a homeolog-specific transcriptome assembly for Arachis hypogaea was assembled and its accuracy was validated. The assembly was used to annotate 21 developmental co-expression networks as tools for gene discovery. Using a set of 8816 putative homeologous gene pairs, homeolog expression bias was documented, and although bias was mostly balanced, there were striking differences in expression bias in a tissue-specific context. Over 9000 alterative splicing events and over 6000 non-coding RNAs were further identified and profiled in a developmental context. Together, this work represents a major new resource for cultivated peanut and will be integrated into peanutbase.org as an available resource for all peanut researchers.