Like many other crops, the cultivated peanut (Arachis hypogaea L.) is of hybrid origin and has a polyploid genome that contains essentially complete sets of chromosomes from two ancestral species. Here we report the genome sequence of peanut and show that after its polyploid origin, the genome has evolved through mobile-element activity, deletions and by the flow of genetic information between corresponding ancestral chromosomes (that is, homeologous recombination). Uniformity of patterns of homeologous recombination at the ends of chromosomes favors a single origin for cultivated peanut and its wild counterpart A. monticola. However, through much of the genome, homeologous recombination has created diversity. Using new polyploid hybrids made from the ancestral species, we show how this can generate phenotypic changes such as spontaneous changes in the color of the flowers. We suggest that diversity generated by these genetic mechanisms helped to favor the domestication of the polyploid A. hypogaea over other diploid Arachis species cultivated by humans.

Blanchability is an often overlooked, but important trait for peanut breeding. The process of blanching, removing the skin, is an important step in the processing of raw nuts for manufacturing. Under strong genetic control and requiring considerable effort to phenotype, blanchability is conducive for marker-assisted selection. We used QTL sequencing (QTL-seq) to identify two QTLs related to blanchability using previously phenotyped breeding populations. To validate the QTLs, we show that two markers can select for significantly increased blanchability in an independent recombinant inbred line (RIL) population. Two wild introgressions from Arachis cardenasii conferring strong disease resistance were segregated in the population and were found to negatively impact blanchability. Finally, we show that by utilizing highly accurate sequence analysis pipelines, low coverage sequencing can be used to genotype whole populations with increased power and precision. This study highlights the potential to mine breeding data to identify and develop useful markers for genetic improvement programs, and provide powerful tools for breeding for processing and quality traits.

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae. However, less progress has been made for A. flavus. As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus. Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Post-harvest aflatoxin contamination is a challenging issue that affects peanut quality. Aflatoxin is produced by fungi belonging to the Aspergilli group, and is known as an acutely toxic, carcinogenic and immune-suppressing class of mycotoxins. Evidence for several host genetic factors that may impact aflatoxin contamination has been reported, e.g., genes for lipoxygenase (PnLOX1 and PnLOX2/PnLOX3 that showed either positive or negative regulation with Aspergillus infection), ROS, and WRKY (highly associated with or differentially expressed upon infection of maize with A. flavus); however, their roles remain unclear. Therefore, we conducted an RNA-seq experiment to differentiate gene response to the infection by Aspergillus flavus between resistant (ICG 1471) and susceptible (Florida-07) cultivated peanut genotypes. The gene expression profiling analysis was designed to reveal differentially expressed genes in response to the infection (infected vs mock-treated seeds). In addition, the differential expression of the fungal genes was profiled. The study revealed the complexity of the interaction between the fungus and peanut seeds as expression of a large number of genes was altered including some in the process of plant defense to aflatoxin accumulation. Analysis of the experimental data with 9keggseq9, a novel designed tool for KEGG enrichment analysis, showed the importance of alpha-linolenic acid metabolism, protein processing in endoplasmic reticulum, spliceosome, and carbon fixation and metabolism pathways in conditioning resistance to aflatoxin accumulation. In addition, co-expression network analysis was carried out to reveal the correlation of gene expression among peanut and fungal genes. The results showed the importance of WRKY, TIR-NBS-LRR, ethylene, and heat shock proteins in the resistance mechanism.

Aflatoxin contamination is a major economic and food safety concern for the peanut industry that largely could be mitigated by genetic resistance. To screen peanut for aflatoxin resistance, ten genotypes were infected with a green fluorescent protein (GFP)—expressing Aspergillus flavus strain. Percentages of fungal infected area and fungal GFP signal intensity were documented by visual ratings every 8 h for 72 h after inoculation. Significant genotypic differences in fungal growth rates were documented by repeated measures and area under the disease progress curve (AUDPC) analyses. SICIA (Seed Infection Coverage and Intensity Analyzer), an image processing software, was developed to digitize fungal GFP signals. Data from SICIA image analysis confirmed visual rating results validating its utility for quantifying fungal growth. Among the tested peanut genotypes, NC 3033 and GT-C20 supported the lowest and highest fungal growth on the surface of peanut seeds, respectively. Although differential fungal growth was observed on the surface of peanut seeds, total fungal growth in the seeds was not significantly different across genotypes based on a fluorometric GFP assay. Significant differences in aflatoxin B levels were detected across peanut genotypes. ICG 1471 had the lowest aflatoxin level whereas Florida-07 had the highest. Two-year aflatoxin tests under simulated late-season drought also showed that ICG 1471 had reduced aflatoxin production under pre-harvest field conditions. These results suggest that all peanut genotypes support A. flavus fungal growth yet differentially influence aflatoxin production.

The genetic basis of general plant vigor is of major interest to food producers, yet the trait is recalcitrant to genetic mapping because of the number of loci involved, their small effects, and linkage. Observations of heterosis in many crops suggests that recessive, malfunctioning versions of genes are a major cause of poor performance, yet we have little information on the mutational spectrum underlying these disruptions. To address this question, we generated a long-read assembly of a tropicaljaponicarice (Oryza sativa) variety, Carolina Gold, which allowed us to identify structural mutations (>50 bp) and orient them with respect to their ancestral state using the outgroup,Oryza glaberrima. Supporting prior work, we find substantial genome expansion in thesativabranch. While transposable elements (TEs) account for the largest share of size variation, the majority of events are not directly TE-mediated. Tandem duplications are the most common source of insertions and are highly enriched among 50-200bp mutations. To explore the relative impact of various mutational classes on crop fitness, we then track these structural events over the last century of US rice improvement using 101 resequenced varieties. Within this material, a pattern of temporary hybridization between medium and long-grain varieties was followed by recent divergence. During this long-term selection, structural mutations that impact gene exons have been removed at a greater rate than intronic indels and single-nucleotide mutations. These results support the use ofab initioestimates of mutational burden, based on structural data, as an orthogonal predictor in genomic selection.