Comprehensive Transcriptome Analyses Reveal Candidate Genes for Variation in Seed Size/Weight During Peanut (Arachis hypogaea L.) Domestication

Omics-driven advances in the understanding of regulatory landscape of peanut seed development

Peanuts (Arachis hypogaea) are an essential oilseed crop known for their unique developmental process, characterized by aerial flowering followed by subterranean fruit development. This crop is polyploid, consisting of A and B subgenomes, which complicates its genetic analysis. The advent and progression of omics technologies-encompassing genomics, transcriptomics, proteomics, epigenomics, and metabolomics-have significantly advanced our understanding of peanut biology, particularly in the context of seed development and the regulation of seed-associated traits. Following the completion of the peanut reference genome, research has utilized omics data to elucidate the quantitative trait loci (QTL) associated with seed weight, oil content, protein content, fatty acid composition, sucrose content, and seed coat color as well as the regulatory mechanisms governing seed development. This review aims to summarize the advancements in peanut seed development regulation and trait analysis based on reference genome-guided omics studies. It provides an overview of the significant progress made in understanding the molecular basis of peanut seed development, offering insights into the complex genetic and epigenetic mechanisms that influence key agronomic traits. These studies highlight the significance of omics data in profoundly elucidating the regulatory mechanisms of peanut seed development. Furthermore, they lay a foundational basis for future research on trait-related functional genes, highlighting the pivotal role of comprehensive genomic analysis in advancing our understanding of plant biology.

Low-concentration exogenous 3-indoleacetic acid improves fruit-setting rate of Marsdenia tenacissima by inhibiting the expression of embryo abortion-related genes

Marsdenia tenacissima is a medicinal plant characterized by many flowers, few fruits, and a low fruit-setting rate. Exogenous auxins can improve the fruit-setting rate of plants; however, their impacts on M. tenacissima and regulatory mechanisms remain unclear. In this study, we conducted a field experiment to determine the fruit-setting rate, seed-setting rate, fruit size, and changes in transcriptional expression of related genes by spraying 10 and 50 mg·L-1 of 3-indoleacetic acid (IAA). The control plants were sprayed with distilled water. Our results indicated that the fruit-setting rate was 0.15 when treated with 10 mg·L-1 of IAA, which was 2.76-fold higher than that of the control. Compared with that of the control, the number of differentially expressed genes (DEGs) regulated by 10 mg·L-1 of IAA was 28.6-fold higher than that regulated by 50 mg·L-1 of IAA. These DEGs were closely related to hormone metabolism and fruit development. By transcriptome analysis, spraying 10 mg·L-1 of IAA increased the expressions of STP6, MYB17, and LAX3 and reduced those of CXE18, ILR1-like 3, and SAUR50; this possibly affected the ovule, embryo, and fruit development, thereby elevating the fruit-setting rate of M. tenacissima. Our results indicated that low IAA concentration increased the fruit-setting rate of M. tenacissima, providing theoretical and practical support for promoting the seed yield of M. tenacissima.

Global Transcriptome and Co-Expression Network Analyses Revealed Hub Genes Controlling Seed Size/Weight and/or Oil Content in Peanut

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin-proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut.

Identification of QTL for kernel weight and size and analysis of the pentatricopeptide repeat (PPR) gene family in cultivated peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is an important oilseed crop worldwide. Improving its yield is crucial for sustainable peanut production to meet increasing food and industrial requirements. Deciphering the genetic control underlying peanut kernel weight and size, which are essential components of peanut yield, would facilitate high-yield breeding. A high-density single nucleotide polymorphism (SNP)-based linkage map was constructed using a recombinant inbred lines (RIL) population derived from a cross between the variety Yuanza9102 and a germplasm accession wt09-0023. Kernel weight and size quantitative trait loci (QTLs) were co-localized to a 0.16 Mb interval on Arahy07 using inclusive composite interval mapping (ICIM). Analysis of SNP, and Insertion or Deletion (INDEL) markers in the QTL interval revealed a gene encoding a pentatricopeptide repeat (PPR) superfamily protein as a candidate closely linked with kernel weight and size in cultivated peanut. Examination of the PPR gene family indicated a high degree of collinearity of PPR genes between A. hypogaea and its diploid progenitors, Arachis duranensis and Arachis ipaensis. The candidate PPR gene, Arahy.JX1V6X, displayed a constitutive expression pattern in developing seeds. These findings lay a foundation for further fine mapping of QTLs related to kernel weight and size, as well as validation of candidate genes in cultivated peanut.

Dynamic DNA methylation modification in peanut seed development

Cytosine methylation is an important epigenetic modification involved in regulation of plant development. However, the epigenetic mechanisms governing peanut seed development remain unclear. Herein, we generated DNA methylation profiles of developmental seeds of peanut H2014 and its smaller seed mutant H1314 at 15 and 60 days after pegging (DAP, S1, S4). Accompanying seed development, globally elevated methylation was observed in both lines. The mutant had a higher methylation level of 31.1% than wild type at S4, and 27.1-35.9% of the differentially methylated regions (DMRs) between the two lines were distributed in promoter or genic regions at both stages. Integrated methylome and transcriptome analysis revealed important methylation variations closely associated with seed development. Furthermore, some genes showed significantly negative correlation of expression with the methylation level within promoter or gene body. The results provide insights into the roles of DNA methylation in peanut seed development.

Multi-Omics Profiling Identifies Candidate Genes Controlling Seed Size in Peanut

Seed size is the major yield component and a key target trait that is selected during peanut breeding. However, the mechanisms that regulate peanut seed size are unknown. Two peanut mutants with bigger seed size were isolated in this study by ⁶⁰Co treatment of a common peanut landrace, Huayu 22, and were designated as the "big seed" mutant lines (hybs). The length and weight of the seed in hybs were about 118% and 170% of those in wild-type (WT), respectively. We adopted a multi-omics approach to identify the genomic locus underlying the hybs mutants. We performed whole genome sequencing (WGS) of WT and hybs mutants and identified thousands of large-effect variants (SNPs and indels) that occurred in about four hundred genes in hybs mutants. Seeds from both WT and hybs lines were sampled 20 days after flowering (DAF) and were used for RNA-Seq analysis; the results revealed about one thousand highly differentially expressed genes (DEGs) in hybs compared to WT. Using a method that combined large-effect variants with DEGs, we identified 45 potential candidate genes that shared gene product mutations and expression level changes in hybs compared to WT. Among the genes, two candidate genes encoding cytochrome P450 superfamily protein and NAC transcription factors may be associated with the increased seed size in hybs. The present findings provide new information on the identification and functional research into candidate genes responsible for the seed size phenotype in peanut.

The Key to the Future Lies in the Past: Insights from Grain Legume Domestication and Improvement Should Inform Future Breeding Strategies

Crop domestication is a co-evolutionary process that has rendered plants and animals significantly dependent on human interventions for survival and propagation. Grain legumes have played an important role in the development of Neolithic agriculture some 12,000 years ago. Despite being early companions of cereals in the origin and evolution of agriculture, the understanding of grain legume domestication has lagged behind that of cereals. Adapting plants for human use has resulted in distinct morpho-physiological changes between the wild ancestors and domesticates, and this distinction has been the focus of several studies aimed at understanding the domestication process and the genetic diversity bottlenecks created. Growing evidence from research on archeological remains, combined with genetic analysis and the geographical distribution of wild forms, has improved the resolution of the process of domestication, diversification and crop improvement. In this review, we summarize the significance of legume wild relatives as reservoirs of novel genetic variation for crop breeding programs. We describe key legume features, which evolved in response to anthropogenic activities. Here, we highlight how whole genome sequencing and incorporation of omics-level data have expanded our capacity to monitor the genetic changes accompanying these processes. Finally, we present our perspective on alternative routes centered on de novo domestication and re-domestication to impart significant agronomic advances of novel crops over existing commodities. A finely resolved domestication history of grain legumes will uncover future breeding targets to develop modern cultivars enriched with alleles that improve yield, quality and stress tolerance.

Transcriptome Analysis and Gene Expression Profiling of the Peanut Small Seed Mutant Identified Genes Involved in Seed Size Control

Seed size is a key factor affecting crop yield and a major agronomic trait concerned in peanut (Arachis hypogaea L.) breeding. However, little is known about the regulation mechanism of peanut seed size. In the present study, a peanut small seed mutant1 (ssm1) was identified through irradiating peanut cultivar Luhua11 (LH11) using ⁶⁰Coγ ray. Since the globular embryo stage, the embryo size of ssm1 was significantly smaller than that of LH11. The dry seed weight of ssm1 was only 39.69% of the wild type LH14. The seeds were wrinkled with darker seed coat. The oil content of ssm1 seeds were also decreased significantly. Seeds of ssm1 and LH11 were sampled 10, 20, and 40 days after pegging (DAP) and were used for RNA-seq. The results revealed that genes involved in plant hormones and several transcription factors related to seed development were differentially expressed at all three stages, especially at DAP10 and DAP20. Genes of fatty acid biosynthesis and late embryogenesis abundant protein were significantly decreased to compare with LH11. Interestingly, the gene profiling data suggested that PKp2 and/or LEC1 could be the key candidate genes leading to the small seed phenotype of the mutant. Our results provide valuable clues for further understanding the mechanisms underlying seed size control in peanut.

T. MP, Das S, Pawar PAM, Kumar A, Tripathi K, Kumar RR, Yadava DK, Kumar S, Dikshit HK. Comparative transcriptome analysis, unfolding the pathways regulating the seed-size trait in cultivated lentil (Lens culinaris Medik.)

Market class, cooking time, quality, and milled grain yield are largely influenced by the seed size and shape of the lentil (Lens culinaris Medik.); thus, they are considered to be important quality traits. To unfold the pathways regulating seed size in lentils, a transcriptomic approach was performed using large-seeded (L4602) and small-seeded (L830) genotypes. The study has generated nearly 375 million high-quality reads, of which 98.70% were properly aligned to the reference genome. Among biological replicates, very high similarity in fragments per kilobase of exon per million mapped fragments values (R > 0.9) showed the consistency of RNA-seq results. Various differentially expressed genes associated mainly with the hormone signaling and cell division pathways, transcription factors, kinases, etc. were identified as having a role in cell expansion and seed growth. A total of 106,996 unigenes were used for differential expression (DE) analysis. String analysis identified various modules having certain key proteins like Ser/Thr protein kinase, seed storage protein, DNA-binding protein, microtubule-associated protein, etc. In addition, some growth and cell division–related micro-RNAs like miR3457 (cell wall formation), miR1440 (cell proliferation and cell cycles), and miR1533 (biosynthesis of plant hormones) were identified as having a role in seed size determination. Using RNA-seq data, 5254 EST-SSR primers were generated as a source for future studies aiming for the identification of linked markers. In silico validation using Genevestigator^® was done for the Ser/Thr protein kinase, ethylene response factor, and Myb transcription factor genes. It is of interest that the xyloglucan endotransglucosylase gene was found differentially regulated, suggesting their role during seed development; however, at maturity, no significant differences were recorded for various cell wall parameters including cellulose, lignin, and xylose content. This is the first report on lentils that has unfolded the key seed size regulating pathways and unveiled a theoretical way for the development of lentil genotypes having customized seed sizes.

The fate of drought-related genes after polyploidization in Arachis hypogaea cv. Tifrunner

Drought stress affects plant growth and development. Cultivated peanut (Arachis hypogaea) was formed by a cross between A. duranensis and A. ipaensis. The drought tolerance of A. duranensis and A. ipaensis is reportedly stronger than that of cultivated peanut. However, there has been little study of drought tolerance genes in Arachis. In this study, we compared drought tolerance genes between A. hypogaea cv. Tifrunner and its diploid donors. We have observed that polyploidization does not generate more drought tolerance genes in A. hypogaea cv. Tifrunner but promotes the loss of many ancient drought tolerance genes. Although putative drought tolerance genes occurred on gene duplication events in A. hypogaea cv. Tifrunner, most copies lacked drought tolerance. These findings suggest that the loss of drought tolerance genes in A. hypogaea cv. Tifrunner could possibly result in weaker drought tolerance. In addition, we have observed that the three Arachis species stochastically lost putative drought tolerance genes. The evolution of drought tolerance genes could possibly have correlated with environmental changes. Our results enhance the current understanding of drought tolerance and polyploidy evolution in Arachis species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01198-0.

Transcriptomic Analysis Reveals the Regulatory Networks and Hub Genes Controlling the Unsaturated Fatty Acid Contents of Developing Seed in Soybean

Soybean [Glycine max (L.) Merr.] is one of the most important crops, which produces about 25% of the world's edible oil. The nutritional value of soybean oil depends mostly on the relative contents of three unsaturated fatty acids (UFAs), i.e., oleic acid, linoleic acid (LA), and linolenic acid. However, the biosynthetic mechanism of UFAs remains largely unknown, and there are few studies on RNA-seq analysis of developing seeds. To identify the candidate genes and related pathways involved in the regulation of UFA contents during seed development in soybean, two soybean lines with different UFA profiles were selected from 314 cultivars and landraces originated from Southern China, and RNA-seq analysis was performed in soybean seeds at three developmental stages. Using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, a series of genes and pathways related to fatty acid metabolism were identified, and 40 days after flowering (DAF) was found to be the crucial period in the formation of UFA profiles. Further, weighted gene co-expression network analysis identified three modules with six genes whose functions were highly associated with the contents of oleic and LA. The detailed functional investigation of the networks and hub genes could further improve the understanding of the underlying molecular mechanism of UFA contents and might provide some ideas for the improvement in fatty acids profiles in soybean.

Genome-wide structural and functional features of single nucleotide polymorphisms revealed from the whole genome resequencing of 179 accessions of Arachis

Peanut being an important food, oilseed and fodder crop worldwide, its genetic improvement currently relies on genomics-assisted breeding (GAB). Since the level of marker polymorphism is limited in peanut, the availability of a large number of DNA markers is the prerequisite for GAB. Therefore, we detected 4,309,724 single nucleotide polymorphisms (SNPs) from the whole genome re-sequencing (WGRS) data of 178 peanut accessions along with the reference genome sequence of Tifrunner. SNPs were analyzed for their structural and functional features to conclude on their utility and employability in genetic and genomic studies. ISATGR278-18, a synthetic amphidiploid, showed the highest number of SNPs (2,505,266), while PI_628538 recorded the lowest number (19,058) of SNPs. A03 showed the highest number of SNPs, while B08 recorded the lowest number of SNPs. The number of accessions required to record 50% of the total SNPs varied from 11 to 13 across the chromosomes. The rate of transitions was more than that of transversions. Among the various chromosomal contexts, intergenic and intronic regions carried more SNPs than the exonic regions. SNP impact analysis indicated 2488 SNPs with high impact due to gain of stop codons, variations in splice acceptors and splice donors, and loss of start codons. Of the 4,309,723 SNPs, 46,087 had the highest polymorphic information content (PIC) of 0.375. As an illustration of application, the drought-tolerant accession C76-16 was compared with A72 (an accession with high-stress rating) to identify 637,833 SNPs, of which 418 had high impact substitutions. Overall, these structural and functional features of the SNPs will be of immense importance for their utility in genetic and genomic studies in peanut.

Phenotypic Characteristics and Transcriptome of Cucumber Male Flower Development Under Heat Stress

Cucumber (Cucumis sativus L.) is an important vegetable crop, which is thermophilic not heat resistant. High-temperature stress always results in sterility at reproductive stage. In the present study, we evaluate the male flower developmental changes under normal (CK) and heat stress (HS) condition. After HS, the activities of peroxidase (POD) and superoxide dismutase (SOD) and the contents of malondialdehyde (MDA) were increased. In addition, the pollen fertility was significantly decreased; and abnormal tapetum and microspore were observed by paraffin section. Transcriptome analysis results presented that total of 5828 differentially expressed genes (DEGs) were identified after HS. Among these DEGs, 20 DEGs were found at four stages, including DNA binding transcription factor, glycosyltransferase, and wound-responsive family protein. The gene ontology term of carbohydrate metabolic process was significantly enriched in all anther stages, and many saccharides and starch synthase-related genes, such as invertase, sucrose synthase, and starch branching enzyme, were significantly different expressed in HS compared with CK. Furthermore, co-expression network analysis showed a module (midnightblue) strongly consistent with HS, and two hub genes (CsaV3_6G004180 and CsaV3_5G034860) were found with a high degree of connectivity to other genes. Our results provide comprehensive understandings on male flower development in cucumber under HS.