Single-nucleus RNA and ATAC sequencing analyses provide molecular insights into early pod development of peanut fruit.

Peanut (Arachis hypogaea L.) is an important leguminous oil and economic crop that produces flowers aboveground and fruits belowground. Subterranean fruit-pod development, which significantly affects peanut production, involves complex molecular mechanisms that likely require the coordinated regulation of multiple genes in different tissues. To investigate the molecular mechanisms that underlie peanut fruit-pod development, we characterized the anatomical features of early fruit-pod development and integrated single-nucleus RNA-sequencing (snRNA-seq) and single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) data at the single-cell level. We identified distinct cell types, such as meristem, embryo, vascular tissue, cuticular layer, and stele cells within the shell wall. These specific cell types were used to examine potential molecular changes unique to each cell type during pivotal stages of fruit-pod development. snRNA-seq analyses of differentially expressed genes revealed cell-type-specific insights that were not previously obtainable from transcriptome analyses of bulk RNA. For instance, we identified MADS-box genes that contributes to the formation of parenchyma cells and gravity-related genes that are present in the vascular cells, indicating an essential role for the vascular cells in peg gravitropism. Overall, our single-nucleus analysis provides comprehensive and novel information on specific cell types, gene expression, and chromatin accessibility during the early stages of fruit-pod development. This information will enhance our understanding of the mechanisms that underlie fruit-pod development in peanut and contribute to efforts aimed at improving peanut production.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis.

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

High-density bin-based genetic map reveals a 530-kb chromosome segment derived from wild peanut contributing to late leaf spot resistance.

KEY MESSAGE: Twenty-eight QTLs for LLS disease resistance were identified using an amphidiploid constructed mapping population, a favorable 530-kb chromosome segment derived from wild species contributes to the LLS resistance. Late leaf spot (LLS) is one of the major foliar diseases of peanut, causing serious yield loss and affecting the quality of kernel and forage. Some wild Arachis species possess higher resistance to LLS as compared with cultivated peanut; however, ploidy level differences restrict utilization of wild species. In this study, a synthetic amphidiploid (Ipadur) of wild peanuts with high LLS resistance was used to cross with Tifrunner to construct TI population. In total, 200 recombinant inbred lines were collected for whole-genome resequencing. A high-density bin-based genetic linkage map was constructed, which includes 4,809 bin markers with an average inter-bin distance of 0.43 cM. The recombination across cultivated and wild species was unevenly distributed, providing a novel recombination landscape for cultivated-wild Arachis species. Using phenotyping data collected across three environments, 28 QTLs for LLS disease resistance were identified, explaining 4.35-20.42% of phenotypic variation. The major QTL located on chromosome 14, qLLS14.1, could be consistently detected in 2021 Jiyang and 2022 Henan with 20.42% and 12.12% PVE, respectively. A favorable 530-kb chromosome segment derived from Ipadur was identified in the region of qLLS14.1, in which 23 disease resistance proteins were located and six of them showed significant sequence variations between Tifrunner and Ipadur. Allelic variation analysis indicating the 530-kb segment of wild species might contribute to the disease resistance of LLS. These associate genomic regions and candidate resistance genes are of great significance for peanut breeding programs for bringing durable resistance through pyramiding such multiple LLS resistance loci into peanut cultivars.

Corrigendum: Triticeae crop genome biology: an endless frontier.

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

Genome-wide analysis of PIN genes in cultivated peanuts (Arachis hypogaea L.): identification, subcellular localization, evolution, and expression patterns.

BACKGROUND: Auxin is an important hormone in plants and the PIN-FORMED (PIN) genes are essential to auxin distribution in growth and developmental processes of plants. Peanut is an influential cash crop, but research into PIN genes in peanuts remains limited. RESULTS: In this study, 16 PIN genes were identified in the genome of cultivated peanut, resolving into four subfamilies. All PIN genes were predicted to be located in the plasma membrane and a subcellular location experiment confirmed this prediction for eight of them. The gene structure, cis-elements in the promoter, and evolutionary relationships were elucidated, facilitating our understanding of peanut PINs and their evolution. In addition, the expression patterns of these PINs in various tissues were analyzed according to a previously published transcriptome dataset and qRT-PCR, which gave us a clear understanding of the temporal and spatial expression of PIN genes in different growth stages and different tissues. The expression trend of homologous genes was similar. AhPIN2A and AhPIN2B exhibited predominant expression in roots. AhPIN1A-1 and AhPIN1B-1 displayed significant upregulation following peg penetration, suggesting a potential close association with peanut pod development. Furthermore, we presented the gene network and gene ontology enrichment of these PINs. Notably, AhABCB19 exhibited a co-expression relationship with AhPIN1A and AhPIN1B-1, with all three genes displaying higher expression levels in peanut pegs and pods. These findings reinforce their potential role in peanut pod development. CONCLUSIONS: This study details a comprehensive analysis of PIN genes in cultivated peanuts and lays the foundation for subsequent studies of peanut gene function and phenotype.

Triticeae crop genome biology: an endless frontier.

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

BarleyExpDB: an integrative gene expression database for barley.

BACKGROUND: RNA-sequencing (RNA-seq) has been widely used to study the dynamic expression patterns of transcribed genes, which can lead to new biological insights. However, processing and analyzing these huge amounts of histological data remains a great challenge for wet labs and field researchers who lack bioinformatics experience and computational resources. RESULTS: We present BarleyExpDB, an easy-to-operate, free, and web-accessible database that integrates transcriptional profiles of barley at different growth and developmental stages, tissues, and stress conditions, as well as differential expression of mutants and populations to build a platform for barley expression and visualization. The expression of a gene of interest can be easily queried by searching by known gene ID or sequence similarity. Expression data can be displayed as a heat map, along with functional descriptions as well as Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Proteins Families Database, and Simple Modular Architecture Research Tool annotations. CONCLUSIONS: BarleyExpDB will serve as a valuable resource for the barley research community to leverage the vast publicly available RNA-seq datasets for functional genomics research and crop molecular breeding.

An R-R-type MYB transcription factor promotes non-climacteric pepper fruit carotenoid pigment biosynthesis.

Carotenoids are major accessory pigments in the chloroplast, and they also act as phytohormones and volatile compound precursors to influence plant development and confer characteristic colours, affecting both the aesthetic and nutritional value of fruits. Carotenoid pigmentation in ripening fruits is highly dependent on developmental trajectories. Transcription factors incorporate developmental and phytohormone signalling to regulate the biosynthesis process. By contrast to the well-established pathways regulating ripening-related carotenoid biosynthesis in climacteric fruit, carotenoid regulation in non-climacteric fruit is poorly understood. Capsanthin is the primary carotenoid of non-climacteric pepper (Capsicum) fruit; its biosynthesis is tightly associated with fruit ripening, and it confers red pigmentation to the ripening fruit. In the present study, using a coexpression analysis, we identified an R-R-type MYB transcription factor, DIVARICATA1, and demonstrated its role in capsanthin biosynthesis. DIVARICATA1 encodes a nucleus-localised protein that functions primarily as a transcriptional activator. Functional analyses showed that DIVARICATA1 positively regulates carotenoid biosynthetic gene (CBG) transcript levels and capsanthin levels by directly binding to and activating CBG promoter transcription. Furthermore, an association analysis revealed a significant positive association between DIVARICATA1 transcription level and capsanthin content. ABA promotes capsanthin biosynthesis in a DIVARICATA1-dependent manner. Comparative transcriptomic analysis of DIVARICATA1 in Solanaceae plants showed that its function likely differs among species. Moreover, the pepper DIVARICATA1 gene could be regulated by the ripening regulator MADS-RIN. The present study illustrates the transcriptional regulation of capsanthin biosynthesis and offers a target for breeding peppers with high red colour intensity.

BGFD: an integrated multi-omics database of barley gene families.

BACKGROUND: A gene family comprises a group of genes with similar functional domains that play various roles in plant growth, development, and responses to environmental stimuli. Barley (Hordeum vulgare L.) is the fourth most cultivated cereal crop worldwide, and it is an important model species for genetic studies. Systematic identification and annotation of gene families are key for studies of molecular function and evolutionary history. RESULTS: We constructed a multi-omics database containing 5593 genes of 77 gene families called the Barley Gene Family Database (BGFD: http://barleygfdb.com ). BGFD is a free, user-friendly, and web-accessible platform that provides data on barley family genes. BGFD provides intuitive visual displays to facilitate studies of the physicochemical properties, gene structure, phylogenetic relationships, and motif organization of genes. Massive multi-omics datasets have been acquired and processed to generate an atlas of expression pattern profiles and genetic variation in BGFD. The platform offers several practical toolkits to conduct searches, browse, and employ BLAST functions, and the data are downloadable. CONCLUSIONS: BGFD will aid research on the domestication and adaptive evolution of barley; it will also facilitate the screening of candidate genes and exploration of important agronomic traits in barley.

Analysis of the Transcriptional Dynamics of Regulatory Genes During Peanut Pod Development Caused by Darkness and Mechanical Stress.

Peanut is an oil crop with important economic value that is widely cultivated around the world. It blooms on the ground but bears fruit underground. When the peg penetrates the ground, it enters a dark environment, is subjected to mechanical stress from the soil, and develops into a normal pod. When a newly developed pod emerges from the soil, it turns green and stops growing. It has been reported that both darkness and mechanical stress are necessary for normal pod development. In this study, we investigated changes in gene expression during the reverse process of peg penetration: developmental arrest caused by pod (Pattee 3 pods) excavation. Bagging the aerial pods was used to simulate loss of mechanical pressure, while direct exposure of the aerial pods was used to simulate loss of both mechanical pressure and darkness. After the loss of mechanical stress and darkness, the DEGs were significantly enriched in photosynthesis, photosynthesis-antenna proteins, plant-pathogen interaction, DNA replication, and circadian rhythm pathways. The DNA replication pathway was enriched by down-regulated genes, and the other four pathways were enriched by upregulated genes. Upregulated genes were also significantly enriched in protein ubiquitination and calmodulin-related genes, highlighting the important role of ubiquitination and calcium signaling in pod development. Further analysis of DEGs showed that phytochrome A (Phy A), auxin response factor 9 (IAA9), and mechanosensitive ion channel protein played important roles in geocarpy. The expression of these two genes increased in subterranean pods but decreased in aerial pods. Based on a large number of chloroplast-related genes, calmodulin, kinases, and ubiquitin-related proteins identified in this study, we propose two possible signal transduction pathways involved in peanut geocarpy, namely, one begins in chloroplasts and signals down through phosphorylation, and the other begins during abiotic stress and signals down through calcium signaling, phosphorylation, and ubiquitination. Our study provides valuable information about putative regulatory genes for peanut pod development and contributes to a better understanding of the biological phenomenon of geocarpy.

Genomic insights into positive selection during barley domestication.

BACKGROUND: Cultivated barley (Hordeum vulgare) is widely used in animal feed, beverages, and foods and has become a model crop for molecular evolutionary studies. Few studies have examined the evolutionary fates of different types of genes in barley during the domestication process. RESULTS: The rates of nonsynonymous substitution (Ka) to synonymous substitution (Ks) were calculated by comparing orthologous genes in different barley groups (wild vs. landrace and landrace vs. improved cultivar). The rates of evolution, properties, expression patterns, and diversity of positively selected genes (PSGs) and negatively selected genes (NSGs) were compared. PSGs evolved more rapidly, possessed fewer exons, and had lower GC content than NSGs; they were also shorter and had shorter intron, exon, and first exon lengths. Expression levels were lower, the tissue specificity of expression was higher, and codon usage bias was weaker for PSGs than for NSGs. Nucleotide diversity analysis revealed that PSGs have undergone a more severe genetic bottleneck than NSGs. Several candidate PSGs were involved in plant growth and development, which might make them as excellent targets for the molecular breeding of barley. CONCLUSIONS: Our comprehensive analysis of the evolutionary, structural, and functional divergence between PSGs and NSGs in barley provides new insight into the evolutionary trajectory of barley during domestication. Our findings also aid future functional studies of PSGs in barley.

The complete chloroplast genome and phylogenetic analysis of Astragalus sinicus Linne 1767.

Astragalus sinicus Linne 1767 is a traditional winter-growing green manure, that plays an important role in upgrading soil fertility and maintaining crop yield and quality for rice fields. This study reports the complete chloroplast genome of A. sinicus. The chloroplast genome contained 110 complete genes, including 76 protein-coding genes, 4 ribosomal RNA genes, and 30 tRNA genes with 123,830 bp in length and a 34.66% GC content with IR loss. The evolutionary history, referred to as the maximum-likelihood (ML), showed that A. sinicus and Astragalus bhotanensis were most closely related. The chloroplast genome analysis of A. sinicus will serve as a reference for future studies on species evolution, plant conservation, and molecular phylogeny in Astragalus.

Genome-Wide Analysis and Expression Profiles of Ethylene Signal Genes and Apetala2/Ethylene-Responsive Factors in Peanut (Arachis hypogaea L.).

Peanut is an important oil and economic crop widely cultivated in the world. It has special characteristics such as blooming on the ground but bearing fruits underground. During the peg penetrating into the ground, it is subjected to mechanical stress from the soil at the same time. It has been proved that mechanical stress affects plant growth and development by regulating the ethylene signaling-related genes. In this study, we identified some genes related to ethylene signal of peanut, including 10 ethylene sensors, two constitutive triple responses (CTRs), four ethylene insensitive 2 (EIN2s), four ethylene insensitive 3 (EIN3s), six EIN3-binding F-box proteins (EBFs), and 188 Apetala2/ethylene-responsive factors (AP2/ERFs). One hundred and eighty-eight AP2/ERFs were further divided into four subfamilies, 123 ERFs, 56 AP2s, 6 Related to ABI3/VP1 (RAVs), and three Soloists, of them one hundred and seventy AP2/ERF gene pairs were clustered into segmental duplication events in genome of Arachis hypogaea. A total of 134, 138, 97, and 150 AhAP2/ERF genes formed 210, 195, 166, and 525 orthologous gene pairs with Arachis duranensis, Arachis ipaensis, Arabidopsis thaliana, and Glycine max, respectively. Our transcriptome results showed that two EIN3s (Arahy.J729H0 and Arahy.S7XF8N) and one EBFs (Arahy.G4JMEM) were highly expressed when mechanical stress increased. Among the 188 AhAP2/ERF genes, there were 31 genes with the fragments per kilobase of exon model per million mapped fragments (FPKM) ≥ 100 at least one of the 15 samples of Tifrunner. Among them, three AhAP2/ERFs (Arahy.15RATX, Arahy.FAI7YU, and Arahy.452FBF) were specifically expressed in seeds and five AhAP2/ERFs (Arahy.HGAZ7D, Arahy.ZW7540, Arahy.4XS3FZ, Arahy.QGFJ76, and Arahy.AS0C7C) were highly expressed in the tissues, which responded mechanical stress, suggesting that they might sense mechanical stress. Mechanical stress simulation experiment showed that three AhAP2/ERFs (Arahy.QGFJ76, Arahy.AS0C7C, and Arahy.HGAZ7D) were sensitive to mechanical stress changes and they all had the conservative repressor motif (DLNXXP) in the C-terminus, indicated that they might transmit mechanical stress signals through transcriptional inhibition. This study reveals the regulatory landscape of ethylene signal-related genes in peanut, providing valuable information for the mining of target genes for further study.

Identification of cold tolerance QTLs at the bud burst stage in 211 rice landraces by GWAS.

BACKGROUND: Rice is a crop that is very sensitive to low temperature, and its morphological development and production are greatly affected by low temperature. Therefore, understanding the genetic basis of cold tolerance in rice is of great significance for mining favorable genes and cultivating excellent rice varieties. However, there have been limited studies focusing on cold tolerance at the bud burst stage; therefore, considerable attention should be given to the genetic basis of cold tolerance at this stage. RESULTS: In this study, a natural population consisting of 211 rice landraces collected from 15 provinces in China and other countries was used for the first time to evaluate cold tolerance at the bud burst stage. Population structure analysis showed that this population was divided into two groups and was rich in genetic diversity. Our evaluation results confirmed that japonica rice was more tolerant to cold at the bud burst stage than indica rice. A genome-wide association study (GWAS) was performed with the phenotypic data of 211 rice landraces and a 36,727 SNP dataset under a mixed linear model. Twelve QTLs (P < 0.0001) were identified for the seedling survival rate (SR) after treatment at 4 °C, in which there were five QTLs (qSR2-2, qSR3-1, qSR3-2, qSR3-3 and qSR9) that were colocalized with those from previous studies and seven QTLs (qSR2-1, qSR3-4, qSR3-5, qSR3-6, qSR3-7, qSR4 and qSR7) that were reported for the first time. Among these QTLs, qSR9, harboring the most significant SNP, explained the most phenotypic variation. Through bioinformatics analysis, five genes (LOC_Os09g12440, LOC_Os09g12470, LOC_Os09g12520, LOC_Os09g12580 and LOC_Os09g12720) were identified as candidates for qSR9. CONCLUSION: This natural population consisting of 211 rice landraces combined with high-density SNPs will serve as a better choice for identifying rice QTLs/genes in the future, and the detected QTLs associated with cold tolerance at the bud burst stage in rice will be conducive to further mining favorable genes and breeding rice varieties under cold stress.

Genomic and Phenotypic Divergence in Wild Barley Driven by Microgeographic Adaptation.

Microgeographic adaptation is a fundamental driving force of evolution, but the underlying causes remain undetermined. Here, the phenotypic, genomic and transcriptomic variations of two wild barley populations collected from sharply divergent and adjacent micro-geographic sites to identify candidate genes associated with edaphic local adaptation are investigated. Common garden and reciprocal transplant studies show that large phenotypic differentiation and local adaptation to soils occur between these populations. Genetic, phylogenetic and admixture analyses based on population resequencing show that significant genetic divergences occur between basalt and chalk populations. These divergences are consistent with the phenotypic variations observed in the field. Genome sweep analyses reveal 162.7 Mb of selected regions driven by edaphic local adaptation, in which 445 genes identified, including genes associated with root architecture, metal transport/detoxification, and ABA signaling. When the phenotypic, genomic and transcriptomic data are combined, HvMOR, encoding an LBD transcription factor, is determined to be the vital candidate for regulating the root architecture to adapt to edaphic conditions at the microgeographic scale. This study provides new insights into the genetic basis of edaphic adaptation and demonstrates that edaphic factors may contribute to the evolution and speciation of barley.

Frequent intra- and inter-species introgression shapes the landscape of genetic variation in bread wheat.

BACKGROUND: Bread wheat is one of the most important and broadly studied crops. However, due to the complexity of its genome and incomplete genome collection of wild populations, the bread wheat genome landscape and domestication history remain elusive. RESULTS: By investigating the whole-genome resequencing data of 93 accessions from worldwide populations of bread wheat and its diploid and tetraploid progenitors, together with 90 published exome-capture data, we find that the B subgenome has more variations than A and D subgenomes, including SNPs and deletions. Population genetics analyses support a monophyletic origin of domesticated wheat from wild emmer in northern Levant, with substantial introgressed genomic fragments from southern Levant. Southern Levant contributes more than 676 Mb in AB subgenomes and enriched in the pericentromeric regions. The AB subgenome introgression happens at the early stage of wheat speciation and partially contributes to their greater genetic diversity. Furthermore, we detect massive alien introgressions that originated from distant species through natural and artificial hybridizations, resulting in the reintroduction of ~ 709 Mb and ~ 1577 Mb sequences into bread wheat landraces and varieties, respectively. A large fraction of these intra- and inter-introgression fragments are associated with quantitative trait loci of important traits, and selection events are also identified. CONCLUSION: We reveal the significance of multiple introgressions from distant wild populations and alien species in shaping the genetic components of bread wheat, and provide important resources and new perspectives for future wheat breeding.

Transcriptional Dynamics of Grain Development in Barley (Hordeum vulgare L.).

Grain development, as a vital process in the crop's life cycle, is crucial for determining crop quality and yield. However, the molecular basis and regulatory network of barley grain development is not well understood at present. Here, we investigated the transcriptional dynamics of barley grain development through RNA sequencing at four developmental phases, including early prestorage phase (3 days post anthesis (DPA)), late prestorage or transition phase (8 DPA), early storage phase (13 DPA), and levels off stages (18 DPA). Transcriptome profiling found that pronounced shifts occurred in the abundance of transcripts involved in both primary and secondary metabolism during grain development. The transcripts' activity was decreased during maturation while the largest divergence was observed between the transitions from prestorage phase to storage phase, which coincided with the physiological changes. Furthermore, the transcription factors, hormone signal transduction-related as well as sugar-metabolism-related genes, were found to play a crucial role in barley grain development. Finally, 4771 RNA editing events were identified in these four development stages, and most of the RNA editing genes were preferentially expressed at the prestore stage rather than in the store stage, which was significantly enriched in "essential" genes and plant hormone signal transduction pathway. These results suggested that RNA editing might act as a 'regulator' to control grain development. This study systematically dissected the gene expression atlas of barley grain development through transcriptome analysis, which not only provided the potential targets for further functional studies, but also provided insights into the dynamics of gene regulation underlying grain development in barley and beyond.

Characterization of microRNAs and their targets in wild barley (Hordeum vulgare subsp. spontaneum) using deep sequencing.

MicroRNAs (miRNA) are a class of small, endogenous RNAs that play a negative regulatory role in various developmental and metabolic processes of plants. Wild barley (Hordeum vulgare subsp. spontaneum), as the progenitor of cultivated barley (Hordeum vulgare subsp. vulgare), has served as a valuable germplasm resource for barley genetic improvement. To survey miRNAs in wild barley, we sequenced the small RNA library prepared from wild barley using the Illumina deep sequencing technology. A total of 70 known miRNAs and 18 putative novel miRNAs were identified. Sequence analysis revealed that all of the miRNAs identified in wild barley contained the highly conserved hairpin sequences found in barley cultivars. MiRNA target predictions showed that 12 out of 52 miRNA families were predicted to target transcription factors, including 8 highly conserved miRNA families in plants and 4 wheat-barley conserved miRNA families. In addition to transcription factors, other predicted target genes were involved in diverse physiological and metabolic processes and stress defense. Our study for the first time reported the large-scale investigation of small RNAs in wild barley, which will provide essential information for understanding the regulatory role of miRNAs in wild barley and also shed light on future practical utilization of miRNAs for barley improvement.