Dedicated Industrial Oilseed Crops as Metabolic Engineering Platforms for Sustainable Industrial Feedstock Production

Feedstocks for industrial applications ranging from polymers to lubricants are largely derived from petroleum, a non-renewable resource. Vegetable oils with fatty acid structures and storage forms tailored for specific industrial uses offer renewable and potentially sustainable sources of petrochemical-type functionalities. A wide array of industrial vegetable oils can be generated through biotechnology, but will likely require non-commodity oilseed platforms dedicated to specialty oil production for commercial acceptance. Here we show the feasibility of three Brassicaceae oilseeds crambe, camelina, and carinata, none of which are widely cultivated for food use, as hosts for complex metabolic engineering of wax esters for lubricant applications. Lines producing wax esters >20% of total seed oil were generated for each crop and further improved for high temperature oxidative stability by down-regulation of fatty acid polyunsaturation. Field cultivation of optimized wax ester-producing crambe demonstrated commercial utility of these engineered crops and a path for sustainable production of other industrial oils in dedicated specialty oilseeds.

Combinatorial Effects of Fatty Acid Elongase Enzymes on Nervonic Acid Production in Camelina sativa

Very long chain fatty acids (VLCFAs) with chain lengths of 20 carbons and longer provide feedstocks for various applications; therefore, improvement of VLCFA contents in seeds has become an important goal for oilseed enhancement. VLCFA biosynthesis is controlled by a multi-enzyme protein complex referred to as fatty acid elongase, which is composed of β-ketoacyl-CoA synthase (KCS), β-ketoacyl-CoA reductase (KCR), β-hydroxyacyl-CoA dehydratase (HCD) and enoyl reductase (ECR). KCS has been identified as the rate-limiting enzyme, but little is known about the involvement of other three enzymes in VLCFA production. Here, the combinatorial effects of fatty acid elongase enzymes on VLCFA production were assessed by evaluating the changes in nervonic acid content. A KCS gene from Lunaria annua (LaKCS) and the other three elongase genes from Arabidopsis thaliana were used for the assessment. Five seed-specific expressing constructs, including LaKCS alone, LaKCS with AtKCR, LaKCS with AtHCD, LaKCS with AtECR, and LaKCS with AtKCR and AtHCD, were transformed into Camelina sativa. The nervonic acid content in seed oil increased from null in wild type camelina to 6-12% in LaKCS-expressing lines. However, compared with that from the LaKCS-expressing lines, nervonic acid content in mature seeds from the co-expressing lines with one or two extra elongase genes did not show further increases. Nervonic acid content from LaKCS, AtKCR and AtHCD co-expressing line was significantly higher than that in LaKCS-expressing line during early seed development stage, while the ultimate nervonic acid content was not significantly altered. The results from this study thus provide useful information for future engineering of oilseed crops for higher VLCFA production.

Development of a High-Density Genetic Map Based on Specific Length Amplified Fragment Sequencing and Its Application in Quantitative Trait Loci Analysis for Yield-Related Traits in Cultivated Peanut

High-density genetic maps (HDGMs) are very useful for genomic studies and quantitative trait loci (QTL) mapping. However, the low frequency of DNA polymorphisms in peanut has limited the quantity of available markers and hindered the construction of a HDGM. This study generated a peanut genetic map with the highest number of high-quality SNPs based on specific locus amplified fragment sequencing (SLAF-seq) technology and a newly constructed RIL population ("ZH16" × "sd-H1"). The constructed HDGM included 3,630 SNP markers belonging to 2,636 bins on 20 linkage groups (LGs), and it covers 2,098.14 cM in length, with an average marker distance of 0.58 cM. This HDGM was applied for the following collinear comparison, scaffold anchoring and analysis of genomic characterization including recombination rates and segregation distortion in peanut. For QTL mapping of investigated 14 yield-related traits, a total of 62 QTLs were detected on 12 chromosomes across 3 environments, and the co-localization of QTLs was observed for these traits which were significantly correlated on phenotype. Two stable co-located QTLs for seed- and pod-related traits were significantly identified in the chromosomal end of B06 and B07, respectively. The construction of HDGM and QTL analysis for yield-related traits in this study provide useful information for fine mapping and functional analysis of genes as well as molecular marker-assisted breeding.

Integrated Transcriptomics and Metabolomics Analysis Reveal Key Metabolism Pathways Contributing to Cold Tolerance in Peanut

Low temperature (non-freezing) is one of the major limiting factors in peanut (Arachis hypogaea L.) growth, yield, and geographic distribution. Due to the complexity of cold-resistance trait in peanut, the molecular mechanism of cold tolerance and related gene networks were largely unknown. In this study, metabolomic analysis of two peanut cultivars subjected to chilling stress obtained a set of cold-responsive metabolites, including several carbohydrates and polyamines. These substances showed a higher accumulation pattern in cold-tolerant variety SLH than cold-susceptible variety ZH12 under cold stress, indicating their importance in protecting peanut from chilling injuries. In addition, 3,620 cold tolerance genes (CTGs) were identified by transcriptome sequencing, and the CTGs were most significantly enriched in the "phenylpropanoid biosynthesis" pathway. Two vital modules and several novel hub genes were obtained by weighted gene co-expression network analysis (WGCNA). Several key genes involved in soluble sugar, polyamine, and G-lignin biosynthetic pathways were substantially higher and/or responded more quickly in SLH (cold tolerant) than ZH12 (cold susceptible) under low temperature, suggesting they might be crucial contributors during the adaptation of peanut to low temperature. These findings will not only provide valuable resources for study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

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.

Functional Biology and Molecular Mechanisms of Host-Pathogen Interactions for Aflatoxin Contamination in Groundnut (Arachis hypogaea L.) and Maize (Zea mays L.)

Aflatoxins are secondary metabolites produced by soilborne saprophytic fungus Aspergillus flavus and closely related species that infect several agricultural commodities including groundnut and maize. The consumption of contaminated commodities adversely affects the health of humans and livestock. Aflatoxin contamination also causes significant economic and financial losses to producers. Research efforts and significant progress have been made in the past three decades to understand the genetic behavior, molecular mechanisms, as well as the detailed biology of host-pathogen interactions. A range of omics approaches have facilitated better understanding of the resistance mechanisms and identified pathways involved during host-pathogen interactions. Most of such studies were however undertaken in groundnut and maize. Current efforts are geared toward harnessing knowledge on host-pathogen interactions and crop resistant factors that control aflatoxin contamination. This study provides a summary of the recent progress made in enhancing the understanding of the functional biology and molecular mechanisms associated with host-pathogen interactions during aflatoxin contamination in groundnut and maize.

Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.)

The genetic architecture determinants of yield traits in peanut (Arachis hypogaea L.) are poorly understood. In the present study, an effort was made to map quantitative trait loci (QTLs) for yield traits using recombinant inbred lines (RIL). A genetic linkage map was constructed containing 609 loci, covering a total of 1557.48 cM with an average distance of 2.56 cM between adjacent markers. The present map exhibited good collinearity with the physical map of diploid species of Arachis. Ninety-two repeatable QTLs were identified for 11 traits including height of main stem, total branching number, and nine pod- and seed-related traits. Of the 92 QTLs, 15 QTLs were expressed across three environments and 65 QTLs were newly identified. Twelve QTLs for the height of main stem and the pod- and seed-related traits explaining more than 10 % of phenotypic variation showed a great potential for marker-assisted selection in improving these traits. The trait-by-trait meta-analysis revealed 33 consensus QTLs. The consensus QTLs and other QTLs were further integrated into 29 pleiotropic unique QTLs with the confidence interval of 1.86 cM on average. The significant co-localization of QTLs was consistent with the significant phenotypic correlations among these traits. The complexity of the genetic architecture of yield traits was demonstrated. The present QTLs for pod- and seed-related traits could be the most fundamental genetic factors contributing to the yield traits in peanut. The results provide a good foundation for fine mapping, cloning and designing molecular breeding of favorable genes in peanut.

Genome sequencing and comparative genomic analysis of highly and weakly aggressive strains of Sclerotium rolfsii, the causal agent of peanut stem rot

BACKGROUND

Stem rot caused by Sclerotium rolfsii is a very important soil-borne disease of peanut. S. rolfsii is a necrotrophic plant pathogenic fungus with an extensive host range and worldwide distribution. It can infect peanut stems, roots, pegs and pods, leading to varied yield losses. S. rolfsii strains GP3 and ZY collected from peanut in different provinces of China exhibited a significant difference in aggressiveness on peanut plants by artificial inoculation test. In this study, de-novo genome sequencing of these two distinct strains was performed aiming to reveal the genomic basis of difference in aggressiveness.

RESULTS

Scleotium rolfsii strains GP3 and ZY, with weak and high aggressiveness on peanut plants, exhibited similar growth rate and oxalic acid production in laboratory. The genomes of S. rolfsii strains GP3 and ZY were sequenced by Pacbio long read technology and exhibited 70.51 Mb and 70.61 Mb, with contigs of 27 and 23, and encoded 17,097 and 16,743 gene models, respectively. Comparative genomic analysis revealed that the pathogenicity-related gene repertoires, which might be associated with aggressiveness, differed between GP3 and ZY. There were 58 and 45 unique pathogen-host interaction (PHI) genes in GP3 and ZY, respectively. The ZY strain had more carbohydrate-active enzymes (CAZymes) in its secretome than GP3, especially in the glycoside hydrolase family (GH), the carbohydrate esterase family (CBM), and the polysaccharide lyase family (PL). GP3 and ZY also had different effector candidates and putative secondary metabolite synthetic gene clusters. These results indicated that differences in PHI, secreted CAZymes, effectors and secondary metabolites may play important roles in aggressive difference between these two strains.

CONCLUSIONS

The data provided a further understanding of the S. rolfsii genome. Genomic comparison provided clues to the difference in aggressiveness of S. rolfsii strains.

Identification of genomic regions and diagnostic markers for resistance to aflatoxin contamination in peanut (Arachis hypogaea L.)

Background

Aflatoxin contamination caused by Aspergillus flavus is a major constraint to peanut industry worldwide due to its toxicological effects to human and animals. Developing peanut varieties with resistance to seed infection and/or aflatoxin accumulation is the most effective and economic strategy for reducing aflatoxin risk in food chain. Breeding for resistance to aflatoxin in peanut is a challenging task for breeders because the genetic basis is still poorly understood. To identify the quantitative trait loci (QTLs) for resistance to aflatoxin contamination in peanut, a recombinant inbred line (RIL) population was developed from crossing Zhonghua 10 (susceptible) with ICG 12625 (resistant). The percent seed infection index (PSII), the contents of aflatoxin B₁ (AFB₁) and aflatoxin B₂ (AFB₂) of RILs were evaluated by a laboratory kernel inoculation assay.

Results

Two QTLs were identified for PSII including one major QTL with 11.32–13.00% phenotypic variance explained (PVE). A total of 12 QTLs for aflatoxin accumulation were detected by unconditional analysis, and four of them (qAFB1A07 and qAFB1B06.1 for AFB₁, qAFB2A07 and qAFB2B06 for AFB₂) exhibited major and stable effects across multiple environments with 9.32–21.02% PVE. Furthermore, not only qAFB1A07 and qAFB2A07 were co-localized in the same genetic interval on LG A07, but qAFB1B06.1 was also co-localized with qAFB2B06 on LG B06. Conditional QTL mapping also confirmed that there was a strong interaction between resistance to AFB₁ and AFB₂ accumulation. Genotyping of RILs revealed that qAFB1A07 and qAFB1B06.1 interacted additively to improve the resistance to both AFB₁ and AFB₂ accumulation. Additionally, validation of the two markers was performed in diversified germplasm collection and four accessions with resistance to aflatoxin accumulation were identified.

Conclusions

Single major QTL for resistance to PSII and two important co-localized intervals associated with major QTLs for resistance to AFB₁ and AFB₂. Combination of these intervals could improve the resistance to aflatoxin accumulation in peanut. SSR markers linked to these intervals were identified and validated. The identified QTLs and associated markers exhibit potential to be applied in improvement of resistance to aflatoxin contamination.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0734-z) contains supplementary material, which is available to authorized users.

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.

Genome-Wide Identification of Peanut KCS Genes Reveals That AhKCS1 and AhKCS28 Are Involved in Regulating VLCFA Contents in Seeds

The peanut (Arachis hypogaea L.) is an important oilseed crop worldwide. Compared to other common edible vegetable oils, peanut oil contains a higher content of saturated fatty acids (SFAs), approximately 20-40% of which are very long chain fatty acids (VLCFAs). To understand the basis for this oil profile, we interrogated genes for peanut β-ketoacyl-CoA synthase (KCS), which is known to be a key enzyme in VLCFA biosynthesis. A total of 30 AhKCS genes were identified in the assembled genome of the peanut. Based on transcriptome data, nine AhKCS genes with high expression levels in developing seeds were cloned and expressed in yeast. All these AhKCSs could produce VLCFAs but result in different profiles, indicating that the AhKCSs catalyzed fatty acid elongation with different substrate specificities. Expression level analysis of these nine AhKCS genes was performed in developing seeds from six peanut germplasm lines with different VLCFA contents. Among these genes, the expression levels of AhKCS1 or AhKCS28 were, 4-10-fold higher than that of any other AhKCS. However, only the expression levels of AhKCS1 and AhKCS28 were significantly and positively correlated with the VLCFA content, suggesting that AhKCS1 and AhKCS28 were involved in the regulation of VLCFA content in the peanut seed. Further subcellular localization analysis indicated that AhKCS1 and AhKCS28 were located at the endoplasmic reticulum (ER). Overexpression of AhKCS1 or AhKCS28 in Arabidopsis increased the contents of VLCFAs in the seed, especially for very long chain saturated fatty acids (VLCSFAs). Taken together, this study suggests that AhKCS1 and AhKCS28 could be key genes in regulating VLCFA biosynthesis in the seed, which could be applied to improve the health-promoting and nutritional qualities of the peanut.

The conserved transcription factors, MYB115 and MYB118, control expression of the newly evolved benzoyloxy glucosinolate pathway in Arabidopsis thaliana

The evolution of plant metabolic diversity is largely driven by gene duplication and ensuing sub-functionalization and/or neo-functionalization to generate new enzymatic activities. However, it is not clear whether the transcription factors (TFs) regulating these new enzyme encoding genes were required to co-evolve with these genes in a similar fashion or if these new genes can be captured by existing conserved TFs to provide the appropriate expression pattern. In this study, we found two conserved TFs, MYB115, and MYB118, co-expressed with the key enzyme encoding genes in the newly evolved benzoyloxy glucosinolate (GLS) pathway. These TFs interacted with the promoters of the GLS biosynthetic genes and negatively influenced their expression. Similarly, the GLS profiles of these two TFs knockouts showed that they influenced the aliphatic GLS accumulation within seed, leaf and flower, while they mainly expressed in seeds. Further studies indicated that they are functionally redundant and epistatically interact to control the transcription of GLS genes. Complementation study confirmed their roles in regulating the aliphatic GLS biosynthesis. These results suggest that the newly evolved enzyme encoding genes for novel metabolites can be regulated by conserved TFs, which helps to improve our model for newly evolved genes regulation.

Transcriptomic profiling reveals pigment regulation during peanut testa development

Although peanut (Arachis hypogaea L.) is one of the most important edible oil crops globally, pigments present in the testa influence both the processing efficiency and the quality of the oil. In peanut, polymeric phenolic compounds are present in the episperm rather than in the endothelium and their levels increase during ripening; therefore, to better understand testa development, and especially the accumulation of pigments, RNA-Seq was applied to elucidate the mechanisms underlying the regulation of peanut testae at three different developmental stages (i.e., at 20 days after flowering - 20DAF - and at 40DAF and 60DAF). A total of 5452 differentially expressed unigenes (DEGs) were obtained encompassing these three stages; comparative results showed that phenylpropanoid biosynthesis, phenylalanine metabolism, flavonoid biosynthesis, and plant hormone signal transduction comprised the principal KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways expressed during peanut testa development. Further studies revealed that the expression patterns of the flavonoid biosynthesis pathway genes PAL, C4H, CHS, and CHI (early biosynthetic genes - EBGs) were consistent with the accumulation of testa pigments. Thus, the results of this study demonstrate that EBGs, as well as the homologs of AtMYB111 (i.e., c35101_g4 and c37398_g2), are likely the principal regulators of testa pigment accumulation; the gene database assembled here is therefore a sequencing resource for future research and provides a foundation for understanding the regulation of pink testa pigmentation in peanuts.

Key Regulators of Sucrose Metabolism Identified through Comprehensive Comparative Transcriptome Analysis in Peanuts

Sucrose content is a crucial indicator of quality and flavor in peanut seed, and there is a lack of clarity on the molecular basis of sucrose metabolism in peanut seed. In this context, we performed a comprehensive comparative transcriptome study on the samples collected at seven seed development stages between a high-sucrose content variety (ICG 12625) and a low-sucrose content variety (Zhonghua 10). The transcriptome analysis identified a total of 8334 genes exhibiting significantly different abundances between the high- and low-sucrose varieties. We identified 28 differentially expressed genes (DEGs) involved in sucrose metabolism in peanut and 12 of these encoded sugars will eventually be exported transporters (SWEETs). The remaining 16 genes encoded enzymes, such as cell wall invertase (CWIN), vacuolar invertase (VIN), cytoplasmic invertase (CIN), cytosolic fructose-bisphosphate aldolase (FBA), cytosolic fructose-1,6-bisphosphate phosphatase (FBP), sucrose synthase (SUS), cytosolic phosphoglucose isomerase (PGI), hexokinase (HK), and sucrose-phosphate phosphatase (SPP). The weighted gene co-expression network analysis (WGCNA) identified seven genes encoding key enzymes (CIN, FBA, FBP, HK, and SPP), three SWEET genes, and 90 transcription factors (TFs) showing a high correlation with sucrose content. Furthermore, upon validation, six of these genes were successfully verified as exhibiting higher expression in high-sucrose recombinant inbred lines (RILs). Our study suggested the key roles of the high expression of SWEETs and enzymes in sucrose synthesis making the genotype ICG 12625 sucrose-rich. This study also provided insights into the molecular basis of sucrose metabolism during seed development and facilitated exploring key candidate genes and molecular breeding for sucrose content in peanuts.

Comprehensive evaluation of Chinese peanut mini-mini core collection and QTL mapping for aflatoxin resistance

BACKGROUND

Aflatoxin contamination caused by Aspergillus fungi has been a serious factor affecting food safety of peanut (Arachis hypogaea L.) because aflatoxins are highly harmful for human and animal health. As three mechanisms of resistance to aflatoxin in peanut including shell infection resistance, seed infection resistance and aflatoxin production resistance exist among naturally evolved germplasm stocks, it is highly crucial to pyramid these three resistances for promoting peanut industry development and protecting consumers' health. However, less research effort has been made yet to investigate the differentiation and genetic relationship among the three resistances in diversified peanut germplasm collections.

RESULTS

In this study, the Chinese peanut mini-mini core collection selected from a large basic collection was systematically evaluated for the three resistances against A. flavus for the first time. The research revealed a wide variation among the diversified peanut accessions for all the three resistances. Totally, 14 resistant accessions were identified, including three with shell infection resistance, seven with seed infection resistance and five with aflatoxin production resistance. A special accession, Zh.h1312, was identified with both seed infection and aflatoxin production resistance. Among the five botanic types of A. hypogaea, the var. vulgaris (Spanish type) belonging to subspecies fastigiata is the only one which possessed all the three resistances. There was no close correlation between shell infection resistance and other two resistances, while there was a significant positive correlation between seed infection and toxin production resistance. All the three resistances had a significant negative correlation with pod or seed size. A total of 16 SNPs/InDels associated with the three resistances were identified through genome-wide association study (GWAS). Through comparative analysis, Zh.h1312 with seed infection resistance and aflatoxin production resistance was also revealed to possess all the resistance alleles of associated loci for seed infection index and aflatoxin content.

CONCLUSIONS

This study provided the first comprehensive understanding of differentiation of aflatoxin resistance in diversified peanut germplasm collection, and would further contribute to the genetic enhancement for resistance to aflatoxin contamination.

Genome-wide analysis of UDP-glycosyltransferase gene family and identification of a flavonoid 7-O-UGT (AhUGT75A) enhancing abiotic stress in peanut (Arachis hypogaea L.)

BACKGROUND

Glycosylation, catalyzed by UDP-glycosyltransferase (UGT), was important for enhancing solubility, bioactivity, and diversity of flavonoids. Peanut (Arachis hypogaea L.) is an important oilseed and cash crop worldwide. In addition to provide high quality of edible oils and proteins, peanut seeds contain a rich source of flavonoid glycosides that benefit human health. However, information of UGT gene family was quite limited in peanut.

RESULTS

In present study, a total of 267 AhUGTs clustered into 15 phylogenetic groups were identified in peanut genome. Group I has greatly expanded to contain the largest number of AhUGT genes. Segmental duplication was the major driving force for AhUGT gene family expansion. Transcriptomic analysis of gene expression profiles in various tissues and under different abiotic stress treatments indicated AhUGTs were involved in peanut growth and abiotic stress response. AhUGT75A (UGT73CG33), located in mitochondria, was characterized as a flavonoid 7-O-UGT by in vitro enzyme assays. The transcript level of AhUGT75A was strongly induced by abiotic stress. Overexpression of AhUGT75A resulted in accumulating less amount of malondialdehyde (MDA) and superoxide, and enhancing tolerance against drought and/or salt stress in transgenic Arabidopsis. These results indicated AhUGT75A played important roles in conferring abiotic stress tolerance through reactive oxygen species scavenging.

CONCLUSIONS

Our research only not provides valuable information for functional characterization of UGTs in peanut, but also gives new insights into potential applications in breeding new cultivars with both desirable stress tolerance and health benefits.

Construction of ddRADseq-Based High-Density Genetic Map and Identification of Quantitative Trait Loci for Trans-resveratrol Content in Peanut Seeds

Resveratrol (trans-3,4',5-trihydroxystilbene) is a natural stilbene phytoalexin which is also found to be good for human health. Cultivated peanut (Arachis hypogaea L.), a worldwide important legume crop, is one of the few sources of human's dietary intake of resveratrol. Although the variations of resveratrol contents among peanut varieties were observed, the variations across environments and its underlying genetic basis were poorly investigated. In this study, the resveratrol content in seeds of a recombination inbred line (RIL) population (Zhonghua 6 × Xuhua 13, 186 progenies) were quantified by high performance liquid chromatography (HPLC) method across four environments. Genotypes, environments and genotype × environment interactions significantly influenced the resveratrol contents in the RIL population. A total of 8,114 high-quality single nucleotide polymorphisms (SNPs) were identified based on double-digest restriction-site-associated DNA sequencing (ddRADseq) reads. These SNPs were clustered into bins using a reference-based method, which facilitated the construction of high-density genetic map (2,183 loci with a total length of 2,063.55 cM) and the discovery of several chromosome translocations. Through composite interval mapping (CIM), nine additive quantitative trait loci (QTL) for resveratrol contents were identified on chromosomes A01, A07, A08, B04, B05, B06, B07, and B10 with 5.07-8.19% phenotypic variations explained (PVE). Putative genes within their confidential intervals might play roles in diverse primary and secondary metabolic processes. These results laid a foundation for the further genetic dissection of resveratrol content as well as the breeding and production of high-resveratrol peanuts.

High-Density Genetic Linkage Map Construction Using Whole-Genome Resequencing for Mapping QTLs of Resistance to Aspergillus flavus Infection in Peanut

The cultivated peanut (Arachis hypogaea L.), which is rich in edible oil and protein, is widely planted around the world as an oil and cash crop. However, aflatoxin contamination seriously affects the quality safety of peanuts, hindering the development of the peanut industry and threatening the health of consumers. Breeding peanut varieties with resistance to Aspergillus flavus infection is important for the control of aflatoxin contamination, and understanding the genetic basis of resistance is vital to its genetic enhancement. In this study, we reported the quantitative trait locus (QTL) mapping of resistance to A. flavus infection of a well-known resistant variety, J11. A mapping population consisting of 200 recombinant inbred lines (RILs) was constructed by crossing a susceptible variety, Zhonghua 16, with J11. Through whole-genome resequencing, a genetic linkage map was constructed with 2,802 recombination bins and an average inter-bin distance of 0.58 cM. Combined with phenotypic data of an infection index in 4 consecutive years, six novel resistant QTLs with 5.03-10.87% phenotypic variances explained (PVE) were identified on chromosomes A05, A08, B01, B03, and B10. The favorable alleles of five QTLs were from J11, while that of one QTL was from Zhonghua 16. The combination of these favorable alleles significantly improved resistance to A. flavus infection. These results could contribute greatly to the understanding of the genetic basis of A. flavus resistance and could be meaningful in the improvement of further resistance in peanuts.

Genome-wide identification of peanut PIF family genes and their potential roles in early pod development

Peanut is typically a geocarpic plant. The developing gynophore ('peg') in air could not swell normally until it buries into soil, indicating light-to-dark conversion is necessary for early pod development in peanut. As the subfamily of basic helix-loop-helix (bHLH) transcription factors, phytochrome interacting factors (PIFs) are key regulators involved in light signaling pathways, and play crucial roles in plant growth and development. In the current study, a total of 14 AhPIFs were identified in cultivated peanut genome (Arachis hypogaea L., AABB), while seven AdPIFs and six AiPIFs were discovered in the two wild diploids (A. duranensis (AA), A. ipaensis (BB)) respectively. Phylogenetic analysis revealed that peanut PIFs were clustered into four distinct clades, and members within the same subgroup had conserved motifs and displayed similar exon-intron distribution patterns. Gene synteny analysis indicated most of the PIFs exhibit one-to-one homology relationship between AA and BB subgenome in A. hypogaea, as well as among the three peanut species. Gene duplication detection showed that segmental duplication and purifying selection contributed to the expansion and evolution of peanut PIF gene family. Transcript profiles combined with subcellular localization analysis suggested AhPIF3A4 and AhPIF3B4 may possibly be involved in regulation of peanut early pod development. This study could further facilitate functional characterization of PIFs in peanut and other legumes.

Genetic dissection of fatty acid components in the Chinese peanut (Arachis hypogaea L.) mini-core collection under multi-environments

Peanut (Arachis hypogaea L.) is an important source of edible oil and protein for human nutrition. The quality of peanut seed oil is mainly determined by the composition of fatty acids, especially the contents of oleic acid and linoleic acid. Improving the composition of fatty acids in the seed oil is one of the main objectives for peanut breeding globally. To uncover the genetic basis of fatty acids and broaden the genetic variation in future peanut breeding programs, this study used genome-wide association studies (GWAS) to identify loci associated with target traits and developed diagnostic marker. The contents of eight fatty acid components of the Chinese peanut mini-core collection were measured under four environments. Using the phenotypic information and over one hundred thousand single nucleotide polymorphisms (SNPs), GWAS were conducted to investigate the genetics basis of fatty acids under multi-environments. Overall, 75 SNPs were identified significant trait associations with fatty acid components. Nineteen associations were repeatedly identified in multiple environments, and 13 loci were co-associated with two or three traits. Three stable major associated loci were identified, including two loci for oleic acid and linoleic acid on chromosome A09 [mean phenotypic variation explained (PVE): 38.5%, 10.35%] and one for stearic acid on B06 (mean PVE: 23%). According to functional annotations, 21 putative candidate genes related to fatty acid biosynthesis were found underlying the three associations. The allelic effect of SNP A09-114690064 showed that the base variation was highly correlated with the phenotypic variation of oleic acid and linoleic acid contents, and a cost-effective Kompetitive allele-Specific PCR (KASP) diagnostic marker was developed. Furthermore, the SNP A09-114690064 was found to change the cis-element CAAT (-) in the promoter of ahFAD2A to YACT (+), leading dozens of times higher expression level. The enhancer-like activity of ahFAD2A promoter was identified that was valuable for enriching the regulation mechanism of ahFAD2A. This study improved our understanding on the genetic architecture of fatty acid components in peanut, and the new effective diagnostic marker would be useful for marker-assisted selection of high-oleic peanut breeding.

Identification of a stable major sucrose-related QTL and diagnostic marker for flavor improvement in peanut

An InDel marker closely linked with a major and stable quantitative trait locus (QTL) on chromosome A08, qSUCA08.2, controlling sucrose content will benefit peanut flavor improvement. Sucrose is the main soluble sugar in mature peanut kernel, and its content is a key determinant of flavor. However, the genetic basis of sucrose content in peanut remains poorly understood, which limits the progress of flavor improvement. In the present study, two genomic regions (qSUCA08a and qSUCB06a) for sucrose content on chromosomes A08 and B06 were identified by QTL-seq in a RIL population derived from a cross between Zhonghua 10 and ICG 12625. In the interval of qSUCB06a, QTL qSUCB06.2 was detected through QTL mapping in a single environment. The qSUCA08a was further dissected into 3 adjacent genomic regions using linkage analysis including a major QTL qSUCA08.2 explaining 5.43-17.84% phenotypic variation across five environments. A 61-bp insertion at position 35,099,320 in the higher sucrose parent ICG 12625 was found in qSUCA08.2. An InDel marker SUC.InDel.A08 based on the insertion/deletion polymorphism was developed and validated within a natural population containing 172 peanut cultivars in two environments. The mean sucrose content of 93 cultivars with ICG 12625 allele was significantly higher than that of 79 cultivars with Zhonghua 10 allele. The qSUCA08.2 corresponding to a 2.11 Mb interval harbored 110 genes. Among these genes, a total of 19 genes were considered as candidate genes including 5 non-synonymous mutation genes and 14 differentially expressed genes during seed development. These results provide new insights into the genetic basis of sucrose regulation in peanut and benefit the breeding program for developing new varieties with excellent flavor.

Genome-wide expression quantitative trait locus analysis in a recombinant inbred line population for trait dissection in peanut

The transcriptome connects genome to the gene function and ultimate phenome in biology. So far, transcriptomic approach was not used in peanut for performing trait mapping in bi-parental populations. In this research, we sequenced the whole transcriptome in immature seeds in a peanut recombinant inbred line (RIL) population and explored thoroughly the landscape of transcriptomic variations and its genetic basis. The comprehensive analysis identified total 49 691 genes in RIL population, of which 92 genes followed a paramutation-like expression pattern. Expression quantitative trait locus (eQTL) analysis identified 1207 local eQTLs and 15 837 distant eQTLs contributing to the whole-genome transcriptomic variation in peanut. There were 94 eQTL hot spot regions detected across the genome with the dominance of distant eQTL. By integrating transcriptomic profile and annotation analyses, we unveiled a putative candidate gene and developed a linked marker InDel02 underlying a major QTL responsible for purple testa colour in peanut. Our result provided a first understanding of genetic basis of whole-genome transcriptomic variation in peanut and illustrates the potential of the transcriptome-aid approach in dissecting important traits in non-model plants.

Fine-mapping and gene candidate analysis for AhRt1, a major dominant locus responsible for testa color in cultivated peanut

AhRt1 controlling red testa color in peanut was fine-mapped to an interval of 580 kb on chromosome A03, and one gene encoding bHLH transcriptional factor was identified as the putative candidate gene. Peanut with red testa has higher nutritional and economic value than the traditional pink testa varieties. Identification of genes controlling red testa color will accelerate the breeding program and facilitate uncovering the genetic mechanism. In this study, in order to identify gene underlying the red testa color in peanut, a F₂ population derived from a cross between a pink testa peanut variety "Fuhua 8" and a red testa variety "Quanhonghua 1" was constructed. The genetic analysis for the F₂ population revealed that the red testa color was controlled by one single dominant locus. This locus, named as AhRt1 (Arachis hypogaea Red Testa 1), was preliminary identified in chromosome A03 by BSA-sequencing analysis. Using a segregation mapping population, AhRt1 was fine-mapped to a 580-kb genomic region by substitution mapping strategy. Gene candidate analysis suggested that one predicted gene encoding bHLH transcriptional factor may be the possible candidate gene for AhRt1. A diagnostic marker closely linked to candidate gene has been developed for validating the fine-mapping result in different populations and peanut germplasm. Our findings will benefit the breeding program for developing new varieties with red testa color and laid foundation for map-based cloning gene responsible for red testa in peanut.

Red fluorescence protein (DsRed2) promotes the screening efficiency in peanut genetic transformation

Peanut (Arachis hypogaea L.), one of the leading oilseed crops worldwide, is an important source of vegetable oil, protein, minerals and vitamins. Peanut is widely cultivated in Asia, Africa and America, and China is the largest producer and consumer of peanut. Genetic engineering has shown great potential to alter the DNA makeup of an organism which is largely hindered by the low transformation and screening efficiency including in peanut. DsRed2 is a reporter gene widely utilized in genetic transformation to facilitate the screening of transformants, but never used in peanut genetic transformation. In this study, we have demonstrated the potential of the red fluorescence protein DsRed2 as a visual reporter to improve screening efficiency in peanut. DsRed2 was firstly expressed in protoplasts isolated from peanut cultivar Zhonhua 12 by PEG, and red fluorescence was successfully detected. Then, DsRed2 was expressed in peanut plants Zhonghua 12 driven by 35S promoter via Agrobacterium tumefaciens-mediated transformation. Red fluorescence was visually observed in calli and regenerated shoots, as well as in roots, leaves, flowers, fresh pod shells and mature seeds, suggesting that transgenic screening could be initiated at the early stage of transformation, and continued to the progeny. Upon screening with DsRed2, the positive plant rate was increased from 56.9% to 100%. The transgenic line was then used as the male parent to be crossed with Zhonghua 24, and the hybrid seeds showed red fluorescence as well, indicating that DsRed2 could be applied to hybrid plant identification very efficiently. DsRed2 was also expressed in hairy roots of Huayu 23 via Agrobacterium rhizogenes-mediated transformation, and the transgenic roots were easily selected by red fluorescence. In summary, the DsRed2 is an ideal reporter to achieve maximum screening efficiency and accuracy in peanut genetic transformation.

Visualizing the Distribution of Lipids in Peanut Seeds by MALDI Mass Spectrometric Imaging

Peanut (also called groundnut, Arachis hypogaea L.) seeds are used for producing edible oils and functional foods, and offer a rich source of lipids, proteins and carbohydrates. However, the location of these metabolites has not yet been firmly established. In the present study, the matrix-assisted laser desorption/ionization mass spectrometric imaging (MALDI-MSI) technique was applied to investigate spatial distribution of lipids and other key components in seeds of three peanut cultivars (ZH9, KQBH, HP). A total of 103 metabolites, including 34 lipid compounds, were putatively identified by MALDI-MSI. The abundance and spatial distribution of glycerolipids (GLs) and glycerophospholipids (GPs) were compared among the three peanut cultivars. All the identified lysophosphatidylcholine (LPC), phosphatidylethanolamine (PE) and phosphatidylcholines (PCs) were distributed mainly in the inner part of seeds. The visualization of phosphatidic acids (PAs) and triacylglycerols (TGs) revealed a dramatic metabolic heterogeneity between the different tissues making up the seed. The non-homogeneous spatial distribution of metabolites appeared to be related to the different functions of particular tissue regions. These results indicated that MALDI-MSI could be useful for investigating the lipids of foodstuffs from a spatial perspective. The present study may contribute to the development of oil crops with higher oil yields, and to improvement of food processing.