High oil and protein content make tetraploid peanut a leading oil and food legume. Here we report a high-quality peanut genome sequence, comprising 2.54 Gb with 20 pseudomolecules and 83,709 protein-coding gene models. We characterize gene functional groups implicated in seed size evolution, seed oil content, disease resistance and symbiotic nitrogen fixation. The peanut B subgenome has more genes and general expression dominance, temporally associated with long-terminal-repeat expansion in the A subgenome that also raises questions about the A-genome progenitor. The polyploid genome provided insights into the evolution of Arachis hypogaea and other legume chromosomes. Resequencing of 52 accessions suggests that independent domestications formed peanut ecotypes. Whereas 0.42–0.47 million years ago (Ma) polyploidy constrained genetic variation, the peanut genome sequence aids mapping and candidate-gene discovery for traits such as seed size and color, foliar disease resistance and others, also providing a cornerstone for functional genomics and peanut improvement.

Background: MicroRNAs (miRNAs) are a new class of small, endogenous RNAs that play a regulatory role in the cell by negatively affecting gene expression at the post-transcriptional level. miRNAs have been shown to control numerous genes involved in various biological and metabolic processes. There have been extensive studies on discovering miRNAs and analyzing their functions in model species, such as Arabidopsis and rice. Increasing investigations have been performed on important agricultural crops including soybean, conifers, and Phaselous vulgaris but no studies have been reported on discovering peanut miRNAs using a cloning strategy. Results: In this study, we employed the next generation high through-put Solexa sequencing technology to clone and identify both conserved and species-specific miRNAs in peanuts. Next generation high through-put Solexa sequencing showed that peanuts have a complex small RNA population and the length of small RNAs varied, 24-nt being the predominant length for a majority of the small RNAs. Combining the deep sequencing and bioinformatics, we discovered 14 novel miRNA families as well as 75 conserved miRNAs in peanuts. All 14 novel peanut miRNAs are considered to be species-specific because no homologs have been found in other plant species except ahy-miRn1, which has a homolog in soybean. qRT-PCR analysis demonstrated that both conserved and peanut-specific miRNAs are expressed in peanuts. Conclusions: This study led to the discovery of 14 novel and 22 conserved miRNA families from peanut. These results show that regulatory miRNAs exist in agronomically important peanuts and may play an important role in peanut growth, development, and response to environmental stress.

WRKY proteins are plant specific transcription factors involved in various developmental and physiological processes, especially in biotic and abiotic stress resistance. Although previous studies suggested that WRKY proteins in soybean (Glycine max var. Williams 82) involved in both abiotic and biotic stress responses, the global information of WRKY proteins in the latest version of soybean genome (Wm82.a2v1) and their response to dehydration and salt stress have not been reported. In this study, we identified 176 GmWRKY proteins from soybean Wm82.a2v1 genome. These proteins could be classified into three groups, namely group I (32 proteins), group II (120 proteins), and group III (24 proteins). Our results showed that most GmWRKY genes were located on Chromosome 6, while chromosome 11, 12, and 20 contained the least number of this gene family. More GmWRKY genes were distributed on the ends of chromosomes to compare with other regions. The cis-acting elements analysis suggested that GmWRKY genes were transcriptionally regulated upon dehydration and salt stress. RNA-seq data analysis indicated that three GmWRKY genes responded negatively to dehydration, and 12 genes positively responded to salt stress at 1, 6, and 12 h, respectively. We confirmed by qRT-PCR that the expression of GmWRKY47 and GmWRKY 58 genes was decreased upon dehydration, and the expression of GmWRKY92, 144 and 165 genes was increased under salt treatment.

DNA methylation plays important roles in genome protection, regulation of gene expression and was associated with plants development. Plant DNA methylation pattern was mediated by cytosine-5 DNA methyltransferases and demethylase. Although the genomes of AA and BB wild peanuts have been fully sequence, these two gene families have not been studied. In this study we report the identification and analysis of putative cytosine-5 DNA methyltransferases (C5-MTases) and demethylase in AA and BB wild peanuts. Cytosine-5 DNA methyltransferases in AA and BB wild peanuts could be classified in known MET, CMT and DRM2 groups based on their domain organization. This result was supported by the gene and protein structural characteristics and phylogenetic analysis. We found that some wild peanut DRM2 numbers didn’t contain UBA domain which was different from other plants such as Arabidopsis, maize, soybean. Five DNA demethylase were found in AA genome and five in BB genome. The selective pressure analysis showed that wild peanut C5-MTases gene mainly underwent purifying selection but many positive selection sites can be detected. Conversely, DNA demethylase genes mainly underwent positive selection during evolution. Additionally, the expression dynamic of cytosine-5 DNA methyltransferases and demethylase genes in different cultivated peanut tissues were analyzed. Expression result showed that cold, heat or drought stress could influence the expression level of C5-MTases and DNA demethylase genes in cultivated peanut. These results are useful for better understanding the complexity of these two gene families, and will facilitate epigenetic studies in peanut.

Summary: Whole‐genome resequencing (WGRS) of mapping populations has facilitated development of high‐density genetic maps essential for fine mapping and candidate gene discovery for traits of interest in crop species. Leaf spots, including early leaf spot (ELS) and late leaf spot (LLS), and Tomato spotted wilt virus (TSWV) are devastating diseases in peanut causing significant yield loss. We generated WGRS data on a recombinant inbred line population, developed a SNP‐based high‐density genetic map, and conducted fine mapping, candidate gene discovery and marker validation for ELS, LLS and TSWV. The first sequence‐based high‐density map was constructed with 8869 SNPs assigned to 20 linkage groups, representing 20 chromosomes, for the ‘T’ population (Tifrunner × GT‐C20) with a map length of 3120 cM and an average distance of 1.45 cM. The quantitative trait locus (QTL) analysis using high‐density genetic map and multiple season phenotyping data identified 35 main‐effect QTLs with phenotypic variation explained (PVE) from 6.32% to 47.63%. Among major‐effect QTLs mapped, there were two QTLs for ELS on B05 with 47.42% PVE and B03 with 47.38% PVE, two QTLs for LLS on A05 with 47.63% and B03 with 34.03% PVE and one QTL for TSWV on B09 with 40.71% PVE. The epistasis and environment interaction analyses identified significant environmental effects on these traits. The identified QTL regions had disease resistance genes including R‐genes and transcription factors. KASP markers were developed for major QTLs and validated in the population and are ready for further deployment in genomics‐assisted breeding in peanut.

WRKY, an important transcription factor family, is widely distributed in the plant kingdom. Many reports focused on analysis of phylogenetic relationship and biological function of WRKY protein at the whole genome level in different plant species. However, little is known about WRKY proteins in the genome of Arachis species and their response to salicylic acid (SA) and jasmonic acid (JA) treatment. In this study, we identified 77 and 75 WRKY proteins from the two wild ancestral diploid genomes of cultivated tetraploid peanut, Arachis duranensis and Arachis ipaënsis, using bioinformatics approaches. Most peanut WRKY coding genes were located on A. duranensis chromosome A6 and A. ipaënsis chromosome B3, while the least number of WRKY genes was found in chromosome 9. The WRKY orthologous gene pairs in A. duranensis and A. ipaënsis chromosomes were highly syntenic. Our analysis indicated that segmental duplication events played a major role in AdWRKY and AiWRKY genes, and strong purifying selection was observed in gene duplication pairs. Furthermore, we translate the knowledge gained from the genome-wide analysis result of wild ancestral peanut to cultivated peanut to reveal that gene activities of specific cultivated peanut WRKY gene were changed due to SA and JA treatment. Peanut WRKY7, 8 and 13 genes were down-regulated, whereas WRKY1 and 12 genes were up-regulated with SA and JA treatment. These results could provide valuable information for peanut improvement.