Genome of an allotetraploid wild peanut Arachis monticola: a de novo assembly

Homologous mapping yielded a comprehensive predicted protein-protein interaction network for peanut (Arachis hypogaea L.)

BACKGROUND

Protein-protein interactions are the primary means through which proteins carry out their functions. These interactions thus have crucial roles in life activities. The wide availability of fully sequenced animal and plant genomes has facilitated establishment of relatively complete global protein interaction networks for some model species. The genomes of cultivated and wild peanut (Arachis hypogaea L.) have also been sequenced, but the functions of most of the encoded proteins remain unclear.

RESULTS

We here used homologous mapping of validated protein interaction data from model species to generate complete peanut protein interaction networks for A. hypogaea cv. 'Tifrunner' (282,619 pairs), A. hypogaea cv. 'Shitouqi' (256,441 pairs), A. monticola (440,470 pairs), A. duranensis (136,363 pairs), and A. ipaensis (172,813 pairs). A detailed analysis was conducted for a putative disease-resistance subnetwork in the Tifrunner network to identify candidate genes and validate functional interactions. The network suggested that DX2UEH and its interacting partners may participate in peanut resistance to bacterial wilt; this was preliminarily validated with overexpression experiments in peanut.

CONCLUSION

Our results provide valuable new information for future analyses of gene and protein functions and regulatory networks in peanut.

Chloroplast and whole-genome sequencing shed light on the evolutionary history and phenotypic diversification of peanuts

Cultivated peanut (Arachis hypogaea L.) is a widely grown oilseed crop worldwide; however, the events leading to its origin and diversification are not fully understood. Here by combining chloroplast and whole-genome sequence data from a large germplasm collection, we show that the two subspecies of A. hypogaea (hypogaea and fastigiata) likely arose from distinct allopolyploidization and domestication events. Peanut genetic clusters were then differentiated in relation to dissemination routes and breeding efforts. A combination of linkage mapping and genome-wide association studies allowed us to characterize genes and genomic regions related to main peanut morpho-agronomic traits, namely flowering pattern, inner tegument color, growth habit, pod/seed weight and oil content. Together, our findings shed light on the evolutionary history and phenotypic diversification of peanuts and might be of broad interest to plant breeders.

Identification and application of a candidate gene AhAftr1 for aflatoxin production resistance in peanut seed (Arachis hypogaea L.)

INTRODUCTION

Peanut is susceptible to infection of Aspergillus fungi and conducive to aflatoxin contamination, hence developing aflatoxin-resistant variety is highly meaningful. Identifying functional genes or loci conferring aflatoxin resistance and molecular diagnostic marker are crucial for peanut breeding.

OBJECTIVES

This work aims to (1) identify candidate gene for aflatoxin production resistance, (2) reveal the related resistance mechanism, and (3) develop diagnostic marker for resistance breeding program.

METHODS

Resistance to aflatoxin production in a recombined inbred line (RIL) population derived from a high-yielding variety Xuhua13 crossed with an aflatoxin-resistant genotype Zhonghua 6 was evaluated under artificial inoculation for three consecutive years. Both genetic linkage analysis and QTL-seq were conducted for QTL mapping. The candidate gene was further fine-mapped using a secondary segregation mapping population and validated by transgenic experiments. RNA-Seq analysis among resistant and susceptible RILs was used to reveal the resistance pathway for the candidate genes.

RESULTS

The major effect QTL qAFTRA07.1 for aflatoxin production resistance was mapped to a 1.98 Mbp interval. A gene, AhAftr1 (Arachis hypogaea Aflatoxin resistance 1), was detected structure variation (SV) in leucine rich repeat (LRR) domain of its production, and involved in disease resistance response through the effector-triggered immunity (ETI) pathway. Transgenic plants with overexpression of AhAftr1(ZH6) exhibited 57.3% aflatoxin reduction compared to that of AhAftr1(XH13). A molecular diagnostic marker AFTR.Del.A07 was developed based on the SV. Thirty-six lines, with aflatoxin content decrease by over 77.67% compared to the susceptible control Zhonghua12 (ZH12), were identified from a panel of peanut germplasm accessions and breeding lines through using AFTR.Del.A07.

CONCLUSION

Our findings would provide insights of aflatoxin production resistance mechanisms and laid meaningful foundation for further breeding programs.

Advances in the evolution research and genetic breeding of peanut

Peanut is an important cash crop used in oil, food and feed in our country. The rapid development of sequencing technology has promoted the research on the related aspects of peanut genetic breeding. This paper reviews the research progress of peanut origin and evolution, genetic breeding, molecular markers and their applications, genomics, QTL mapping and genome selection techniques. The main problems of molecular genetic breeding in peanut research worldwide include: the narrow genetic resources of cultivated species, unstable genetic transformation and unclear molecular mechanism of important agronomic traits. Considering the severe challenges regarding the supply of edible oil, and the main problems in peanut production, the urgent research directions of peanut are put forward: The de novo domestication and the exploitation of excellent genes from wild resources to improve modern cultivars; Integration of multi-omics data to enhance the importance of big data in peanut genetics and breeding; Cloning the important genes related to peanut agronomic traits and analyzing their fine regulation mechanisms; Precision molecular design breeding and using gene editing technology to accurately improve the key traits of peanut.

Chromosome-level genome assembly of Guide Black-Fur sheep (Ovis aries)

Guide Black-Fur sheep (GD) is a breed of Tibetan sheep (Ovis aries) that lives in the Qinghai-Tibetan plateau region at an altitude of over 4,000 m. However, a lack of genomic information has made it difficult to understand the high-altitude adaptation of these sheep. We sequenced and assembled the GD reference genome using PacBio, Hi-C, and Illumina sequencing technologies. The final assembled genome size was 2.73 Gb, with a contig N50 of 20.30 Mb and a scaffold N50 of 107.63 Mb. The genome is predicted to contain 20,759 protein-coding genes, of which 98.42 have functional annotations. Repeat elements account for approximately 52.2% of the genomic landscape. The completeness of the GD genome assembly is highlighted by a BUSCO score of 93.1%. This high-quality genome assembly provides a critical resource for future molecular breeding and genetic improvement of Tibetan sheep.

WRKY transcription factors modulate flowering time in four Arachis species: a bioinformatics analysis

BACKGROUND

WRKY proteins are important transcription factors (TFs) in plants, involved in growth and development and responses to environmental changes. Although WRKY TFs have been studied at the genome level in Arachis genus, including oil crop and turfgrass, their regulatory networks in controlling flowering time remain unclear. The aim of this study was to predict the molecular mechanisms of WRKY TFs regulation flowering time in Arachis genus at the genome level using bioinformatics approaches.

RESULTS

The flowering-time genes of Arachis genus were retrieved from the flowering-time gene database. The regulatory networks between WRKY TFs and downstream genes in Arachis genus were predicted using bioinformatics tools. The results showed that WRKY TFs were involved in aging, autonomous, circadian clock, hormone, photoperiod, sugar, temperature, and vernalization pathways to modulate flowering time in Arachis duranensis, Arachis ipaensis, Arachis monticola, and Arachis hypogaea cv. Tifrunner. The WRKY TF binding sites in homologous flowering-time genes exhibited asymmetric evolutionary pattern, indicating that the WRKY TFs interact with other transcription factors to modulate flowering time in the four Arachis species. Protein interaction network analysis showed that WRKY TFs interacted with FRUITFULL and APETALA2 to modulate flowering time in the four Arachis species. WRKY TFs implicated in regulating flowering time had low expression levels, whereas their interaction proteins had varying expression patterns in 22 tissues of A. hypogaea cv. Tifrunner. These results indicate that WRKY TFs exhibit antagonistic or synergistic interactions with the associated proteins.

CONCLUSIONS

This study reveals complex regulatory networks through which WRKY TFs modulate flowering time in the four Arachis species using bioinformatics approaches.

The complete chloroplast genome of Arachis lutescens Krapov. & Rigoni (Fabaceae)

Arachis lutescens Krapov. & Rigoni 1958 is an important species due to their potentially extensive applications for cultivated peanut breeding. The whole chloroplast genome of A. lutescens was successfully assembled and annotated for the first time. The complete chloroplast genome of A. lutescens is a typically circular structure of 156,398 bp with a GC content of 36.3%. It comprises a large single-copy (LSC) region of 85,950 bp, a small single-copy (SSC) region of 18,800 bp, and two inverted repeat regions (IRs) of 25,824 bp, each. The plastome of A. lutescens contains a total of 125 genes, including 81 protein-coding genes, 36 tRNAs, and eight rRNAs. The phylogenetic analysis strongly supports the close relationship between A. lutescens and cultivated peanut clades. This study contributes to our understanding of the molecular characteristics and evolutionary relationships of this plant species.

Designing future peanut: the power of genomics-assisted breeding

KEY MESSAGE

Integrating GAB methods with high-throughput phenotyping, genome editing, and speed breeding hold great potential in designing future smart peanut cultivars to meet market and food supply demands. Cultivated peanut (Arachis hypogaea L.), a legume crop greatly valued for its nourishing food, cooking oil, and fodder, is extensively grown worldwide. Despite decades of classical breeding efforts, the actual on-farm yield of peanut remains below its potential productivity due to the complicated interplay of genotype, environment, and management factors, as well as their intricate interactions. Integrating modern genomics tools into crop breeding is necessary to fast-track breeding efficiency and rapid progress. When combined with speed breeding methods, this integration can substantially accelerate the breeding process, leading to faster access of improved varieties to farmers. Availability of high-quality reference genomes for wild diploid progenitors and cultivated peanuts has accelerated the process of gene/quantitative locus discovery, developing markers and genotyping assays as well as a few molecular breeding products with improved resistance and oil quality. The use of new breeding tools, e.g., genomic selection, haplotype-based breeding, speed breeding, high-throughput phenotyping, and genome editing, is probable to boost genetic gains in peanut. Moreover, renewed attention to efficient selection and exploitation of targeted genetic resources is also needed to design high-quality and high-yielding peanut cultivars with main adaptation attributes. In this context, the combination of genomics-assisted breeding (GAB), genome editing, and speed breeding hold great potential in designing future improved peanut cultivars to meet market and food supply demands.

A genomic variation map provides insights into peanut diversity in China and associations with 28 agronomic traits

Peanut (Arachis hypogaea L.) is an important allotetraploid oil and food legume crop. China is one of the world's largest peanut producers and consumers. However, genomic variations underlying the migration and divergence of peanuts in China remain unclear. Here we reported a genome-wide variation map based on the resequencing of 390 peanut accessions, suggesting that peanuts might have been introduced into southern and northern China separately, forming two cultivation centers. Selective sweep analysis highlights asymmetric selection between the two subgenomes during peanut improvement. A classical pedigree from South China offers a context for the examination of the impact of artificial selection on peanut genome. Genome-wide association studies identified 22,309 significant associations with 28 agronomic traits, including candidate genes for plant architecture and oil biosynthesis. Our findings shed light on peanut migration and diversity in China and provide valuable genomic resources for peanut improvement.

Coordinated Lipid Mobilization during Seed Development and Germination in Peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is one of the most important oil crops in the world due to its lipid-rich seeds. Lipid accumulation and degradation play crucial roles in peanut seed maturation and seedling establishment, respectively. Here, we utilized lipidomics and transcriptomics to comprehensively identify lipids and the associated functional genes that are important in the development and germination processes of a large-seed peanut variety. A total of 332 lipids were identified; triacylglycerols (TAGs) and diacylglycerols were the most abundant during seed maturation, constituting 70.43 and 16.11%, respectively, of the total lipids. Significant alterations in lipid profiles were observed throughout seed maturation and germination. Notably, TAG (18:1/18:1/18:2) and (18:1/18:2/18:2) peaked at 23386.63 and 23392.43 nmol/g, respectively, at the final stage of seed development. Levels of hydroxylated TAGs (HO-TAGs) increased significantly during the initial stage of germination. Accumulation patterns revealed an inverse relationship between free fatty acids and TAGs. Lipid degradation was determined to be regulated by diacylglycerol acyltransferase, triacylglycerol lipase, and associated transcription factors, predominantly yielding oleic acid, linoleic acid, and linolenic acid. Collectively, the results of this study provide valuable insights into lipid dynamics during the development and germination of large-seed peanuts, gene resources, and guiding future research into lipid accumulation in an economically important crop.

WRKY transcription factors in Arachis hypogaea and its donors: From identification to function prediction

WRKY transcription factors (TFs) play important roles in plant growth and development and responses to abiotic and biotic stresses. Since the initial isolation of a WRKY TF in Ipomoea batatas in 1994, WRKY TFs have been identified in plants, protozoa, and fungi. Peanut (Arachis hypogaea) is a key oil and protein crop for humans and a forage source for animal consumption. Several Arachis genomes have been sequenced and genome-wide WRKY TFs have been identified. In this review, we summarized WRKY TFs and their functions in A. hypogaea and its donors. We also standardized the nomenclature for Arachis WRKY TFs to ensure uniformity. We determined the evolutionary relationships between Arachis and Arabidopsis thaliana WRKY (AtWRKY) TFs using a phylogenetic analysis. Biological functions and regulatory networks of Arachis WRKY TFs were predicted using AtWRKY TFs. Thus, this review paves the way for studies of Arachis WRKY TFs.

PSW1, an LRR receptor kinase, regulates pod size in peanut

Pod size is a key agronomic trait that greatly determines peanut yield, the regulatory genes and molecular mechanisms that controlling peanut pod size are still unclear. Here, we used quantitative trait locus analysis to identify a peanut pod size regulator, POD SIZE/WEIGHT1 (PSW1), and characterized the associated gene and protein. PSW1 encoded leucine-rich repeat receptor-like kinase (LRR-RLK) and positively regulated pod stemness. Mechanistically, this allele harbouring a 12-bp insertion in the promoter and a point mutation in the coding region of PSW1 causing a serine-to-isoleucine (S618I) substitution substantially increased mRNA abundance and the binding affinity of PSW1 for BRASSINOSTEROID INSENSITIVE1-ASSOCIATED RECEPTOR KINASE 1 (BAK1). Notably, PSW1HapII (super-large pod allele of PSW1) expression led to up-regulation of a positive regulator of pod stemness PLETHORA 1 (PLT1), thereby resulting in larger pod size. Moreover, overexpression of PSW1HapII increased seed/fruit size in multiple plant species. Our work thus discovers a conserved function of PSW1 that controls pod size and provides a valuable genetic resource for breeding high-yield crops.

Mapping Quantitative Trait Loci (QTLs) for Hundred-Pod and Hundred-Seed Weight under Seven Environments in a Recombinant Inbred Line Population of Cultivated Peanut (Arachis hypogaea L.)

The cultivated peanut (Arachis hypogaea L.) is a significant oil and cash crop globally. Hundred-pod and -seed weight are important components for peanut yield. To unravel the genetic basis of hundred-pod weight (HPW) and hundred-seed weight (HSW), in the current study, a recombinant inbred line (RIL) population with 188 individuals was developed from a cross between JH5 (JH5, large pod and seed weight) and M130 (small pod and seed weight), and was utilized to identify QTLs for HPW and HSW. An integrated genetic linkage map was constructed by using SSR, AhTE, SRAP, TRAP and SNP markers. This map consisted of 3130 genetic markers, which were assigned to 20 chromosomes, and covered 1998.95 cM with an average distance 0.64 cM. On this basis, 31 QTLs for HPW and HSW were located on seven chromosomes, with each QTL accounting for 3.7-10.8% of phenotypic variance explained (PVE). Among these, seven QTLs were detected under multiple environments, and two major QTLs were found on B04 and B08. Notably, a QTL hotspot on chromosome A08 contained seven QTLs over a 2.74 cM genetic interval with an 0.36 Mb physical map, including 18 candidate genes. Of these, Arahy.D52S1Z, Arahy.IBM9RL, Arahy.W18Y25, Arahy.CPLC2W and Arahy.14EF4H might play a role in modulating peanut pod and seed weight. These findings could facilitate further research into the genetic mechanisms influencing pod and seed weight in cultivated peanut.

Fine mapping of qAHPS07 and functional studies of AhRUVBL2 controlling pod size in peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.) is an important oil and cash crop. Pod size is one of the major traits determining yield and commodity characteristic of peanut. Fine mapping of quantitative trait locus (QTL) and identification of candidate genes associated with pod size are essential for genetic improvement and molecular breeding of peanut varieties. In this study, a major QTL related to pod size, qAHPS07, was fine mapped to a 36.46 kb interval on chromosome A07 using F2 , recombinant inbred line (RIL) and secondary F2 populations. qAHPS07 explained 38.6%, 23.35%, 37.48%, 25.94% of the phenotypic variation for single pod weight (SPW), pod length (PL), pod width (PW) and pod shell thickness (PST), respectively. Whole genome resequencing and gene expression analysis revealed that a RuvB-like 2 protein coding gene AhRUVBL2 was the most likely candidate for qAHPS07. Overexpression of AhRUVBL2 in Arabidopsis led to larger seeds and plants than the wild type. AhRUVBL2-silenced peanut seedlings represented small leaves and shorter main stems. Three haplotypes were identified according to three SNPs in the promoter of AhRUVBL2 among 119 peanut accessions. Among them, SPW, PW and PST of accessions carrying Hap_ATT represent 17.6%, 11.2% and 26.3% higher than those carrying Hap_GAC，respectively. In addition, a functional marker of AhRUVBL2 was developed. Taken together, our study identified a key functional gene of peanut pod size, which provides new insights into peanut pod size regulation mechanism and offers practicable markers for the genetic improvement of pod size-related traits in peanut breeding.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Evolution of NLR genes in genus Arachis reveals asymmetric expansion of NLRome in wild and domesticated tetraploid species

Arachis hypogaea is an allotetraploid crop widely grown in the world. Wild relatives of genus Arachis are the rich source of genetic diversity and high levels of resistance to combat pathogens and climate change. The accurate identification and characterization of plant resistance gene, nucleotide binding site leucine rich repeat receptor (NLRs) substantially contribute to the repertoire of resistances and improve production. In the current study, we have studied the evolution of NLR genes in genus Arachis and performed their comparative genomics among four diploids (A. duranensis, A. ipaensis, A. cardenasii, A. stenosperma) and two tetraploid (wild: A. monticola and domesticated: A. hypogaea) species. In total 521, 354, 284, 794, 654, 290 NLR genes were identified from A. cardenasii, A. stenosperma and A. duranensis, A. hypogaea, A. monticola and A. ipaensis respectively. Phylogenetic analysis and classification of NLRs revealed that they belong to 7 subgroups and specific subgroups have expanded in each genome leading towards divergent evolution. Gene gain and loss, duplication assay reveals that wild and domesticated tetraploids species have shown asymmetric expansion of NLRome in both sub-genome (AA and BB). A-subgenome of A. monticola exhibited significant contraction of NLRome while B-subgenome shows expansion and vice versa in case of A. hypogaea probably due to distinct natural and artificial selection pressure. In addition, diploid species A. cardenasii revealed the largest repertoire of NLR genes due to higher frequency of gene duplication and selection pressure. A. cardenasii and A. monticola can be regarded as putative resistance resources for peanut breeding program for introgression of novel resistance genes. Findings of this study also emphasize the application neo-diploids and polyploids due to higher quantitative expression of NLR genes. To the best of our knowledge, this is the first study that studied the effect of domestication and polyploidy on the evolution of NLR genes in genus Arachis to identify genomic resources for improving resistance of polyploid crop with global importance on economy and food security.

Fine mapping of a QTL and identification of candidate genes associated with cold tolerance during germination in peanut (Arachis hypogaea L.) on chromosome B09 using whole genome re-sequencing

Low temperatures significantly affect the growth and yield of peanuts. Temperatures lower than 12 °C are generally detrimental for the germination of peanuts. To date, there has been no report on precise information on the quantitative trait loci (QTL) for cold tolerance during the germination in peanuts. In this study, we developed a recombinant inbred line (RIL) population comprising 807 RILs by tolerant and sensitive parents. Phenotypic frequencies of germination rate low-temperature conditions among RIL population showed normally distributed in five environments. Then, we constructed a high density SNP-based genetic linkage map through whole genome re-sequencing (WGRS) technique and identified a major quantitative trait locus (QTL), qRGRB09, on chromosome B09. The cold tolerance-related QTLs were repeatedly detected in all five environments, and the genetic distance was 6.01 cM (46.74 cM - 61.75 cM) after taking a union set. To further confirm that qRGRB09 was located on chromosome B09, we developed Kompetitive Allele Specific PCR (KASP) markers for the corresponding QTL regions. A regional QTL mapping analysis, which was conducted after taking the intersection of QTL intervals of all environments into account, confirmed that qRGRB09 was between the KASP markers, G22096 and G220967 (chrB09:155637831-155854093), and this region was 216.26 kb in size, wherein a total of 15 annotated genes were detected. This study illustrates the relevance of WGRS-based genetic maps for QTL mapping and KASP genotyping that facilitated QTL fine mapping of peanuts. The results of our study also provided useful information on the genetic architecture underlying cold tolerance during germination in peanuts, which in turn may be useful for those engaged in molecular studies as well as crop improvement in the cold-stressed environment.

Investigation of NLR Genes Reveals Divergent Evolution on NLRome in Diploid and Polyploid Species in Genus Trifolium

Crop wild relatives contain a greater variety of phenotypic and genotypic diversity compared to their domesticated counterparts. Trifolium crop species have limited genetic diversity to cope with biotic and abiotic stresses due to artificial selection for consumer preferences. Here, we investigated the distribution and evolution of nucleotide-binding site leucine-rich repeat receptor (NLR) genes in the genus of Trifolium with the objective to identify reference NLR genes. We identified 412, 350, 306, 389 and 241 NLR genes were identified from Trifolium. subterraneum, T. pratense, T. occidentale, subgenome-A of T. repens and subgenome-B of T. repens, respectively. Phylogenetic and clustering analysis reveals seven sub-groups in genus Trifolium. Specific subgroups such as G4-CNL, CCG10-CNL and TIR-CNL show distinct duplication patterns in specific species, which suggests subgroup duplications that are the hallmarks of their divergent evolution. Furthermore, our results strongly suggest the overall expansion of NLR repertoire in T. subterraneum is due to gene duplication events and birth of gene families after speciation. Moreover, the NLRome of the allopolyploid species T. repens has evolved asymmetrically, with the subgenome -A showing expansion, while the subgenome-B underwent contraction. These findings provide crucial background data for comprehending NLR evolution in the Fabaceae family and offer a more comprehensive analysis of NLR genes as disease resistance genes.

Genetic mapping of AhVt1, a novel genetic locus that confers the variegated testa color in cultivated peanut (Arachis hypogaea L.) and its utilization for marker-assisted selection

INTRODUCTION

Peanut (Arachis hypogaea L.) is an important cash crop worldwide. Compared with the ordinary peanut with pure pink testa, peanut with variegated testa color has attractive appearance and a higher market value. In addition, the variegated testa represents a distinct regulation pattern of anthocyanin accumulation in integument cells.

METHODS

In order to identify the genetic locus underlying variegated testa color in peanut, two populations were constructed from the crosses between Fuhua 8 (pure-pink testa) and Wucai (red on white variegated testa), Quanhonghua 1 (pure-red testa) and Wucai, respectively. Genetic analysis and bulked sergeant analysis sequencing were applied to detect and identify the genetic locus for variegated testa color. Marker-assisted selection was used to develop new variegated testa peanut lines.

RESULTS

As a result, all the seeds harvested from the F1 individuals of both populations showed the variegated testa type with white trace. Genetic analysis revealed that the pigmentation of colored region in red on white variegated testa was controlled by a previous reported gene AhRt1, while the formation of white region (un-pigmented region) in variegated testa was controlled by another single genetic locus. This locus, named as AhVt1 (Arachis hypogaea Variegated Testa 1), was preliminary mapped on chromosome 08 through bulked sergeant analysis sequencing. Using a secondary mapping population derived from the cross between Fuhua 8 and Wucai, AhVt1 was further mapped to a 1.89-Mb genomic interval by linkage analysis, and several potential genes associated with the uneven distribution of anthocyanin, such as MADS-box, MYB, and Chalcone synthase-like protein, were harbored in the region. Moreover, the molecular markers closely linked to the AhVt1 were developed, and the new variegated testa peanut lines were obtained with the help of marker-assisted selection.

CONCLUSION

Our findings will accelerate the breeding program for developing new peanut varieties with "colorful" testa colors and laid a foundation for map-based cloning of gene responsible for variegated testa.

Dynamic Evolution of NLR Genes in Dalbergioids

Dalbergioid is a large group within the family Fabaceae that consists of diverse plant species distributed in distinct biogeographic realms. Here, we have performed a comprehensive study to understand the evolution of the nucleotide-binding leucine-rich repeats (NLRs) gene family in Dalbergioids. The evolution of gene families in this group is affected by a common whole genome duplication that occurred approximately 58 million years ago, followed by diploidization that often leads to contraction. Our study suggests that since diploidization, the NLRome of all groups of Dalbergioids is expanding in a clade-specific manner with fewer exceptions. Phylogenetic analysis and classification of NLRs revealed that they belong to seven subgroups. Specific subgroups have expanded in a species-specific manner, leading to divergent evolution. Among the Dalbergia clade, the expansion of NLRome in six species of the genus Dalbergia was observed, with the exception of Dalbergia odorifera, where a recent contraction of NLRome occurred. Similarly, members of the Pterocarpus clade genus Arachis revealed a large-scale expansion in the diploid species. In addition, the asymmetric expansion of NLRome was observed in wild and domesticated tetraploids after recent duplications in the genus Arachis. Our analysis strongly suggests that whole genome duplication followed by tandem duplication after divergence from a common ancestor of Dalbergioids is the major cause of NLRome expansion. To the best of our knowledge, this is the first ever study to provide insight toward the evolution of NLR genes in this important tribe. In addition, accurate identification and characterization of NLR genes is a substantial contribution to the repertoire of resistances among members of the Dalbergioids species.

Deciphering evolutionary dynamics of WRKY genes in Arachis species

BACKGROUND

Cultivated peanut (Arachis hypogaea), a progeny of the cross between A. duranensis and A. ipaensis, is an important oil and protein crop from South America. To date, at least six Arachis genomes have been sequenced. WRKY transcription factors (TFs) play crucial roles in plant growth, development, and response to abiotic and biotic stresses. WRKY TFs have been identified in A. duranensis, A. ipaensis, and A. hypogaea cv. Tifrunner; however, variations in their number and evolutionary patterns across various Arachis spp. remain unclear.

RESULTS

WRKY TFs were identified and compared across different Arachis species, including A. duranensis, A. ipaensis, A. monticola, A. hypogaea cultivars (cv.) Fuhuasheng, A. hypogaea cv. Shitouqi, and A. hypogaea cv. Tifrunner. The results showed that the WRKY TFs underwent dynamic equilibrium between diploid and tetraploid peanut species, characterized by the loss of old WRKY TFs and retention of the new ones. Notably, cultivated peanuts inherited more conserved WRKY orthologs from wild tetraploid peanuts than their wild diploid donors. Analysis of the W-box elements and protein-protein interactions revealed that different domestication processes affected WRKY evolution across cultivated peanut varieties. WRKY TFs of A. hypogaea cv. Fuhuasheng and Shitouqi exhibited a similar domestication process, while those of cv. Tifrunner of the same species underwent a different domestication process based on protein-protein interaction analysis.

CONCLUSIONS

This study provides new insights into the evolution of WRKY TFs in Arachis spp.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

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.

Advancing Grain Legumes Domestication and Evolution Studies with Genomics

Grain legumes were domesticated in parallel with cereals in several regions of the world and formed the economic basis of early farming cultures. Since then, legumes have played a vital role in human and animal diets and in fostering agrobiodiversity. Increasing grain legume cultivation will be crucial to safeguard nutritional security and the resilience of agricultural ecosystems across the globe. A better understanding of the molecular underpinnings of domestication and crop evolution of grain legumes may be translated into practical approaches in modern breeding programs to stabilize yield, which is threatened by evolving pathogens and changing climates. During recent decades, domestication research in all crops has greatly benefited from the fast progress in genomic technologies. Yet still, many questions surrounding the domestication and diversification of legumes remain unanswered. In this review, we assess the potential of genomic approaches in grain legume research. We describe the centers of origin and the crucial domestication traits of grain legumes. In addition, we survey the effect of domestication on both above-ground and below-ground traits that have economic importance. Finally, we discuss open questions in grain legume domestication and diversification and outline how to bridge the gap between the preservation of historic crop diversity and their utilization in modern plant breeding.

Comparison of ONT and CCS sequencing technologies on the polyploid genome of a medicinal plant showed that high error rate of ONT reads are not suitable for self-correction

BACKGROUND

Many medicinal plants are known for their complex genomes with high ploidy, heterozygosity, and repetitive content which pose severe challenges for genome sequencing of those species. Long reads from Oxford nanopore sequencing technology (ONT) or Pacific Biosciences Single Molecule, Real-Time (SMRT) sequencing offer great advantages in de novo genome assembly, especially for complex genomes with high heterozygosity and repetitive content. Currently, multiple allotetraploid species have sequenced their genomes by long-read sequencing. However, we found that a considerable proportion of these genomes (7.9% on average, maximum 23.7%) could not be covered by NGS (Next Generation Sequencing) reads (uncovered region by NGS reads, UCR) suggesting the questionable and low-quality of those area or genomic areas that can't be sequenced by NGS due to sequencing bias. The underlying causes of those UCR in the genome assembly and solutions to this problem have never been studied.

METHODS

In the study, we sequenced the tetraploid genome of Veratrum dahuricum (Turcz.) O. Loes (VDL), a Chinese medicinal plant, with ONT platform and assembled the genome with three strategies in parallel. We compared the qualities, coverage, and heterozygosity of the three ONT assemblies with another released assembly of the same individual using reads from PacBio circular consensus sequencing (CCS) technology, to explore the cause of the UCR.

RESULTS

By mapping the NGS reads against the three ONT assemblies and the CCS assembly, we found that the coverage of those ONT assemblies by NGS reads ranged from 49.15 to 76.31%, much smaller than that of the CCS assembly (99.53%). And alignment between ONT assemblies and CCS assembly showed that most UCR can be aligned with CCS assembly. So, we conclude that the UCRs in ONT assembly are low-quality sequences with a high error rate that can't be aligned with short reads, rather than genomic regions that can't be sequenced by NGS. Further comparison among the intermediate versions of ONT assemblies showed that the most probable origin of those errors is a combination of artificial errors introduced by "self-correction" and initial sequencing error in long reads. We also found that polishing the ONT assembly with CCS reads can correct those errors efficiently.

CONCLUSIONS

Through analyzing genome features and reads alignment, we have found the causes for the high proportion of UCR in ONT assembly of VDL are sequencing errors and additional errors introduced by self-correction. The high error rates of ONT-raw reads make them not suitable for self-correction prior to allotetraploid genome assembly, as the self-correction will introduce artificial errors to > 5% of the UCR sequences. We suggest high-precision CCS reads be used to polish the assembly to correct those errors effectively for polyploid genomes.

De novo genes in Arachis hypogaea cv. Tifrunner: systematic identification, molecular evolution, and potential contributions to cultivated peanut

De novo genes are derived from non-coding sequences, and they can play essential roles in organisms. Cultivated peanut (Arachis hypogaea) is a major oil and protein crop derived from a cross between Arachis duranensis and Arachis ipaensis. However, few de novo genes have been documented in Arachis. Here, we identified 381 de novo genes in A. hypogaea cv. Tifrunner based on comparison with five closely related Arachis species. There are distinct differences in gene expression patterns and gene structures between conserved and de novo genes. The identified de novo genes originated from ancestral sequence regions associated with metabolic and biosynthetic processes, and they were subsequently integrated into existing regulatory networks. De novo paralogs and homoeologs were identified in A. hypogaea cv. Tifrunner. De novo paralogs and homoeologs with conserved expression have mismatching cis-acting elements under normal growth conditions. De novo genes potentially have pluripotent functions in responses to biotic stresses as well as in growth and development based on quantitative trait locus data. This work provides a foundation for future research examining gene birth processes and gene function in Arachis and related taxa.

Genome-Wide Identification and Expression of FAR1 Gene Family Provide Insight Into Pod Development in Peanut (Arachis hypogaea)

The far-red-impaired response 1 (FAR1) transcription family were initially identified as important factors for phytochrome A (phyA)-mediated far-red light signaling in Arabidopsis; they play crucial roles in controlling the growth and development of plants. The reported reference genome sequences of Arachis, including A. duranensis, A. ipaensis, A. monticola, and A. hypogaea, and its related species Glycine max provide an opportunity to systematically perform a genome-wide identification of FAR1 homologous genes and investigate expression patterns of these members in peanut species. Here, a total of 650 FAR1 genes were identified from four Aarchis and its closely related species G. max. Of the studied species, A. hypogaea contained the most (246) AhFAR1 genes, which can be classified into three subgroups based on phylogenic relationships. The synonymous (Ks) and non-synonymous (Ka) substitution rates, phylogenetic relationship and synteny analysis of the FAR1 family provided deep insight into polyploidization, evolution and domestication of peanut AhFAR1 genes. The transcriptome data showed that the AhFAR1 genes exhibited distinct tissue- and stage-specific expression patterns in peanut. Three candidate genes including Ahy_A10g049543, Ahy_A06g026579, and Ahy_A10g048401, specifically expressed in peg and pod, might participate in pod development in the peanut. The quantitative real-time PCR (qRT-PCR) analyses confirmed that the three selected genes were highly and specifically expressed in the peg and pod. This study systematically analyzed gene structure, evolutionary characteristics and expression patterns of FAR1 gene family, which will provide a foundation for the study of genetic and biological function in the future.

AhNPR3 regulates the expression of WRKY and PR genes, and mediates the immune response of the peanut (Arachis hypogaea L.)

Systemic acquired resistance is an essential immune response that triggers a broad-spectrum disease resistance throughout the plant. In the present study, we identified a peanut lesion mimic mutant m14 derived from an ethyl methane sulfonate-mutagenized mutant pool of peanut cultivar "Yuanza9102." Brown lesions were observed in the leaves of an m14 mutant from seedling stage to maturity. Using MutMap together with bulked segregation RNA analysis approaches, a G-to-A point mutation was identified in the exon region of candidate gene Arahy.R60CUW, which is the homolog of AtNPR3 (Nonexpresser of PR genes) in Arabidopsis. This point mutation caused a transition from Gly to Arg within the C-terminal transactivation domain of AhNPR3A. The mutation of AhNPR3A showed no effect in the induction of PR genes when treated with salicylic acid. Instead, the mutation resulted in upregulation of WRKY genes and several PR genes, including pathogenesis-related thaumatin- and chitinase-encoding genes, which is consistent with the resistant phenotype of m14 to leaf spot disease. Further study on the AhNPR3A gene will provide valuable insights into understanding the molecular mechanism of systemic acquired resistance in peanut. Moreover, our results indicated that a combination of MutMap and bulked segregation RNA analysis is an effective method for identifying genes from peanut mutants.

BSA‑seq and genetic mapping reveals AhRt2 as a candidate gene responsible for red testa of peanut

The candidate recessive gene AhRt2 responsible for red testa of peanut was identified through combined BSA-seq and linkage mapping approaches. The testa color of peanuts (Arachis hypogaea L.) is an important trait, and those with red testa are particularly popular owing to the high-anthocyanin content. However, the identification of genes underlying the regulation of the red testa trait in peanut are rarely reported. In order to fine map red testa gene, two F_2:4 populations were constructed through the cross of YZ9102 (pink testa) with ZH12 (red testa) and ZH2 (red testa). Genetic analysis indicated that red testa was controlled by a single recessive gene named as AhRt2 (Red testa gene 2). Using BSA-seq approach, AhRt2 was preliminary identified on chromosome 12, which was further mapped to a 530-kb interval using 220 recombinant lines through linkage mapping. Furthermore, functional annotation, expression profiling, and the analyses of sequence variation confirmed that the anthocyanin reductase namely (Arahy.IK60LM) was the most likely candidate gene for AhRt2. It was found that a SNP in the third exon of AhRt2 altered the encoding amino acids, and was associated with red testa in peanut. In addition, a closely linked molecular marker linked with red testa trait in peanut was also developed for future studies. Our results provide valuable insight into the molecular mechanism underlying peanut testa color and present significant diagnostic marker resources for marker-assisted selected breeding in peanut.

Dissection of valine-glutamine genes and their responses to drought stress in Arachis hypogaea cv. Tifrunner

Valine-glutamine sequences (VQs) interact with WRKY transcription factors (TFs), forming VQ-WRKY protein complexes crucial for plant development and response to environmental changes. Cultivated peanut (Arachis hypogaea) is a tetraploid from A. duranensis and A. ipaensis cross. The Arachis spp. WRKY TFs have been identified, but Arachis VQs are largely unknown. This study identified VQs in A. duranensis, A. ipaensis, A. monticola, A. hypogaea cv. Fuhuasheng, A. hypogaea cv. Shitouqi, and A. hypogaea cv. Tifrunner. The study analyzed the homologous relationships between VQs in these Arachis spp. The VQ drought-tolerant genes were detected and VQ-WRKY interactions were determined in A. hypogaea cv. Tifrunner. The results showed that tetraploid Arachis spp. retained duplicated VQs, but lost ancestral VQs after allopolyploidization. The number of VQs in A. monticola, A. hypogaea cv. Fuhuasheng, and A. hypogaea cv. Shitouqi increased relative to their diploid ancestors. RNA-seq and quantitative real-time PCR experiments confirmed that three AhTVQs tolerate drought stress in A. hypogaea cv. Tifrunner. However, evidence of VQ-WRKY interaction for drought stress response is lacking in A. hypogaea cv. Tifrunner. Nevertheless, this study identified VQ-WRKY interactions, which possibly have multiple functions in A. hypogaea cv. Tifrunner. Altogether, this study dissected Arachis VQs, providing insights into Arachis VQ evolution and drought function.

Insights into the Genomic Architecture of Seed and Pod Quality Traits in the U.S. Peanut Mini-Core Diversity Panel

Traits such as seed weight, shelling percent, percent sound mature kernels, and seed dormancy determines the quality of peanut seed. Few QTL (quantitative trait loci) studies using biparental mapping populations have identified QTL for seed dormancy and seed grade traits. Here, we report a genome-wide association study (GWAS) to detect marker–trait associations for seed germination, dormancy, and seed grading traits in peanut. A total of 120 accessions from the U.S. peanut mini-core collection were evaluated for seed quality traits and genotyped using Axiom SNP (single nucleotide polymorphism) array for peanut. We observed significant variation in seed quality traits in different accessions and different botanical varieties. Through GWAS, we were able to identify multiple regions associated with sound mature kernels, seed weight, shelling percent, seed germination, and dormancy. Some of the genomic regions that were SNP associated with these traits aligned with previously known QTLs. For instance, QTL for seed dormancy has been reported on chromosome A05, and we also found SNP on the same chromosome associated with seed dormancy, explaining around 20% of phenotypic variation. In addition, we found novel genomic regions associated with seed grading, seed germination, and dormancy traits. SNP markers associated with seed quality and dormancy identified here can accelerate the selection process. Further, exploring the function of candidate genes identified in the vicinity of the associated marker will assist in understanding the complex genetic network that governs seed quality.

Insights into the Ancient Adaptation to Intertidal Environments by Red Algae Based on a Genomic and Multiomics Investigation of Neoporphyra haitanensis

Colonization of land from marine environments was a major transition for biological life on Earth, and intertidal adaptation was a key evolutionary event in the transition from marine- to land-based lifestyles. Multicellular intertidal red algae exhibit the earliest, systematic, and successful adaptation to intertidal environments, with Porphyra sensu lato (Bangiales, Rhodophyta) being a typical example. Here, a chromosome-level 49.67 Mb genome for Neoporphyra haitanensis comprising 9,496 gene loci is described based on metagenome-Hi-C-assisted whole-genome assembly, which allowed the isolation of epiphytic bacterial genome sequences from a seaweed genome for the first time. The compact, function-rich N. haitanensis genome revealed that ancestral lineages of red algae share common horizontal gene transfer events and close relationships with epiphytic bacterial populations. Specifically, the ancestor of N. haitanensis obtained unique lipoxygenase family genes from bacteria for complex chemical defense, carbonic anhydrases for survival in shell-borne conchocelis lifestyle stages, and numerous genes involved in stress tolerance. Combined proteomic, transcriptomic, and metabolomic analyses revealed complex regulation of rapid responses to intertidal dehydration/rehydration cycling within N. haitanensis. These adaptations include rapid regulation of its photosynthetic system, a readily available capacity to utilize ribosomal stores, increased methylation activity to rapidly synthesize proteins, and a strong anti-oxidation system to dissipate excess redox energy upon exposure to air. These novel insights into the unique adaptations of red algae to intertidal lifestyles inform our understanding of adaptations to intertidal ecosystems and the unique evolutionary steps required for intertidal colonization by biological life.

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.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

Genome-wide identification and expression analysis of auxin response factors in peanut (Arachis hypogaea L.)

Auxin response factors (ARFs) are transcription factors that regulate the expression of auxin response genes, and have important functions in plant growth and development. In this study, available genome data for peanut (Arachis hypogaea L.) were used to identify AhARF genes. In total, 61 AhARFs and 23 AtARFs were divided into six groups (I-VI). Molecular structural analysis revealed that the protein members of AhARF contain at least two domains, the B3 domain and the Auxin-resp domain, and that some have a C-terminal dimerisation domain. Screening of the transcriptome data of 22 tissues of A. hypogaea cv. Tifrunner in a public database showed high expression levels of AhARF2 and AhARF6. AhARF6 was expressed more highly in the stem and branch than in the root and leaf of the wild species Arachis monticola (A. mon) and cultivated species H103. After treatment with exogenous auxin (NAA), the expression of AhARF6 was inhibited, and this inhibition was greater in A. mon than in H103. The transcriptome map revealed that the expression of AhARF6 was higher in the larger pods of H8107 and ZP06 than in the medium pods of H103 and small pods of A. mon. Moreover, AhARF6-5 was proven to be localised in the nucleus, consistent with the location of AtARF6. These results suggest that AhARF6 may play an important role in pod development in peanut.

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.

Genotyping-by-Sequencing Based Genetic Mapping Identified Major and Consistent Genomic Regions for Productivity and Quality Traits in Peanut

With an objective of identifying the genomic regions for productivity and quality traits in peanut, a recombinant inbred line (RIL) population developed from an elite variety, TMV 2 and its ethyl methane sulfonate (EMS)-derived mutant was phenotyped over six seasons and genotyped with genotyping-by-sequencing (GBS), Arachis hypogaea transposable element (AhTE) and simple sequence repeats (SSR) markers. The genetic map with 700 markers spanning 2,438.1 cM was employed for quantitative trait loci (QTL) analysis which identified a total of 47 main-effect QTLs for the productivity and oil quality traits with the phenotypic variance explained (PVE) of 10-52% over the seasons. A common QTL region (46.7-50.1 cM) on Ah02 was identified for the multiple traits, such as a number of pods per plant (NPPP), pod weight per plant (PWPP), shelling percentage (SP), and test weight (TW). Similarly, a QTL (7.1-18.0 cM) on Ah16 was identified for both SP and protein content (PC). Epistatic QTL (epiQTL) analysis revealed intra- and inter-chromosomal interactions for the main-effect QTLs and other genomic regions governing these productivity traits. The markers identified by a single marker analysis (SMA) mapped to the QTL regions for most of the traits. Among the five potential candidate genes identified for PC, SP and oil quality, two genes (Arahy.7A57YA and Arahy.CH9B83) were affected by AhMITE1 transposition, and three genes (Arahy.J5SZ1I, Arahy.MZJT69, and Arahy.X7PJ8H) involved functional single nucleotide polymorphisms (SNPs). With major and consistent effects, the genomic regions, candidate genes, and the associated markers identified in this study would provide an opportunity for gene cloning and genomics-assisted breeding for increasing the productivity and enhancing the quality of peanut.

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.

Subgenome evolution in allotetraploid plants

Polyploidization is a well-known speciation and adaptation mechanism. Traces of former polyploidization events were discovered within many genomes, and especially in plants. Allopolyploidization by interspecific hybridization between two species is common. Among hybrid plants, many are domesticated species of agricultural interest and many of their genomes and of their presumptive parents have been sequenced. Hybrid genomes remain challenging to analyse because of the presence of multiple subgenomes. The genomes of hybrids often undergo rearrangement and degradation over time. Based on 10 hybrid plant genomes from six different genera, with hybridization dating from 10,000 to 5 million years ago, we assessed subgenome degradation, subgenomic intermixing and biased subgenome fractionation. The restructuring of hybrid genomes does not proceed proportionally with the age of the hybrid. The oldest hybrids in our data set display completely different fates: whereas the subgenomes of the tobacco plant Nicotiana benthamiana are in an advanced stage of degradation, the subgenomes of quinoa (Chenopodium quinoa) are exceptionally well conserved by structure and sequence. We observed statistically significant biased subgenome fractionation in seven out of 10 hybrids, which had different ages and subgenomic intermixing levels. Hence, we conclude that no correlation exists between biased fractionation and subgenome intermixing. Lastly, domestication may encourage or hinder subgenome intermixing, depending on the evolutionary context. In summary, comparative analysis of hybrid genomes and their presumptive parents allowed us to determine commonalities and differences between their evolutionary fates. In order to facilitate the future analysis of further hybrid genomes, we automated the analysis steps within manticore, which is publicly available at https://github.com/MatteoSchiavinato/manticore.git.

The chromosome-level Hemerocallis citrina Borani genome provides new insights into the rutin biosynthesis and the lack of colchicine

Hemerocallis citrina Borani (huang hua cai in Chinese) is an important horticultural crop whose flower buds are widely consumed as a delicious vegetable in Asia. Here we assembled a high-quality reference genome of H. citrina using single-molecule sequencing and Hi-C technologies. The genome assembly was 3.77 Gb and consisted of 3183 contigs with a contig N50 of 2.09 Mb, which were further clustered into 11 pseudochromosomes. A larger portion (3.25 Gb or 86.20%) was annotated as a repetitive content and 54,295 protein-coding genes were annotated in the genome. Genome evolution analysis showed that H. citrina experienced a recent whole-genome duplication (WGD) event at ~15.73 million years ago (Mya), which was the main factor leading to many multiple copies of orthologous genes. We used this reference genome to predict 20 genes involved in the rutin biosynthesis pathway. Moreover, our metabolomics data revealed neither colchicine nor its precursors in H. citrina, challenging the long-standing belief that this alkaloid causes poisoning by the plant. The results of our disruptive research are further substantiated by our genomic finding that H. citrina does not contain any genes involved in colchicine biosynthesis. The high-quality genome lays a solid foundation for genetic research and molecular breeding of H. citrina.

De novo full length transcriptome analysis of Arachis glabrata provides insights into gene expression dynamics in response to biotic and abiotic stresses

The perennial ornamental peanut Arachis glabrata represents one of the most adaptable wild Arachis species. This study used PacBio combined with BGISEQ-500 RNA-seq technology to study the transcriptome and gene expression dynamics of A. glabrata. Of the total 109,747 unique transcripts obtained, >90,566 transcripts showed significant homology to known proteins and contained the complete coding sequence (CDS). RNA-seq revealed that 1229, 1039, 1671, 3923, 1521 and 1799 transcripts expressed specifically in the root, stem, leaf, flower, peg and pod, respectively. We also identified thousands of differentially expressed transcripts in response to drought, salt, cold and leaf spot disease. Furthermore, we identified 30 polyphenol oxidase encoding genes associated with the quality of forage, making A. glabrata suitable as a forage crop. Our findings presented the first transcriptome study of A. glabrata which will facilitate genetic and genomics studies and lays the groundwork for a deeper understanding of the A. glabrata genome.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

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

Seed size/weight, a key domestication trait, is also an important selection target during peanut breeding. However, the mechanisms that regulate peanut seed development are unknown. We re-sequenced 12 RNA samples from developing seeds of two cultivated peanut accessions (Lines 8106 and 8107) and wild Arachis monticola at 15, 30, 45, and 60 days past flowering (DPF). Transcriptome analyses showed that ∼36,000 gene loci were expressed in each of the 12 RNA samples, with nearly half exhibiting moderate (2 ≤ FPKM < 10) expression levels. Of these genes, 12.2% (4,523) were specifically expressed during seed development, mainly at 15 DPF. Also, ∼12,000 genes showed significant differential expression at 30, 45, and/or 60 DPF within each of the three peanut accessions, accounting for 31.8-34.1% of the total expressed genes. Using a method that combined comprehensive transcriptome analysis and previously mapped QTLs, we identified several candidate genes that encode transcription factor TGA7, topless-related protein 2, IAA-amino acid hydrolase ILR1-like 5, and putative pentatricopeptide repeat-containing (PPR) protein. Based on sequence variations identified in these genes, SNP markers were developed and used to genotype both 30 peanut landraces and a genetic segregated population, implying that EVM0025654 encoding a PPR protein may be associated with the increased seed size/weight of the cultivated accessions in comparison with the allotetraploid wild peanut. Our results provide additional knowledge for the identification and functional research into candidate genes responsible for the seed size/weight phenotype in peanut.

AhABI4s Negatively Regulate Salt-Stress Response in Peanut

Soil salinity is one of the major factors that limit the area of cultivable land and yield potential of crops. The ability of salt tolerance varies with plant species. Peanut (Arachis hypogaea L.) is a moderately salt-sensitive and economically important crop, however, their biological processes involved in salt-stress response remain unclear. In this study, we investigated the role of A. hypogaea L. ABSCISIC ACID INSENSITIVE 4s (AhABI4s) in salt tolerance and elucidated its mode of action in peanuts. The results showed that the downregulation of AhABI4s via whole plant virus-induced gene silencing has enhanced the survival rate, biomass accumulation, and root/shoot ratio of peanut seedlings in response to salt-stress. Transcriptomics, quantitative proteomics, and phosphoproteomic analyses were performed using AhABI4s-silenced and Mock plants. The expression pattern of 15,247 genes, 1,900 proteins, and 2,620 phosphorylation sites were affected by silencing of AhABI4s in peanut leaf and root after sodium chloride (NaCl) treatment. Among them, 63 potential downstream target genes of ABI4 changed consistently at both transcription and translation levels, and the protein/phosphorylation levels of 31 ion transporters/channels were also affected. Electrophoretic mobility shift assays (EMSA) showed that ABI4 was able to bind to the promoters of HSP70, fructokinase (FRK), and pyruvate kinase (PK) coding genes in vitro. In addition, we also detected a binding preference of AhABI4 for CACT(G/T)GCA motif in the promoters of down-regulated genes in peanut leaf. Collectively, the potential downstream targets which were regulated at the levels of transcription and translation, binding preference, and in vivo phosphorylation sites that had been revealed in this study will provide new insight into the AhABI4s-mediated salt tolerance regulation mechanism in peanuts.

Construction of Pseudomolecules for the Chinese Chestnut (Castanea mollissima) Genome

The Chinese chestnut (Castanea mollissima Bl.) is a woody nut crop with a high ecological value. Although many cultivars have been selected from natural seedlings, elite lines with comprehensive agronomic traits and characters remain rare. To explore genetic resources with aid of whole genome sequence will play important roles in modern breeding programs for chestnut. In this study, we generated a high-quality C. mollissima genome assembly by combining 90× Pacific Biosciences long read and 170× high-throughput chromosome conformation capture data. The assembly was 688.93 Mb in total, with a contig N50 of 2.83 Mb. Most of the assembled sequences (99.75%) were anchored onto 12 chromosomes, and 97.07% of the assemblies were accurately anchored and oriented. A total of 33,638 protein-coding genes were predicted in the C. mollissima genome. Comparative genomic and transcriptomic analyses provided insights into the genes expressed in specific tissues, as well as those associated with burr development in the Chinese chestnut. This highly contiguous assembly of the C. mollissima genome provides a valuable resource for studies aiming at identifying and characterizing agronomical-important traits, and will aid the design of breeding strategies to develop more focused, faster, and predictable improvement programs.

Comprehensive genomic characterization of NAC transcription factor family and their response to salt and drought stress in peanut

BACKGROUND

Peanut is one of the most important oil crop species worldwide. NAC transcription factor (TF) genes play important roles in the salt and drought stress responses of plants by activating or repressing target gene expression. However, little is known about NAC genes in peanut.

RESULTS

We performed a genome-wide characterization of NAC genes from the diploid wild peanut species Arachis duranensis and Arachis ipaensis, which included analyses of chromosomal locations, gene structures, conserved motifs, expression patterns, and cis-acting elements within their promoter regions. In total, 81 and 79 NAC genes were identified from A. duranensis and A. ipaensis genomes. Phylogenetic analysis of peanut NACs along with their Arabidopsis and rice counterparts categorized these proteins into 18 distinct subgroups. Fifty-one orthologous gene pairs were identified, and 46 orthologues were found to be highly syntenic on the chromosomes of both A. duranensis and A. ipaensis. Comparative RNA sequencing (RNA-seq)-based analysis revealed that the expression of 43 NAC genes was up- or downregulated under salt stress and under drought stress. Among these genes, the expression of 17 genes in cultivated peanut (Arachis hypogaea) was up- or downregulated under both stresses. Moreover, quantitative reverse transcription PCR (RT-qPCR)-based analysis revealed that the expression of most of the randomly selected NAC genes tended to be consistent with the comparative RNA-seq results.

CONCLUSION

Our results facilitated the functional characterization of peanut NAC genes, and the genes involved in salt and drought stress responses identified in this study could be potential genes for peanut improvement.