Integrated transcriptomic and metabolomic analyses reveals anthocyanin biosynthesis in leaf coloration of quinoa (Chenopodium quinoa Willd.).

BACKGROUND: Quinoa leaves demonstrate a diverse array of colors, offering a potential enhancement to landscape aesthetics and the development of leisure-oriented sightseeing agriculture in semi-arid regions. This study utilized integrated transcriptomic and metabolomic analyses to investigate the mechanisms underlying anthocyanin synthesis in both emerald green and pink quinoa leaves. RESULTS: Integrated transcriptomic and metabolomic analyses indicated that both flavonoid biosynthesis pathway (ko00941) and anthocyanin biosynthesis pathway (ko00942) were significantly associated with anthocyanin biosynthesis. Differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were analyzed between the two germplasms during different developmental periods. Ten DEGs were verified using qRT-PCR, and the results were consistent with those of the transcriptomic sequencing. The elevated expression of phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), 4-coumarate CoA ligase (4CL) and Hydroxycinnamoyltransferase (HCT), as well as the reduced expression of flavanone 3-hydroxylase (F3H) and Flavonol synthase (FLS), likely cause pink leaf formation. In addition, bHLH14, WRKY46, and TGA indirectly affected the activities of CHS and 4CL, collectively regulating the levels of cyanidin 3-O-(3'', 6''-O-dimalonyl) glucoside and naringenin. The diminished expression of PAL, 4CL, and HCT decreased the formation of cyanidin-3-O-(6"-O-malonyl-2"-O-glucuronyl) glucoside, leading to the emergence of emerald green leaves. Moreover, the lowered expression of TGA and WRKY46 indirectly regulated 4CL activity, serving as another important factor in maintaining the emerald green hue in leaves N1, N2, and N3. CONCLUSION: These findings establish a foundation for elucidating the molecular regulatory mechanisms governing anthocyanin biosynthesis in quinoa leaves, and also provide some theoretical basis for the development of leisure and sightseeing agriculture.

Transcriptomics-metabolomics joint analysis: New highlight into the triterpenoid saponin biosynthesis in quinoa (Chenopodium quinoa Willd.).

Quinoa (Chenopodium quinoa Willd.) contains various physiologically active substances, including vitamins, polyphenols, flavonoids, phytosterols, and saponins. Research showed that saponins were the protective substances in the outer layer of quinoa seeds to defend against microbes, herbivores, and insects. Because the aglycones of quinoa saponins are triterpenoids, they are called triterpenoid saponins (TSs). In addition, the presence of TS imparted bitterness in quinoa and resulted in anticancer and anti-inflammatory effects. In this study, the seeds of low-saponin quinoa, NT376-2 (N), and high-saponin quinoa, B-12071(B), at 30 and 60 days after flowering (DAF) were used to measure the TS content and evaluated for their transcriptomic and metabolomic profiles. The amounts of TS were found to significantly differ between all possible comparisons: N and B at 30 DAF (N1_vs_B1), N and B at 60 DAF (N2_vs_B2), N at 30 DAF and 60 DAF (N1_vs_N2), and B at 30 DAF and 60 DAF (B1_vs_B2). RNA sequencing (RNA-seq) was used to screen differentially expressed genes (DEGs) and revealed 14,703 upregulated DEGs and 26,267 downregulated DEGs in the four comparison groups. The 311 overlapping DEGs found in the four comparisons were used for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to screen for DEGs related to TS biosynthesis in quinoa. Metabolomics analysis identified acetyl-CoA, 1-hydroxy-2-methyl-2-butenyl-4-diphosphate, farnesal, and (S)-2,3-epoxysqualene as the key differentially accumulated metabolites (DAMs). Transcriptomics-metabolomics joint analysis showed that triterpenoid biosynthesis and terpenoid backbone biosynthesis were the enriched pathways of TS biosynthesis; farnesal were the key DAMs shared in the four comparison groups and associated with 10 key candidate DEGs related to TS biosynthesis in quinoa. These results provided important references for in-depth research on the metabolic mechanism of TS in quinoa.

Corrigendum: Multi-Omics and miRNA Interaction Joint Analysis Highlight New Insights Into Anthocyanin Biosynthesis in Peanuts (Arachis hypogaea L.).

[This corrects the article DOI: 10.3389/fpls.2022.818345.].

Multi-Omics and miRNA Interaction Joint Analysis Highlight New Insights Into Anthocyanin Biosynthesis in Peanuts (Arachis hypogaea L.).

Peanut (Arachis hypogaea L.) is one of the most important economic and oil crops in the world. At present, peanut varieties with rich anthocyanin in testa are rare in the market, but the selection and breeding of varieties with the related traits has always attracted the attention of breeders. In this study, two peanut varieties with the pink and purple testa, G110 (G) and Z18-40 (Z) were used to conduct interaction joint analysis of multi-omics and miRNA-target gene. The anthocyanin content of Z18-40 was 7.49-8.62-folds higher than G110 on 30 DAF (days after flowering) and 45 DAF via Ultraviolet-visible Spectrophotometer (UV-5800, Shanghai, China). And then, a total of 14 candidate genes related with the anthocyanin biosynthesis were identified for correlation in different comparison groups (R 2 ≥ 0.80), among of a novel gene Ah21440 related with hydroxycinnamoyl transferase (HCT) biosynthesis was identified. In addition, Cyanidin 3-O-glucoside (Kuromanin, pmb0550) was the only common differentially accumulated metabolite (DAM) identified using multi-omics joint analysis in G1_vs_G2, Z1_vs_Z2, G1_vs_Z1, and G2_vs_Z2, respectively. Correlation analysis of miRNA-target genes and DEGs in the transcriptome shows that, AhmiR2950, AhmiR398, AhmiR50, and AhmiR51 regulated to HCT and chalcone biosynthesis related candidate genes (Ah21440, AhCHS, AhCHI). Lastly, all of 14 candidate genes and 4 differentially expressed miRNAs were validated using quantitative real-time PCR (qRT-PCR), which trends were consistent with that of the former transcriptome data. The results provide important reference for in-depth research on the anthocyanin metabolism mechanism in peanut testa.

Transcriptome and Co-expression Network Analyses Reveal Differential Gene Expression and Pathways in Response to Severe Drought Stress in Peanut (Arachis hypogaea L.).

Drought is one of the major abiotic stress factors limiting peanut production. It causes the loss of pod yield during the pod formation stage. Here, one previously identified drought-tolerant cultivar, "L422" of peanut, was stressed by drought (35 ± 5%) at pod formation stage for 5, 7, and 9 days. To analyze the drought effects on peanut, we conducted physiological and transcriptome analysis in leaves under well-watered (CK1, CK2, and CK3) and drought-stress conditions (T1, T2, and T3). By transcriptome analysis, 3,586, 6,730, and 8,054 differentially expressed genes (DEGs) were identified in "L422" at 5 days (CK1 vs T1), 7 days (CK2 vs T2), and 9 days (CK3 vs T3) of drought stress, respectively, and 2,846 genes were common DEGs among the three-time points. Furthermore, the result of weighted gene co-expression network analysis (WGCNA) revealed one significant module that was closely correlated between drought stress and physiological data. A total of 1,313 significantly up-/down-regulated genes, including 61 transcription factors, were identified in the module at three-time points throughout the drought stress stage. Additionally, six vital metabolic pathways, namely, "MAPK signaling pathway-plant," "flavonoid biosynthesis," "starch and sucrose metabolism," "phenylpropanoid biosynthesis," "glutathione metabolism," and "plant hormone signal transduction" were enriched in "L422" under severe drought stress. Nine genes responding to drought tolerance were selected for quantitative real-time PCR (qRT-PCR) verification and the results agreed with transcriptional profile data, which reveals the reliability and accuracy of transcriptome data. Taken together, these findings could lead to a better understanding of drought tolerance and facilitate the breeding of drought-resistant peanut cultivars.

Transcriptomic and metabolomic joint analysis reveals distinct flavonoid biosynthesis regulation for variegated testa color development in peanut (Arachis hypogaea L.).

Peanut is one of the important oil and economic crops, among which the variegated testa peanut is a unique member. The molecular mechanisms underlying the pigment synthesis in variegated testa are still unclear. Differentially expressed genes (DEGs) in the flavonoid metabolism pathway in pigmented areas indicated that there were 27 DEGs highly related to the synthesis of variegated testa color among 1,050 DEGs. Of these 27, 13 were up-regulated and 14 were down-regulated, including 3 PALs, 1 C4H, 2 CHSs, 1 F3H, 1 F3'H, 2 DFRs, 2 LARs, 2 IAAs, 4 bHLHs, and 9 MYBs. GO (Gene Ontology) analysis indicated that DEGs were similarly enriched in three branches. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis suggested flavonoid biosynthesis is the most direct metabolic pathway for the synthesis of testa variegation. The liquid chromatography-tandem mass spectrometry (LC-MS/MS) results showed that cyanidin and delphinidin were the primary metabolites that caused the color differences between the pigmented and the non-pigmented areas. Through the verification of 20 DEGs via qPCR, the results were consistent with transcriptome sequencing in four comparison groups. The results in this study lay the foundation for revealing the molecular regulation mechanisms of flavonoid synthesis in variegated testa peanut.

Analysis of Microbial Diversity and Community Structure of Peanut Pod and Its Surrounding Soil in Peanut Rot Epidemic Area.

Fungal communities are associated with healthy peanut crops and good crop production, through the regulation of pod rot disease. Rotted peanut pods and their surrounding soil samples were collected from locations in northern China. Fungal species were identified by next-generation sequencing, using the conserved sequences of their internal transcribed spacer regions. Results showed that rotted pod samples were rich in the phyla Ascomycota and Basidiomycota, and soil samples also contained these, plus Chytridiomycota and Zygomycota. There were regional variations in the species of fungi related to peanut pod rot and its surrounding soil, between locations. Fungal species of Cryptococcus and Fusarium were less abundant in soil samples than in rotted pod samples, and were the main pathogenic fungi identified in our study. Soil total carbon, nitrogen, and potassium had a strong influence on the fungal community, and total phosphorous and calcium ions, together with soil pH, had a modest influence. Only Mycosphaerella and Gibberella were not significantly affected by these factors. These findings may be of some help to control pod rot disease and reduce the production loss of peanut crops.

Molecular cloning and expression characterization of flavonol synthase genes in peanut (Arachis hypogaea).

Flavonol is an important functional bioactive substance in peanut seeds, and plays important roles responding to abiotic stress. The flavonol content is closely related to the activity and regulation of gene expression patterns of flavonol synthase (FLS). In this study, eight FLS genes, AhFLSs were cloned and their expression characterization in different peanut organ and seedling under different abiotic stress were conducted. The results showed that the expressions levels of AhFLSs were differed in all assayed peanut organs and seedlings under abiotic stress treatments. Expression levels of AhFLS2, AhFLS3, AhFLS4, and AhFLS6 were higher than those of other AhFLSs. The flavonol contents of peanut organs and seedlings under different abiotic stress were also determined using high performance liquid chromatography (HPLC). Dried mature peanut seeds were the organ tissue with the highest flavonol content, and flavonol content increased with seed development. Under abiotic stress treatments, the types of flavonols induced differed among stress treatments. Correlation analysis results suggested that eight AhFLS genes may have different functions in peanut. Moreover, changes in the expression of the eight genes appear to has substrate preference. These results can lay the foundation for the study of improving nutritional value of peanut seed and resistance of peanut plant.

Identification and expression analysis of WRKY gene family under drought stress in peanut (Arachis hypogaea L.).

WRKY transcription factors play crucial roles in regulation mechanism leading to the adaption of plants to the complex environment. In this study, AhWRKY family was comprehensively analyzed using bioinformatic approaches in combination with transcriptome sequencing data of the drought-tolerant peanut variety 'L422'. A total of 158 AhWRKY genes were identified and named according to their distribution on the chromosomes. Based on the structural features and phylogenetic analysis of AhWRKY proteins, the AhWRKY family members were classified into three (3) groups, of which group II included five (5) subgroups. Results of structure and conserved motifs analysis for the AhWRKY genes confirmed the accuracy of the clustering analysis. In addition, 12 tandem and 136 segmental duplication genes were identified. The results indicated that segmental duplication events were the main driving force in the evolution of AhWRKY family. Collinearity analysis found that 32 gene pairs existed between Arachis hypogaea and two diploid wild ancestors (Arachis duranensis and Arachis ipaensis), which provided valuable clues for phylogenetic characteristics of AhWRKY family. Furthermore, 19 stress-related cis-acting elements were found in the promoter regions. During the study of gene expression level of AhWRKY family members in response to drought stress, 73 differentially expressed AhWRKY genes were obtained to have been influenced by drought stress. These results provide fundamental insights for further study of WRKY genes in peanut drought resistance.

Construction of High-Density Genetic Map and Mapping Quantitative Trait Loci for Growth Habit-Related Traits of Peanut (Arachis hypogaea L.).

Plant growth habit is an important and complex agronomic trait and is associated with yield, disease resistance, and mechanized harvesting in peanuts. There are at least two distinct growth habits (erect and prostrate) and several intermediate forms existing in the peanut germplasm. A recombinant inbred line population containing 188 individuals was developed from a cross of "Jihua 5" and "M130" for genetically dissecting the architecture of the growth habit. A new high-density genetic linkage map was constructed by using specific locus amplified fragment sequencing technology. The map contains 2,808 single-nucleotide polymorphism markers distributed on 20 linkage groups with a total length of 1,308.20 cM and an average inter-marker distance of 0.47 cM. The quantitative trait locus (QTL) analysis of the growth habit-related traits was conducted based on phenotyping data from seven environments. A total of 39 QTLs for growth habit-related traits was detected on 10 chromosomes explaining 4.55-27.74% of the phenotypic variance, in which 6 QTLs were for lateral branch angle, 8 QTLs were for extent radius, 7 QTLs were for the index of plant type, 11 QTLs were for main stem height, and 7 QTLs were for lateral branch length. Among these QTLs, 12 were co-localized on chromosome B05 spanning an approximately 0.17 Mb physical interval in comparison with the allotetraploid reference genome of "Tifrunner." Analysis of the co-localized genome region has shown that the putative genes are involved in light and hormones and will facilitate peanut growth habit molecular breeding and study of peanut domestication.

Selection of suitable reference genes for abiotic stress-responsive gene expression studies in peanut by real-time quantitative PCR

Background: Because of its strong specificity and high accuracy, real-time quantitative PCR (RT-qPCR) has been a widely used method to study the expression of genes responsive to stress. It is crucial to have a suitable set of reference genes to normalize target gene expression in peanut under different conditions using RT-qPCR. In this study, 11 candidate reference genes were selected and examined under abiotic stresses (drought, salt, heavy metal, and low temperature) and hormone (SA and ABA) conditions as well as across different organ types. Three statistical algorithms (geNorm, NormFinder and BestKeeper) were used to evaluate the expression stabilities of reference genes, and the comprehensive rankings of gene stability were generated. Results: The results indicated that ELF1B and YLS8 were the most stable reference genes under PEG-simulated drought treatment. For high-salt treatment using NaCl, YLS8 and GAPDH were the most stable genes. Under CdCl2 treatment, UBI1 and YLS8 were suitable as stable reference genes. UBI1, ADH3, and ACTIN11 were sufficient for gene expression normalization in low-temperature experiment. All the 11 candidate reference genes showed relatively high stability under hormone treatments. For organs subset, UBI1, GAPDH, and ELF1B showed the maximum stability. UBI1 and ADH3 were the top two genes that could be used reliably in all the stress conditions assessed. Furthermore, the necessity of the reference genes screened was further confirmed by the expression pattern of AnnAhs. Conclusions: The results perfect the selection of stable reference genes for future gene expression studies in peanut and provide a list of reference genes that may be used in the future.

[Ecological adaptability evaluation of peanut cultivars based on biomass and nutrient accumulation].

To identify the good peanut cultivars with the properties of high yield, high nutrient use efficiency and wide adaptability, 19 selected peanut cultivars were planted in the low champaign area and piedmont plain area of Hebei Province. By using principal component analysis, the adaptability of these 19 cultivars was evaluated for different ecological regions through comparing their 16 main traits including biomass and nutrient parameters. According to the critical value of principal component (>1.0), the 16 biomass and nutrient characteristics were integrated into 4 principal components which accounted for 85% of the original information. The results indicated that there were obvious differences in yield and nutrient use efficiency for the peanut cultivars in different ecological regions. The 19 peanut cultivars were classified into 2 groups according to their ecological adaptability, and the cultivars from the group with wide adaptability could further be divided into 3 categories according to their yield and nutrient use efficiency. Among these cultivars, Yuhua 9719, Jihua 0212-4, Weihua 10, Yuhua 15, Puhua 28 and Jihua 10 were selected as the better peanut cultivars with the properties of high yield, high nutrient use efficiency and wide adaptability.

Molecular cloning and characterization of annexin genes in peanut (Arachis hypogaea L.).

Annexin, Ca(2+) or phospholipid binding proteins, with many family members are distributed throughout all tissues during plant growth and development. Annexins participate in a number of physiological processes, such as exocytosis, cell elongation, nodule formation in legumes, maturation and stress response. Six different full-length cDNAs and two partial-length cDNAs of peanut, (AnnAh1, AnnAh2, AnnAh3, AnnAh5, AnnAh6, AnnAh7, AnnAh4 and AnnAh8) encoding annexin proteins, were isolated and characterized using a RT-PCR/RACE-PCR based strategy. The predicted molecular masses of these annexins were 36.0kDa with acidic pIs of 5.97-8.81. ANNAh1, ANNAh2, ANNAh3, ANNAh5, ANNAh6 and ANNAh7 shared sequence similarity from 35.76 to 66.35% at amino acid level. Phylogenetic analysis revealed their evolutionary relationships with corresponding orthologous sequences in soybean and deduced proteins in various plant species. Real-time quantitative assays indicated that these genes were differentially expressed in various organs. Transcript level analysis for six annexin genes under stress conditions showed that these genes were regulated by drought, salinity, heavy metal stress, low temperature and hormone. Additionally, the prediction of cis-regulatory element suggested that different cis-responsive elements including stress- and hormone-responsive-related elements could respond to various stress conditions. These results indicated that members of AnnAhs family may play important roles in the adaptation of peanut to various environmental stresses.

Resistance switching in electrodeposited magnetite superlattices.

Defect-chemistry magnetite superlattices and compositional superlattices in the magnetite/zinc ferrite system are electrodeposited as epitaxial films onto single-crystal Au(111). The defect-chemistry superlattices have alternating nanolayers with different Fe(III)/Fe(II) ratios, whereas the compositional superlattices have alternating nanolayers with different Zn/Fe ratios. The electrochemical/chemical (EC) nature of the electrodeposition reaction is exploited to deposit the superlattices by pulsing the applied potential during deposition. The defect-chemistry superlattices show low-to-high and high-to-low resistance switching that may be applicable to the fabrication of resistive random access memory (RRAM).