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

Response mechanism of carbon metabolism of Pinus massoniana to gradient high temperature and drought stress

BACKGROUND

The carbon metabolism pathway is of paramount importance for the growth and development of plants, exerting a pivotal regulatory role in stress responses. The exacerbation of drought impacts on the plant carbon cycle due to global warming necessitates comprehensive investigation into the response mechanisms of Masson Pine (Pinus massoniana Lamb.), an exemplary pioneer drought-tolerant tree, thereby establishing a foundation for predicting future forest ecosystem responses to climate change.

RESULTS

The seedlings of Masson Pine were utilized as experimental materials in this study, and the transcriptome, metabolome, and photosynthesis were assessed under varying temperatures and drought intensities. The findings demonstrated that the impact of high temperature and drought on the photosynthetic rate and transpiration rate of Masson Pine seedlings was more pronounced compared to individual stressors. The analysis of transcriptome data revealed that the carbon metabolic pathways of Masson Pine seedlings were significantly influenced by high temperature and drought co-stress, with a particular impact on genes involved in starch and sucrose metabolism. The metabolome analysis revealed that only trehalose and Galactose 1-phosphate were specifically associated with the starch and sucrose metabolic pathways. Furthermore, the trehalose metabolic heat map was constructed by integrating metabolome and transcriptome data, revealing a significant increase in trehalose levels across all three comparison groups. Additionally, the PmTPS1, PmTPS5, and PmTPPD genes were identified as key regulatory genes governing trehalose accumulation.

CONCLUSIONS

The combined effects of high temperature and drought on photosynthetic rate, transpiration rate, transcriptome, and metabolome were more pronounced than those induced by either high temperature or drought alone. Starch and sucrose metabolism emerged as the pivotal carbon metabolic pathways in response to high temperature and drought stress in Masson pine. Trehalose along with PmTPS1, PmTPS5, and PmTPPD genes played crucial roles as metabolites and key regulators within the starch and sucrose metabolism.

Molecular tagging of seed size using MITE markers in an induced large seed mutant with higher cotyledon cell size in groundnut

A large seed mutant, TG 89 having a 76.7% increment in hundred kernel weight in comparison to its parent TG 26, was isolated from an electron beam-induced mutagenized population. Studies based on environmental scanning electron microscopy of both parent and mutant revealed that the mutant seed cotyledon had significantly bigger cell size than parent. A mapping population with 122 F2 plants derived from the mutant and a distant normal seed genotype (ICGV 15007) was utilized to map the QTL associated with higher HKW. Bulk segregant analysis revealed putative association of three markers with this mutant large seed trait. Further, genotyping of F2 individuals with polymorphic markers detected 14 linkage groups with a map distance of 1053 cM. QTL analysis revealed a significant additive major QTL for the mutant large seed trait on linkage group A05 explaining 12.7% phenotypic variation for the seed size. This QTL was located between flanking markers AhTE333 and AhTE810 having a map interval of 4.7 cM which corresponds to 90.65 to 107.24 Mbp in A05 chromosome, respectively. Within this genomic fragment, an ortholog of the BIG SEEDS 1 gene was found at 102,476,137 bp. Real-time PCR revealed down-regulation of this BIG SEEDS 1 gene in the mutant indicating a loss of function mutation giving rise to a large seed phenotype. This QTL was validated in 11 advanced breeding lines having large seed size from this mutant but with varied genetic backgrounds. This validation showcased a highly promising selection accuracy of 90.9% for the marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03909-0.

Identification of QTLs associated with very-long chain fatty acid (VLCFA) content via linkage mapping and BSA-seq in peanut

KEY MESSAGE

Three major QTLs qA01, qB04.1 and qB05 for VLCFA content and their corresponding allele-specific markers will benefit peanut low VLCFA breeding, and a candidate gene Arahy.IF1JV3 was predicted. Peanut is a globally significant oilseed crop worldwide, and contains a high content (20%) of saturated fatty acid (SFA) in its seeds. As high level SFA intake in human dietary may increase the cardiovascular disease risk, reducing the SFA content in peanut is crucial for improving its nutritional quality. Half of the SFAs in peanut are very long-chain fatty acids (VLCFA), so reducing the VLCFA content is a feasible strategy to decrease the total SFA content. Luoaowan with extremely low VLCFA (4.80%) was crossed with Jihua16 (8.00%) to construct an F2:4 population. Three major QTLs including qA01, qB04.1 and qB05 for VLCFA content were detected with 4.43 ~ 14.32% phenotypic variation explained through linkage mapping. Meanwhile, three genomic regions on chromosomes B03, B04 and B05 were identified via BSA-seq approach. Two co-localized intervals on chromosomes B04 (100.10 ~ 103.97 Mb) and B05 (6.39 ~ 10.90 Mb) were identified. With markers developed based on SNP/InDel variations in qA01 between the two parents, the remaining interval was refined to 103.58 ~ 111.14 Mb. A candidate gene Arahy.IF1JV3 encoding a β-ketoacyl-CoA synthase was found in qA01, and its expression level in Luoaowan was significantly lower than that in Jihua16. Allele-specific markers targeting qA01, qB04.1 and qB05 were developed and validated in F4 population, and an elite line with high oleic, low VLCFA (5.05%) and low SFA (11.48%) contents was selected. This study initially revealed the genetic mechanism of VLCFA content, built a marker-assisted selection system for low VLCFA breeding, and provided an effective method to decrease the SFA content in peanut.

Genetic control of pod morphological traits and pod edibility in a common bean RIL population

KEY MESSAGE

QTL mapping, association analysis, and colocation study with previously reported QTL revealed three main regions controlling pod morphological traits and two loci for edible pod characteristics on the common bean chromosomes Pv01 and Pv06. Bean pod phenotype is a complex characteristic defined by the combination of different traits that determine the potential use of a genotype as a snap bean. In this study, the TUM RIL population derived from a cross between 'TU' (dry) and 'Musica' (snap) was used to investigate the genetic control of pod phenotype. The character was dissected into pod morphological traits (PMTs) and edible pod characteristics (EPC). The results revealed 35 QTL for PMTs located on seven chromosomes, suggesting a strong QTL colocation on chromosomes Pv01 and Pv06. Some QTL were colocated with previously reported QTL, leading to the mapping of 15 consensus regions associated with bean PMTs. Analysis of EPC of cooked beans revealed that two major loci with epistatic effect, located on chromosomes Pv01 and Pv06, are involved in the genetic control of this trait. An association study using a subset of the Spanish Diversity Panel (snap vs. non-snap) detected 23 genomic regions, with three regions being mapped at a position similar to those of two loci identified in the TUM population. The results demonstrated the relevant roles of Pv01 and Pv06 in the modulation of bean pod phenotype. Gene ontology enrichment analysis revealed a significant overrepresentation of genes regulating the phenylpropanoid metabolic process and auxin response in regions associated with PMTs and EPC, respectively. Both biological functions converged in the lignin biosynthetic pathway, suggesting the key role of the pathway in the genetic control of bean pod phenotype.

The determination of peanut (Arachis hypogaea L.) pod-sizes during the rapid-growth stage by phytohormones

BACKGROUND

Pod size is an important yield target trait for peanut breeding. However, the molecular mechanism underlying the determination of peanut pod size still remains unclear.

RESULTS

In this study, two peanut varieties with contrasting pod sizes were used for comparison of differences on the transcriptomic and endogenous hormonal levels. Developing peanut pods were sampled at 10, 15, 20, 25 and 30 days after pegging (DAP). Our results showed that the process of peanut pod-expansion could be divided into three stages: the gradual-growth stage, the rapid-growth stage and the slow-growth stage. Cytological analysis confirmed that the faster increase of cell-number during the rapid-growth stage was the main reason for the formation of larger pod size in Lps. Transcriptomic analyses showed that the expression of key genes related to the auxin, the cytokinin (CK) and the gibberellin (GA) were mostly up-regulated during the rapid-growth stage. Meanwhile, the cell division-related differentially expressed genes (DEGs) were mostly up-regulated at 10DAP which was consistent with the cytological-observation. Additionally, the absolute quantification of phytohormones were carried out by liquid-chromatography coupled with the tandem-mass-spectrometry (LC-MS/MS), and results supported the findings from comparative transcriptomic studies.

CONCLUSIONS

It was speculated that the differential expression levels of TAA1 and ARF (auxin-related), IPT and B-ARR (CK-related), KAO, GA20ox and GA3ox (GA-related), and certain cell division-related genes (gene-LOC112747313 and gene-LOC112754661) were important participating factors of the determination-mechanism of peanut pod sizes. These results were informative for the elucidation of the underlying regulatory network in peanut pod-growth and would facilitate further identification of valuable target genes.

Comparative Multi-Omics Analysis Reveals Lignin Accumulation Affects Peanut Pod Size

Pod size is one of the important factors affecting peanut yield. However, the metabolites relating to pod size and their biosynthesis regulatory mechanisms are still unclear. In the present study, two peanut varieties (Tif and Lps) with contrasting pod sizes were used for a comparative metabolome and transcriptome analysis. Developing peanut pods were sampled at 10, 20 and 30 days after pegging (DAP). A total of 720 metabolites were detected, most of which were lipids (20.3%), followed by phenolic acids (17.8%). There were 43, 64 and 99 metabolites identified as differentially accumulated metabolites (DAMs) at 10, 20 and 30 DAP, respectively, and flavonoids were the major DAMs between Tif and Lps at all three growth stages. Multi-omics analysis revealed that DAMs and DEGs (differentially expressed genes) were significantly enriched in the phenylpropanoid biosynthesis (ko00940) pathway, the main pathway of lignin biosynthesis, in each comparison group. The comparisons of the metabolites in the phenylpropanoid biosynthesis pathway accumulating in Tif and Lps at different growth stages revealed that the accumulation of p-coumaryl alcohol (H-monolignol) in Tif was significantly greater than that in Lps at 30 DAP. The differential expression of gene-LOC112771695, which is highly correlated with p-coumaryl alcohol and involved in the biosynthesis of monolignols, between Tif and Lps might explain the differential accumulation of p-coumaryl alcohol. The content of H-lignin in genetically diverse peanut varieties demonstrated that H-lignin content affected peanut pod size. Our findings would provide insights into the metabolic factors influencing peanut pod size and guidance for the genetic improvement of the peanut.

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.

Transcriptome analysis of pod mutant reveals plant hormones are important regulators in controlling pod size in peanut (Arachis hypogaea L.)

Pod size is an important yield-influencing trait in peanuts. It is affected by plant hormones and identifying the genes related to these hormones may contribute to pod-related trait improvements in peanut breeding programs. However, there is limited information on the molecular mechanisms of plant hormones that regulate pod size in peanuts. We identified a mutant with an extremely small pod (spm) from Yuanza 9102 (WT) by ⁶⁰Co γ-radiation mutagenesis. The length and width of the natural mature pod in spm were only 71.34% and 73.36% of those in WT, respectively. We performed comparative analyses for morphological characteristics, anatomy, physiology, and global transcriptome between spm and WT pods. Samples were collected at 10, 20, and 30 days after peg elongation into the soil, representing stages S1, S2, and S3, respectively. The differences in pod size between WT and spm were seen at stage S1 and became even more striking at stages S2 and S3. The cell sizes of the pods were significantly smaller in spm than in WT at stages S1, S2, and S3. These results suggested that reduced cell size may be one of the important contributors for the small pod in spm. The contents of indole-3-acetic acid (IAA), gibberellin (GA), and brassinosteroid (BR) were also significantly lower in spm pods than those in WT pods at all three stages. RNA-Seq analyses showed that 1,373, 8,053, and 3,358 differently expressed genes (DEGs) were identified at stages S1, S2, and S3, respectively. Functional analyses revealed that a set of DEGs was related to plant hormone biosynthesis, plant hormone signal transduction pathway, and cell wall biosynthesis and metabolism. Furthermore, several hub genes associated with plant hormone biosynthesis and signal transduction pathways were identified through weighted gene co-expression network analysis. Our results revealed that IAA, GA, and BR may be important regulators in controlling pod size by regulating cell size in peanuts. This study provides helpful information for the understanding of the complex mechanisms of plant hormones in controlling pod size by regulating the cell size in peanuts and will facilitate the improvement of peanut breeding.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development–related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development–related traits in ricebeans. A total of 51 genes encoding SCF^TIR1/AFB, Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean–derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

Understanding Omics Driven Plant Improvement and de novo Crop Domestication: Some Examples

In the current era, one of biggest challenges is to shorten the breeding cycle for rapid generation of a new crop variety having high yield capacity, disease resistance, high nutrient content, etc. Advances in the "-omics" technology have revolutionized the discovery of genes and bio-molecules with remarkable precision, resulting in significant development of plant-focused metabolic databases and resources. Metabolomics has been widely used in several model plants and crop species to examine metabolic drift and changes in metabolic composition during various developmental stages and in response to stimuli. Over the last few decades, these efforts have resulted in a significantly improved understanding of the metabolic pathways of plants through identification of several unknown intermediates. This has assisted in developing several new metabolically engineered important crops with desirable agronomic traits, and has facilitated the de novo domestication of new crops for sustainable agriculture and food security. In this review, we discuss how "omics" technologies, particularly metabolomics, has enhanced our understanding of important traits and allowed speedy domestication of novel crop plants.

Physical mapping of repetitive oligonucleotides facilitates the establishment of a genome map-based karyotype to identify chromosomal variations in peanut

BACKGROUND

Chromosomal variants play important roles in crop breeding and genetic research. The development of single-stranded oligonucleotide (oligo) probes simplifies the process of fluorescence in situ hybridization (FISH) and facilitates chromosomal identification in many species. Genome sequencing provides rich resources for the development of oligo probes. However, little progress has been made in peanut due to the lack of efficient chromosomal markers. Until now, the identification of chromosomal variants in peanut has remained a challenge.

RESULTS

A total of 114 new oligo probes were developed based on the genome-wide tandem repeats (TRs) identified from the reference sequences of the peanut variety Tifrunner (AABB, 2n = 4x = 40) and the diploid species Arachis ipaensis (BB, 2n = 2x = 20). These oligo probes were classified into 28 types based on their positions and overlapping signals in chromosomes. For each type, a representative oligo was selected and modified with green fluorescein 6-carboxyfluorescein (FAM) or red fluorescein 6-carboxytetramethylrhodamine (TAMRA). Two cocktails, Multiplex #3 and Multiplex #4, were developed by pooling the fluorophore conjugated probes. Multiplex #3 included FAM-modified oligo TIF-439, oligo TIF-185-1, oligo TIF-134-3 and oligo TIF-165. Multiplex #4 included TAMRA-modified oligo Ipa-1162, oligo Ipa-1137, oligo DP-1 and oligo DP-5. Each cocktail enabled the establishment of a genome map-based karyotype after sequential FISH/genomic in situ hybridization (GISH) and in silico mapping. Furthermore, we identified 14 chromosomal variants of the peanut induced by radiation exposure. A total of 28 representative probes were further chromosomally mapped onto the new karyotype. Among the probes, eight were mapped in the secondary constrictions, intercalary and terminal regions; four were B genome-specific; one was chromosome-specific; and the remaining 15 were extensively mapped in the pericentric regions of the chromosomes.

CONCLUSIONS

The development of new oligo probes provides an effective set of tools which can be used to distinguish the various chromosomes of the peanut. Physical mapping by FISH reveals the genomic organization of repetitive oligos in peanut chromosomes. A genome map-based karyotype was established and used for the identification of chromosome variations in peanut following comparisons with their reference sequence positions.

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.

Transcriptomic analyses reveal the expression and regulation of genes associated with resistance to early leaf spot in peanut

Objective

Early leaf spot (ELS) caused by Cercospora arachidicola (Hori) is a serious foliar disease in peanut worldwide, which causes considerable reduction of yield. Identification of resistance genes is important for both conventional and molecular breeding. Few resistance genes have been identified and the mechanism of defense responses to this pathogen remains unknown.

Results

We detected several genes involved in disease resistance to ELS through transcriptome analysis. Using RNA-seq technology, one hundred thirty-three differentially expressed genes (DEGs) were identified between resistant and susceptible lines. Among these DEGs, coiled coil-nucleotide binding-leucine rich repeat (NLR) type resistance genes were identified as duplicated R genes on the chromosome B2. Peanut phytoalexin deficient 4 (PAD4) regulator of effector-triggered immunity mediated by NLR resistance proteins and polyphenol oxidase (PPO) genes play important roles in early leaf spot resistance. Our study provides the useful information on plant response to C. arachidicola infection in peanut. The results suggest that a few major genes and several factors mediate the resistance to ELS disease, showing the characteristics of quantitative trait in defense responses.

Unraveling the mechanisms of resistance to Sclerotium rolfsii in peanut (Arachis hypogaea L.) using comparative RNA-Seq analysis of resistant and susceptible genotypes

Stem rot, a devastating fungal disease of peanut, is caused by Sclerotium rolfsii. RNA-sequencing approaches have been used to unravel the mechanisms of resistance to stem rot in peanut over the course of fungal infection in resistant (NRCG-CS85) and susceptible (TG37A) genotypes under control conditions and during the course of infection. Out of about 290 million reads, nearly 251 million (92.22%) high-quality reads were obtained and aligned to the Arachis duranensis and Arachis ipaensis genomes with the average mapping of 78.91% and 78.61%, respectively. In total, about 48.6% of genes were commonly regulated, while approximately 21.8% and 29.6% of uniquely regulated genes from A. duranensis and A. ipaensis genomes, respectively, were identified. Several annotated transcripts, such as receptor-like kinases, jasmonic acid pathway enzymes, and transcription factors (TFs), including WRKY, Zinc finger protein, and C2-H2 zinc finger, showed higher expression in resistant genotypes upon infection. These transcripts have a known role in channelizing the downstream of pathogen perception. The higher expression of WRKY transcripts might have induced the systemic acquired resistance (SAR) by the activation of the jasmonic acid defense signaling pathway. Furthermore, a set of 30 transcripts involved in the defense mechanisms were validated with quantitative real-time PCR. This study suggested PAMP-triggered immunity as a probable mechanism of resistance, while the jasmonic acid signaling pathway was identified as a possible defense mechanism in peanut. The information generated is of immense importance in developing more effective ways to combat the stem rot disease in peanut.

Transcriptome and metabolome reveal redirection of flavonoids in a white testa peanut mutant

BACKGROUND

Coat color determines both appearance and nutrient quality of peanut. White seed coat in peanut can enhance the processing efficiency and quality of peanut oil. An integrative analysis of transcriptomes, metabolomes and histocytology was performed on wsc mutant and its wild type to investigate the regulatory mechanisms underlying color pigmentation.

RESULT

Metabolomes revealed flavonoids were redirected in wsc, while multi-omics analyses of wsc mutant seeds and testae uncovered WSC influenced the flavonoids biosynthesis in testa as well as suberin formation, glycolysis, the TCA cycle and amino acid metabolism. The mutation also enhanced plant hormones synthesis and signaling. Further, co-expression analysis showed that FLS genes co-expressed with MBW complex member genes. Combining tissue expression patterns, genetic analyses, and the annotation of common DEGs for these three stages revealed that three testa specific expressed candidate genes, Araip.M7RY3, Aradu.R8PMF and Araip.MHR6K were likely responsible for the white testa phenotype. WSC might be regulated expression competition between FLS and DFR by controlling hormone synthesis and signaling as well as the MBW complex.

CONCLUSIONS

The results of this study therefore provide both candidate genes and novel approaches that can be applied to improve peanut with desirable seed coat color and flavonoid quality.

Genome-wide systematic characterization of bZIP transcription factors and their expression profiles during seed development and in response to salt stress in peanut

BACKGROUND

Plant basic leucine zipper (bZIP) transcription factors play crucial roles in plant growth, development, and abiotic stress responses. However, systematic investigation and analyses of the bZIP gene family in peanut are lacking in spite of the availability of the peanut genome sequence.

RESULTS

In this study, we identified 50 and 45 bZIP genes from Arachis duranensis and A. ipaensis genomes, respectively. Phylogenetic analysis showed that Arachis bZIP genes were classified into nine groups, and these clusters were supported by several group-specific features, including exon/intron structure, intron phases, MEME motifs, and predicted binding site structure. We also identified possible variations in DNA-binding-site specificity and dimerization properties among different Arachis bZIPs by inspecting the amino acid residues at some key sites. Our analysis of the evolutionary history analysis indicated that segmental duplication, rather than tandem duplication, contributed greatly to the expansion of this gene family, and that most Arachis bZIPs underwent strong purifying selection. Through RNA-seq and quantitative real-time PCR (qRT-PCR) analyses, the co-expressed, differentially expressed and several well-studied homologous bZIPs were identified during seed development stages in peanut. We also used qRT-PCR to explore changes in bZIP gene expression in response to salt-treatment, and many candidate bZIPs in groups A, B, and S were proven to be associated with the salt-stress response.

CONCLUSIONS

This study have conducted a genome-wide identification, characterization and expression analysis of bZIP genes in Arachis genomes. Our results provide insights into the evolutionary history of the bZIP gene family in peanut and the funcntion of Arachis bZIP genes during seed development and in response to salt stress.