Improved pea reference genome and pan-genome highlight genomic features and evolutionary characteristics

Perspectives and opportunities in forensic human, animal, and plant integrative genomics in the Pangenome era

The Human Pangenome Reference Consortium, the Chinese Pangenome Consortium, and other plant and animal pangenome projects have announced the completion of pilot work aimed at constructing high-quality, haplotype-resolved reference graph genomes representative of global ethno-linguistically different populations or different plant and animal species. These graph-based, gapless pangenome references, which are enriched in terms of genomic diversity, completeness, and contiguity, have the potential for enhancing long-read sequencing (LRS)-based genomic research, as well as improving mappability and variant genotyping on traditional short-read sequencing platforms. We comprehensively discuss the advancements in pangenome-based genomic integrative genomic discoveries across forensic-related species (humans, animals, and plants) and summarize their applications in variant identification and forensic genomics, epigenetics, transcriptomics, and microbiome research. Recent developments in multiplexed array sequencing have introduced a highly efficient and programmable technique to overcome the limitations of short forensic marker lengths in LRS platforms. This technique enables the concatenation of short RNA transcripts and DNA fragments into LRS-optimal molecules for sequencing, assembly, and genotyping. The integration of new pangenome reference coordinates and corresponding computational algorithms will benefit forensic integrative genomics by facilitating new marker identification, accurate genotyping, high-resolution panel development, and the updating of statistical algorithms. This review highlights the necessity of integrating LRS-based platforms, pangenome-based study designs, and graph-based pangenome references in short-read mapping and LRS-based innovations to achieve precision forensic science.

Boron controls apical dominance in Pea (Pisum sativum) via promoting polar auxin transport

Plant architecture and subsequent productivity are determined by the shoot apical dominance, which is disturbed by the deficiency of boron, one of the essential trace elements for plant growth and reproduction. However, the mechanism by which B controls shoot apical dominance or axillary bud outgrows under B deficiency is still unclear. This work aimed to investigate the mechanistic basis of this process, with focus on the interaction between B and polar auxin transport. Adopting an all-buds phenotyping methodology and employing several complementary approaches, we found that boron deficiency inhibited plant growth and changed the shoot architecture, resulting in the outgrowth of axillary buds at nodes 1-3. This was related to the auxin accumulation in shoot apical parts buds under B deficiency. Applying N-1-naphthylphthalamic acid to inhibit auxin transport from the shoot apex promoted the outgrowth of axillary buds in boron-sufficient (+B) plants. In decapitated plants, the application of exogenous auxin to the shoot apex only inhibited the outgrowth of axillary buds in +B plants. At higher auxin doses, the toxic effect of IAA was observed in the lower part of the shoot, which was more severe in +B plants than in B-deprived (-B) plants. Furthermore, the expression of PsPIN3 was significantly downregulated under -B conditions. These results indicate that B deficiency inhibits PAT from the apical bud through the main stem to the lower parts, leading to an increase of auxin level in the apical bud, which inhibits the growth of apical buds while stimulating the outgrowth of axillary buds.

Near-complete telomere-to-telomere de novo genome assembly in Egyptian clover (Trifolium alexandrinum)

Egyptian clover (Trifolium alexandrinum L.), also known as berseem clover, is an important forage crop to semi-arid conditions that was domesticated in ancient Egypt in 5,5000 BCE and introduced and well adapted to numerous countries including India, Pakistan, Turkey, and Mediterranean region. Despite its agricultural importance, genomic research on Egyptian clover has been limited to developing efficient modern breeding programs. In the present study, we constructed near-complete telomere-to-telomere-level genome assemblies for 2 Egyptian clover cultivars, Helaly and Fahl. Initial assemblies were established by using highly fidelity long-read technology. To extend sequence contiguity, we developed a gap-targeted sequencing (GAP-Seq) method, in which contig ends are targeted for sequencing to obtain long reads bridging 2 contigs. The total length of the resultant chromosome-level assemblies was 547.7 Mb for Helaly and 536.3 Mb for Fahl. These differences in sequence length can be attributed to the expansion of DNA transposons. Population genomic analysis using single-nucleotide polymorphisms revealed genomic regions highly differentiated between 2 cultivars and increased genetic uniformity within each cultivar. Gene ontologies associated with metabolic and biosynthetic processes and developmental processes were enriched in these genomic regions, indicating that these genes may determine the unique characteristics of each cultivar. Comprehensive genomic resources can provide valuable insights into genetic improvements in Egyptian clover and legume genomics.

Promoter of Vegetable Pea PsPIP2-4 Responds to Abiotic Stresses in Transgenic Tobacco

Plasma membrane intrinsic proteins (PIPs), one sub-family of aquaporins (AQPs), are responsible for plant abiotic stress responses. However, little information is currently available about the stress responsiveness of the PIP promoter in vegetable pea. In the present study, one novel promoter of PsPIP2-4 which shared high similarity to the PIP2-type AQPs from other plants, was isolated. Quantitative real-time PCR (qRT-PCR) assays suggested that PsPIP2-4 was predominantly expressed in leaves and abundantly induced by abiotic stress treatments (polyethylene glycol (PEG) 6000, NaCl, and methyl jasmonate (MeJA)). Further, the promoter activity of PsPIP2-4 was verified in transgenic tobacco plants. Beta-glucuronidase (GUS) staining driven by the PsPIP2-4 promoter confirmed that it was mainly detected in the leaves of transgenic seedlings, especially in the guard cells. Exposure of transgenic seedlings to various environmental stimuli proved that the promoter activity of PsPIP2-4 was abundantly strengthened by osmotic, salt, and MeJA stresses. This research provides one stress-inducible promoter enabling targeted gene expression under abiotic stresses and demonstrates its usefulness in the genetic improvement of plant stress resistance.

Genome-wide mapping of main histone modifications and coordination regulation of metabolic genes under salt stress in pea (Pisum sativum L)

Pea occupy a key position in modern biogenetics, playing multifaceted roles as food, vegetable, fodder, and green manure. However, due to the complex nature of its genome and the prolonged unveiling of high-quality genetic maps, research into the molecular mechanisms underlying pea development and stress responses has been significantly delayed. Furthermore, the exploration of its epigenetic modification profiles and associated regulatory mechanisms remains uncharted. This research conducted a comprehensive investigation of four specific histone marks, namely H3K4me3, H3K27me3, H3K9ac, and H3K9me2, and the transcriptome in pea under normal conditions, and established a global map of genome-wide regulatory elements, chromatin states, and dynamics based on these major modifications. Our analysis identified epigenomic signals across ~82.6% of the genome. Each modification exhibits distinct enrichment patterns: H3K4me3 is predominantly associated with the gibberellin response pathway, H3K27me3 is primarily associated with auxin and ethylene responses, and H3K9ac is primarily associated with negative regulatory stimulus responses. We also identified a novel bivalent chromatin state (H3K9ac-H3K27me3) in pea, which is related to their development and stress response. Additionally, we unveil that these histone modifications synergistically regulate metabolic-related genes, influencing metabolite production under salt stress conditions. Our findings offer a panoramic view of the major histone modifications in pea, elucidate their interplay, and highlight their transcriptional regulatory roles during salt stress.

Unveiling the Genomic Symphony: Identification Cultivar-Specific Genes and Enhanced Insights on Sweet Sorghum Genomes Through Comprehensive superTranscriptomic Analysis

Sorghum (Sorghum bicolor (L.) Moench) is a multipurpose crop grown for food, fodder, and bioenergy production. Its cultivated varieties, along with their wild counterparts, contribute to the core genetic pool. Despite the availability of several re-sequenced sorghum genomes, a variable portion of sorghum genomes is not reported during reference genome assembly and annotation. The present analysis used 223 publicly available RNA-seq datasets from seven sweet sorghum cultivars to construct superTranscriptome. This approach yielded 45,864 Representative Transcript Assemblies (RTAs) that showcased intriguing Presence/Absence Variation (PAV) across 15 published sorghum genomes. We found 301 superTranscripts were exclusive to sweet sorghum, including 58 de novo genes encoded core and linker histones, zinc finger domains, glucosyl transferases, cellulose synthase, etc. The superTranscriptome added 2,802 new protein-coding genes to the Sweet Sorghum Reference Genome (SSRG), of which 559 code for different transcription factors (TFs). Our analysis revealed that MULE-like transposases were abundant in the sweet sorghum genome and could play a hidden role in the evolution of sweet sorghum. We observed large deletions in the D locus and terminal deletions in four other NAC encoding loci in the SSRG compared to its wild progenitor (353) suggesting non-functional NAC genes contributed to trait development in sweet sorghum. Moreover, superTranscript-based methods for Differential Exon Usage (DEU) and Differential Gene Expression (DGE) analyses were more accurate than those based on the SSRG. This study demonstrates that the superTranscriptome can enhance our understanding of fundamental sorghum mechanisms, improve genome annotations, and potentially even replace the reference genome.

Landscape of structural variants reveals insights for local adaptations in the Asian corn borer

Capturing the genetic diversity of different wild populations is crucial for unraveling the mechanisms of adaptation and establishing links between genome evolution and local adaptation. The Asian corn borer (ACB) moth has undergone natural selection during its adaptative evolution. However, structural variants (SVs), which play significant roles in these adaptation processes, have not been previously identified. Here, we constructed a multi-assembly graph pan-genome to highlight the importance of SVs in local adaptation. Our analysis revealed that the graph pan-genome contained 176.60 Mb (∼37.33%) of unique sequences. Subsequently, we performed an analysis of expression quantitative trait loci (QTLs) to explore the impact of SVs on gene expression regulation. Notably, through QTL mapping analysis, we identified the FTZ-F1 gene as a potential candidate gene associated with the traits of larval development rate. In sum, we explored the impact of SVs on the local adaptation of pests, therefore facilitating accelerated pest management strategies.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Understanding the root of the problem for tackling pea root rot disease

Pea (Pisum sativum), a crop historically significant in the field of genetics, is regaining momentum in sustainable agriculture due to its high protein content and environmental benefits. However, its cultivation faces significant challenges from root rot, a complex disease caused by multiple soil-borne pathogens prevalent across most pea growing regions. This disease leads to substantial yield losses, further complicated by the dynamic interactions among pathogens, soil conditions, weather, and agricultural practices. Recent advancements in molecular diagnostics provide promising tools for the early and precise detection of these pathogens, which is critical for implementing effective disease management strategies. In this review, we explore how the availability of latest pea genomic resources and emerging technologies, such as CRISPR and cell-specific transcriptomics, will enable a deeper understanding of the molecular basis underlying host-pathogen interactions. We emphasize the need for a comprehensive approach that integrates genetic resistance, advanced diagnostics, cultural practices and the role of the soil microbiome in root rot. By leveraging these strategies, it is possible to develop pea varieties that can withstand root rot, ensuring the crop's resilience and its continued importance in global agriculture.

Assessing genome-wide adaptations associated with range expansion in the pink rice borer, Sesamia inferens

Understanding the genetic basis of adaptive evolution following habitat expansion can have important implications for pest management. The pink rice borer (PRB), Sesamia inferens (Walker), is a destructive pest of rice that was historically restricted to regions south of 34° N latitude in China. However, with changes in global climate and farming practices, the distribution of this moth has progressively expanded, encompassing most regions in North China. Here, 3 highly differentiated subpopulations were discovered using high-quality single-nucleotide polymorphism and structural variant datasets across China, corresponding to northern, southern China regions, and the Yunnan-Guizhou Plateau, with significant patterns of isolation by geographic and environmental distances. Our estimates of evolutionary history indicate asymmetric migration with varying population sizes across the 3 subpopulations. Selective sweep analyses estimated strong selection at insect cuticle glycine-rich cuticular protein genes which are associated with enhanced desiccation adaptability in the northern group, and at the histone-lysine-N-methyltransferase gene associated with range expansion and local adaptation in the Shandong population. Our findings have significant implications for the development of effective strategies to control this pest.

A chromosome-scale reference genome of grasspea (Lathyrus sativus)

Grasspea (Lathyrus sativus L.) is an underutilised but promising legume crop with tolerance to a wide range of abiotic and biotic stress factors, and potential for climate-resilient agriculture. Despite a long history and wide geographical distribution of cultivation, only limited breeding resources are available. This paper reports a 5.96 Gbp genome assembly of grasspea genotype LS007, of which 5.03 Gbp is scaffolded into 7 pseudo-chromosomes. The assembly has a BUSCO completeness score of 99.1% and is annotated with 31719 gene models and repeat elements. This represents the most contiguous and accurate assembly of the grasspea genome to date.

Genomic-inferred cross-selection methods for multi-trait improvement in a recurrent selection breeding program

The major drawback to the implementation of genomic selection in a breeding program lies in long-term decrease in additive genetic variance, which is a trade-off for rapid genetic improvement in short term. Balancing increase in genetic gain with retention of additive genetic variance necessitates careful optimization of this trade-off. In this study, we proposed an integrated index selection approach within the genomic inferred cross-selection (GCS) framework to maximize genetic gain across multiple traits. With this method, we identified optimal crosses that simultaneously maximize progeny performance and maintain genetic variance for multiple traits. Using a stochastic simulated recurrent breeding program over a 40-years period, we evaluated different GCS methods along with other factors, such as the number of parents, crosses, and progeny per cross, that influence genetic gain in a pulse crop breeding program. Across all breeding scenarios, the posterior mean variance consistently enhances genetic gain when compared to other methods, such as the usefulness criterion, optimal haploid value, mean genomic estimated breeding value, and mean index selection value of the superior parents. In addition, we provide a detailed strategy to optimize the number of parents, crosses, and progeny per cross that can potentially maximize short- and long-term genetic gain in a public breeding program.

Reference genome sequence and population genomic analysis of peas provide insights into the genetic basis of Mendelian and other agronomic traits

Peas are essential for human nutrition and played a crucial role in the discovery of Mendelian laws of inheritance. In this study, we assembled the genome of the elite vegetable pea cultivar 'Zhewan No. 1' at the chromosome level and analyzed resequencing data from 314 accessions, creating a comprehensive map of genetic variation in peas. We identified 235 candidate loci associated with 57 important agronomic traits through genome-wide association studies. Notably, we pinpointed the causal gene haplotypes responsible for four Mendelian traits: stem length (Le/le), flower color (A/a), cotyledon color (I/i) and seed shape (R/r). Additionally, we discovered the genes controlling pod form (Mendelian P/p) and hilum color. Our study also involved constructing a gene expression atlas across 22 tissues, highlighting key gene modules related to pod and seed development. These findings provide valuable pea genomic information and will facilitate the future genome-informed improvement of pea crops.

Pan-genome and phylogenomic analyses highlight Hevea species delineation and rubber trait evolution

The para rubber tree (Hevea brasiliensis) is the world's sole commercial source of natural rubber, a vital industrial raw material. However, the narrow genetic diversity of this crop poses challenges for rubber breeding. Here, we generate high-quality de novo genome assemblies for three H. brasiliensis cultivars, two H. brasiliensis wild accessions, and three other Hevea species (H. nitida, H. pauciflora, and H. benthamiana). Through analyzing genomes of 94 Hevea accessions, we identify five distinct lineages that do not align with their previous species delineations. We discover multiple accessions with hybrid origins between these lineages, indicating incomplete reproductive isolation between them. Only two out of four wild lineages have been introduced to commercial rubber cultivars. Furthermore, we reveal that the rubber production traits emerged following the development of a large REF/SRPP gene cluster and its functional specialization in rubber-producing laticifers within this genus. These findings would enhance rubber breeding and benefit research communities.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

afila, the origin and nature of a major innovation in the history of pea breeding

The afila (af) mutation causes the replacement of leaflets by a branched mass of tendrils in the compound leaves of pea - Pisum sativum L. This mutation was first described in 1953, and several reports of spontaneous af mutations and induced mutants with a similar phenotype exist. Despite widespread introgression into breeding material, the nature of af and the origin of the alleles used remain unknown. Here, we combine comparative genomics with reverse genetic approaches to elucidate the genetic determinants of af. We also investigate haplotype diversity using a set of AfAf and afaf cultivars and breeding lines and molecular markers linked to seven consecutive genes. Our results show that deletion of two tandemly arranged genes encoding Q-type Cys(2)His(2) zinc finger transcription factors, PsPALM1a and PsPALM1b, is responsible for the af phenotype in pea. Eight haplotypes were identified in the af-harbouring genomic region on chromosome 2. These haplotypes differ in the size of the deletion, covering more or less genes. Diversity at the af locus is valuable for crop improvement and sheds light on the history of pea breeding for improved standing ability. The results will be used to understand the function of PsPALM1a/b and to transfer the knowledge for innovation in related crops.

Genome-Wide Association Studies on Resistance to Pea Weevil: Identification of Novel Sources of Resistance and Associated Markers

Little resistance to the pea weevil insect pest (Bruchus pisorum) is available in pea (Pisum sativum) cultivars, highlighting the need to search for sources of resistance in Pisum germplasm and to decipher the genetic basis of resistance. To address this need, we screened the response to pea weevil in a Pisum germplasm collection (324 accession, previously genotyped) under field conditions over four environments. Significant variation for weevil seed infestation (SI) was identified, with resistance being frequent in P. fulvum, followed by P. sativum ssp. elatius, P. abyssinicum, and P. sativum ssp. humile. SI tended to be higher in accessions with lighter seed color. SI was also affected by environmental factors, being favored by high humidity during flowering and hampered by warm winter temperatures and high evapotranspiration during and after flowering. Merging the phenotypic and genotypic data allowed genome-wide association studies (GWAS) yielding 73 markers significantly associated with SI. Through the GWAS models, 23 candidate genes were found associated with weevil resistance, highlighting the interest of five genes located on chromosome 6. These included gene 127136761 encoding squalene epoxidase; gene 127091639 encoding a transcription factor MYB SRM1; gene 127097033 encoding a 60S ribosomal protein L14; gene 127092211, encoding a BolA-like family protein, which, interestingly, was located within QTL BpLD.I, earlier described as conferring resistance to weevil in pea; and gene 127096593 encoding a methyltransferase. These associated genes offer valuable potential for developing pea varieties resistant to Bruchus spp. and efficient utilization of genomic resources through marker-assisted selection (MAS).

A reference-grade genome of the xerophyte Ammopiptanthus mongolicus sheds light on its evolution history in legumes and drought-tolerance mechanisms

Plants that grow in extreme environments represent unique sources of stress-resistance genes and mechanisms. Ammopiptanthus mongolicus (Leguminosae) is a xerophytic evergreen broadleaf shrub native to semi-arid and desert regions; however, its drought-tolerance mechanisms remain poorly understood. Here, we report the assembly of a reference-grade genome for A. mongolicus, describe its evolutionary history within the legume family, and examine its drought-tolerance mechanisms. The assembled genome is 843.07 Mb in length, with 98.7% of the sequences successfully anchored to the nine chromosomes of A. mongolicus. The genome is predicted to contain 47 611 protein-coding genes, and 70.71% of the genome is composed of repetitive sequences; these are dominated by transposable elements, particularly long-terminal-repeat retrotransposons. Evolutionary analyses revealed two whole-genome duplication (WGD) events at 130 and 58 million years ago (mya) that are shared by the genus Ammopiptanthus and other legumes, but no species-specific WGDs were found within this genus. Ancestral genome reconstruction revealed that the A. mongolicus genome has undergone fewer rearrangements than other genomes in the legume family, confirming its status as a "relict plant". Transcriptomic analyses demonstrated that genes involved in cuticular wax biosynthesis and transport are highly expressed, both under normal conditions and in response to polyethylene glycol-induced dehydration. Significant induction of genes related to ethylene biosynthesis and signaling was also observed in leaves under dehydration stress, suggesting that enhanced ethylene response and formation of thick waxy cuticles are two major mechanisms of drought tolerance in A. mongolicus. Ectopic expression of AmERF2, an ethylene response factor unique to A. mongolicus, can markedly increase the drought tolerance of transgenic Arabidopsis thaliana plants, demonstrating the potential for application of A. mongolicus genes in crop improvement.

Detection of colinear blocks and synteny and evolutionary analyses based on utilization of MCScanX

As different taxa evolve, gene order often changes slowly enough that chromosomal 'blocks' with conserved gene orders (synteny) are discernible. The MCScanX toolkit ( https://github.com/wyp1125/MCScanX ) was published in 2012 as freely available software for the detection of such 'colinear blocks' and subsequent synteny and evolutionary analyses based on genome-wide gene location and protein sequence information. Owing to its simplicity and high efficiency for colinear block detection, MCScanX provides a powerful tool for conducting diverse synteny and evolutionary analyses. Moreover, the detection of colinear blocks has been embraced as an integral step for pangenome graph construction. Here, new application trends of MCScanX are explored, striving to better connect this increasingly used tool to other tools and accelerate insight generation from exponentially growing sequence data. We provide a detailed protocol that covers how to install MCScanX on diverse platforms, tune parameters, prepare input files from data from the National Center for Biotechnology Information, run MCScanX and its visualization and evolutionary analysis tools, and connect MCScanX with external tools, including MCScanX-transposed, Circos and SynVisio. This protocol is easily implemented by users with minimal computational background and is adaptable to new data of interest to them. The data and utility programs for this protocol can be obtained from http://bdx-consulting.com/mcscanx-protocol .

Seed production determines the entrance to dormancy of the inflorescence meristem of Pisum sativum and the end of the flowering period

Flowering plants adjust their reproductive period to maximize the success of the offspring. Monocarpic plants, those with a single reproductive cycle that precedes plant senescence and death, tightly regulate both flowering initiation and flowering cessation. The end of the flowering period involves the arrest of the inflorescence meristem activity, known as proliferative arrest, in what has been interpreted as an evolutionary adaptation to maximize the allocation of resources to seed production and the viability of the progeny. Factors influencing proliferative arrest were described for several monocarpic plant species many decades ago, but only in the last few years studies performed in Arabidopsis have allowed to approach proliferative arrest regulation in a comprehensive manner by studying the physiology, hormone dynamics, and genetic factors involved in its regulation. However, these studies remain restricted to Arabidopsis and there is a need to expand our knowledge to other monocarpic species to propose general mechanisms controlling the process. In this work, we have characterized proliferative arrest in Pisum sativum, trying to parallel available studies in Arabidopsis to maximize this comparative framework. We have assessed quantitatively the role of fruits/seeds in the process, the influence of the positional effect of these fruits/seeds in the behavior of the inflorescence meristem, and the transcriptomic changes in the inflorescence associated with the arrested state of the meristem. Our results support a high conservation of the factors triggering arrest in pea and Arabidopsis, but also reveal differences reinforcing the need to perform similar studies in other species.

Fine mapping and identification of a Fusarium wilt resistance gene FwS1 in pea

KEY MESSAGE

A Fusarium wilt resistance gene FwS1 on pea chromosome 6 was identified and mapped to a 91.4 kb region by a comprehensive genomic-based approach, and the gene Psat6g003960 harboring NB-ARC domain was identified as the putative candidate gene. Pea Fusarium wilt, incited by Fusarium oxysporum f. sp. pisi (Fop), has always been a devastating disease that causes severe yield losses and economic damage in pea-growing regions worldwide. The utilization of pea cultivars carrying resistance gene is the most efficient approach for managing this disease. In order to finely map resistance gene, F2 populations were established through the cross between Shijiadacaiwan 1 (resistant) and Y4 (susceptible). The resistance genetic analysis indicated that the Fop resistance in Shijiadacaiwan 1 was governed by a single dominant gene, named FwS1. Based on the bulked segregant analysis sequencing analyses, the gene FwS1 was initially detected on chromosome 6 (i.e., linking group II, chr6LG2), and subsequent linkage mapping with 589 F2 individuals fine-mapped the gene FwS1 into a 91.4 kb region. The further functional annotation and haplotype analysis confirmed that the gene Psat6g003960, characterized by a NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) domain, was considered as the most promising candidate gene. The encoding amino acids were altered by a "T/C" single-nucleotide polymorphism (SNP) in the first exon of the Psat6g003960, and based on this SNP locus, the molecular marker A016180 was determined to be a diagnostic marker for FwS1 by validating its specificity in both pea accessions and genetic populations with different genetic backgrounds. The FwS1 with diagnostic KASP marker A016180 could facilitate marker-assisted selection in resistance pea breeding in pea. In addition, a comparison of the candidate gene Psat6g003960 in 74SN3B and SJ1 revealed the same sequences. This finding indicated that 74SN3B carried the candidate gene for FwS1, suggesting that FwS1 and Fwf may be closely linked or an identical resistant gene against Fusarium wilt.

Pan-transcriptome reveals a large accessory genome contribution to gene expression variation in yeast

Gene expression is an essential step in the translation of genotypes into phenotypes. However, little is known about the transcriptome architecture and the underlying genetic effects at the species level. Here we generated and analyzed the pan-transcriptome of ~1,000 yeast natural isolates across 4,977 core and 1,468 accessory genes. We found that the accessory genome is an underappreciated driver of transcriptome divergence. Global gene expression patterns combined with population structure showed that variation in heritable expression mainly lies within subpopulation-specific signatures, for which accessory genes are overrepresented. Genome-wide association analyses consistently highlighted that accessory genes are associated with proportionally more variants with larger effect sizes, illustrating the critical role of the accessory genome on the transcriptional landscape within and between populations.

Celebrating Professor Rajeev K. Varshney's transformative research odyssey from genomics to the field on his induction as Fellow of the Royal Society

Professor Rajeev K. Varshney's transformative impact on crop genomics, genetics, and agriculture is the result of his passion, dedication, and unyielding commitment to harnessing the potential of genomics to address the most pressing challenges faced by the global agricultural community. Starting from a small town in India and reaching the global stage, Professor Varshney's academic and professional trajectory has inspired many scientists active in research today. His ground-breaking work, especially his effort to list orphan tropical crops to genomic resource-rich entities, has been transformative. Beyond his scientific achievements, Professor Varshney is recognized by his colleagues as an exemplary mentor, fostering the growth of future researchers, building institutional capacity, and strengthening scientific capability. His focus on translational genomics and strengthening seed system in developing countries for the improvement of agriculture has made a tangible impact on farmers' lives. His skills have been best utilized in roles at leading research centres where he has applied his expertise to deliver a new vision for crop improvement. These efforts have now been recognized by the Royal Society with the award of the Fellowship (FRS). As we mark this significant milestone in his career, we not only celebrate Professor Varshney's accomplishments but also his wider contributions that continue to transform the agricultural landscape.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

Genetic diversity and population structure analysis of a diverse panel of pea (Pisum sativum)

Breeding resilient cultivars with increased tolerance to environmental stress and enhanced resistance to pests and diseases demands pre-breeding efforts that include understanding genetic diversity. This study aimed to evaluate the genetic diversity and population structure of 265 pea accessions. The diversity arrays technology (DArT) genotyping method was employed to identify single-nucleotide polymorphisms (SNPs) and silico markers. After stringent filtering, 6966 SNP and 8,454 silico markers were selected for diversity analysis. Genetic diversity was estimated by grouping accessions based on plant material type, geographic origin, growth habit, and seed color. Generally, diversity estimations obtained using SNPs were similar to those estimated using silico markers. The polymorphism information content (PIC) of the SNP markers ranged from 0.0 to 0.5, with a quarter of them displaying PIC values exceeding 0.4, making them highly informative. Analysis based on plant material type revealed narrow observed heterozygosity (Ho = 0.02-0.03) and expected heterozygosity (He = 0.26-0.31), with landrace accessions exhibiting the highest diversity. Geographic origin-based diversity analysis revealed Ho = 0.02-0.03 and He = 0.22 to 0.30, with European accessions showing the greatest diversity. Moreover, private alleles unique to landrace (4) and European (22) accessions were also identified, which merit further investigation for their potential association with desirable traits. The analysis of molecular variance revealed a highly significant genetic differentiation among accession groups classified by seed color, growth habit, plant material types, and geographic origin (p < 0.01). Principal coordinate analysis and neighbor-joining cluster analysis revealed weak clustering of accessions at different grouping levels. This study underscores the significance of genetic diversity in pea collections, offering valuable insights for targeted breeding and conservation efforts. By leveraging genomic data and exploring untapped genetic resources, pea breeding programs can be fortified to ensure sustainable plant protein production and address future challenges in agriculture.

Mutations of PsPALM1a and PsPALM1b associated with the afila phenotype in Pea

Semi-leafless represents an advantageous plant architecture in pea breeding due to its ability to enhance resistance to lodging and potentially to powdery mildew. The introduction of semi-leafless pea varieties is considered a seminal advancement in pea breeding over the past half-century. The afila (af) mutation leads to the replacement of lateral leaflets by highly branched tendrils; combined with the semi-dwarfing le mutation, it forms the semi-leafless cultivated variety. In this study, we identified that mutations in two tandemly-arrayed genes encoding Cys(2)His(2) zinc finger transcription factors, PsPALM1a and PsPALM1b, were closely associated with the afila phenotype. These two genes may be deleted in the af mutant. In situ hybridization showed that both genes exhibit specific expression in early leaflet primordia. Furthermore, suppression of PsPALM1a/PsPALM1b resulted in a high frequency of conversion of lateral leaflets into tendrils. In conclusion, our study provides genetic evidence demonstrating that mutations in PsPALM1a and PsPALM1b are responsible for the af locus, contributing to a better understanding of compound leaf formation in peas and offering new insights for breeding applications related to afila.

Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

Analysis of pea mutants reveals the conserved role of FRUITFULL controlling the end of flowering and its potential to boost yield

Monocarpic plants have a single reproductive phase in their life. Therefore, flower and fruit production are restricted to the length of this period. This reproductive strategy involves the regulation of flowering cessation by a coordinated arrest of the growth of the inflorescence meristems, optimizing resource allocation to ensure seed filling. Flowering cessation appears to be a regulated phenomenon in all monocarpic plants. Early studies in several species identified seed production as a major factor triggering inflorescence proliferative arrest. Recently, genetic factors controlling inflorescence arrest, in parallel to the putative signals elicited by seed production, have started to be uncovered in Arabidopsis, with the MADS-box gene FRUITFULL (FUL) playing a central role in the process. However, whether the genetic network regulating arrest is also at play in other species is completely unknown. Here, we show that this role of FUL is not restricted to Arabidopsis but is conserved in another monocarpic species with a different inflorescence structure, field pea, strongly suggesting that the network controlling the end of flowering is common to other plants. Moreover, field trials with lines carrying mutations in pea FUL genes show that they could be used to boost crop yield.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

Innovations in functional genomics and molecular breeding of pea: exploring advances and opportunities

The garden pea (Pisum sativum L.) is a significant cool-season legume, serving as crucial food sources, animal feed, and industrial raw materials. The advancement of functional genomics over the past two decades has provided substantial theoretical foundations and progress to pea breeding. Notably, the release of the pea reference genome has enhanced our understanding of plant architecture, symbiotic nitrogen fixation (SNF), flowering time, floral organ development, seed development, and stress resistance. However, a considerable gap remains between pea functional genomics and molecular breeding. This review summarizes the current advancements in pea functional genomics and breeding while highlighting the future challenges in pea molecular breeding.

Population genomics of Agrotis segetum provide insights into the local adaptive evolution of agricultural pests

BACKGROUND

The adaptive mechanisms of agricultural pests are the key to understanding the evolution of the pests and to developing new control strategies. However, there are few studies on the genetic basis of adaptations of agricultural pests. The turnip moth, Agrotis segetum (Lepidoptera: Noctuidae) is an important underground pest that affects a wide range of host plants and has a strong capacity to adapt to new environments. It is thus a good model for studying the adaptive evolution of pest species.

RESULTS

We assembled a high-quality reference genome of A. segetum using PacBio reads. Then, we constructed a variation map of A. segetum by resequencing 98 individuals collected from six natural populations in China. The analysis of the population structure showed that all individuals were divided into four well-differentiated populations, corresponding to their geographical distribution. Selective sweep analysis and environmental association studies showed that candidate genes associated with local adaptation were functionally correlated with detoxification metabolism and glucose metabolism.

CONCLUSIONS

Our study of A. segetum has provided insights into the genetic mechanisms of local adaptation and evolution; it has also produced genetic resources for developing new pest management strategies.

Genome-wide identification, structural characterization and gene expression analysis of the WRKY transcription factor family in pea (Pisum sativum L.)

BACKGROUND

The WRKY gene family is one of the largest families of transcription factors in higher plants, and WRKY transcription factors play important roles in plant growth and development as well as in response to abiotic stresses; however, the WRKY gene family in pea has not been systematically reported.

RESULTS

In this study, 89 pea WRKY genes were identified and named according to the random distribution of PsWRKY genes on seven chromosomes. The gene family was found to have nine pairs of tandem duplicates and 19 pairs of segment duplicates. Phylogenetic analyses of the PsWRKY and 60 Arabidopsis WRKY proteins were performed to determine their homology, and the PsWRKYs were classified into seven subfamilies. Analysis of the physicochemical properties, motif composition, and gene structure of pea WRKYs revealed significant differences in the physicochemical properties within the PsWRKY family; however, their gene structure and protein-conserved motifs were highly conserved among the subfamilies. To further investigate the evolutionary relationships of the PsWRKY family, we constructed comparative syntenic maps of pea with representative monocotyledonous and dicotyledonous plants and found that it was most recently homologous to the dicotyledonous WRKY gene families. Cis-acting element analysis of PsWRKY genes revealed that this gene family can respond to hormones, such as abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin (GA), methyl jasmonate (MeJA), and salicylic acid (SA). Further analysis of the expression of 14 PsWRKY genes from different subfamilies in different tissues and fruit developmental stages, as well as under five different hormone treatments, revealed differences in their expression patterns in the different tissues and fruit developmental stages, as well as under hormone treatments, suggesting that PsWRKY genes may have different physiological functions and respond to hormones.

CONCLUSIONS

In this study, we systematically identified WRKY genes in pea for the first time and further investigated their physicochemical properties, evolution, and expression patterns, providing a theoretical basis for future studies on the functional characterization of pea WRKY genes during plant growth and development.

Fine mapping of Ae-Ps4.5, a major locus for resistance to pathotype III of Aphanomyces euteiches in pea

KEY MESSAGE

QTL mapping and recombinant screening confirmed the major effect of QTL Ae-Ps4.5 on pea resistance to pathotype III of Aphanomyces euteiches and fine-mapped the QTL to a 3.06-Mb interval. Aphanomyces root rot, caused by Aphanomyces euteiches, is the most important disease of pea (Pisum sativum L.) worldwide. The development of pea-resistant varieties is a major challenge to control the disease. Previous linkage studies identified seven main resistance quantitative trait loci (QTL), including the QTL Ae-Ps4.5 associated with partial resistance in US nurseries infested by the pea pathotype III of A. euteiches. This study aimed to confirm the major effect of Ae-Ps4.5 on A. euteiches pathotype III, refine its interval, and identify candidate genes underlying the QTL. QTL mapping on an updated genetic map from the Puget × 90-2079 pea recombinant inbred line population identified Ae-Ps4.5 in a 0.8-cM confidence interval with a high effect (R2 = 89%) for resistance to the Ae109 reference strain of A. euteiches (pathotype III) under controlled conditions. However, the QTL mapping did not detect Ae-Ps4.5 for resistance to the RB84 reference strain of A. euteiches (pathotype I). Screening 224-pea BC5F2 plant progeny derived from three near-isogenic lines (NILs) carrying the 90-2079 allele at Ae-Ps4.5 in the Puget genetic background with 26 SNP markers identified 15 NILs showing recombination in the QTL interval. Phenotyping of the recombinant lines for resistance to the Ae109 strain of A. euteiches reduced the QTL to a physical interval of 3.06 Mb, containing 50 putative annotated genes on the Caméor pea genome V1a among which three candidate genes highlighted. This study provides closely linked SNP markers and putative candidate genes to accelerate pea breeding for resistant varieties to Aphanomyces root rot.

The genomes of 5 underutilized Papilionoideae crops provide insights into root nodulation and disease resistance

BACKGROUND

The Papilionoideae subfamily contains a large amount of underutilized legume crops, which are important for food security and human sustainability. However, the lack of genomic resources has hindered the breeding and utilization of these crops.

RESULTS

Here, we present chromosome-level reference genomes for 5 underutilized diploid Papilionoideae crops: sword bean (Canavalia gladiata), scarlet runner bean (Phaseolus coccineus), winged bean (Psophocarpus tetragonolobus), smooth rattlebox (Crotalaria pallida), and butterfly pea (Clitoria ternatea), with assembled genome sizes of 0.62 Gb, 0.59 Gb, 0.71 Gb, 1.22 Gb, and 1.72 Gb, respectively. We found that the long period of higher long terminal repeat retrotransposon activity is the major reason that the genome size of smooth rattlebox and butterfly pea is enlarged. Additionally, there have been no recent whole-genome duplication (WGD) events in these 5 species except for the shared papilionoid-specific WGD event (∼55 million years ago). Then, we identified 5,328 and 10,434 species-specific genes between scarlet runner bean and common bean, respectively, which may be responsible for their phenotypic and functional differences and species-specific functions. Furthermore, we identified the key genes involved in root-nodule symbiosis (RNS) in all 5 species and found that the NIN gene was duplicated in the early Papilionoideae ancestor, followed by the loss of 1 gene copy in smooth rattlebox and butterfly pea lineages. Last, we identified the resistance (R) genes for plant defenses in these 5 species and characterized their evolutionary history.

CONCLUSIONS

In summary, this study provides chromosome-scale reference genomes for 3 grain and vegetable beans (sword bean, scarlet runner bean, winged bean), along with genomes for a green manure crop (smooth rattlebox) and a food dyeing crop (butterfly pea). These genomes are crucial for studying phylogenetic history, unraveling nitrogen-fixing RNS evolution, and advancing plant defense research.

Telomere-to-telomere genome of the allotetraploid legume Sesbania cannabina reveals transposon-driven subgenome divergence and mechanisms of alkaline stress tolerance

Alkaline soils pose an increasing problem for agriculture worldwide, but using stress-tolerant plants as green manure can improve marginal land. Here, we show that the legume Sesbania cannabina is very tolerant to alkaline conditions and, when used as a green manure, substantially improves alkaline soil. To understand genome evolution and the mechanisms of stress tolerance in this allotetraploid legume, we generated the first telomere-to-telomere genome assembly of S. cannabina spanning ∼2,087 Mb. The assembly included all centromeric regions, which contain centromeric satellite repeats, and complete chromosome ends with telomeric characteristics. Further genome analysis distinguished A and B subgenomes, which diverged approximately 7.9 million years ago. Comparative genomic analysis revealed that the chromosome homoeologs underwent large-scale inversion events (>10 Mb) and a significant, transposon-driven size expansion of the chromosome 5A homoeolog. We further identified four specific alkali-induced phosphate transporter genes in S. cannabina; these may function in alkali tolerance by relieving the deficiency in available phosphorus in alkaline soil. Our work highlights the significance of S. cannabina as a green tool to improve marginal lands and sheds light on subgenome evolution and adaptation to alkaline soils.

Photosynthesis governed by nanoparticulate titanium dioxide. The Pisum sativum L. case study

Wide availability of anthropogenic TiO2 nanoparticles facilitates their penetration into environment and prompts interactions with plants. They alter plants growth and change their nutritional status. In particular, metabolic processes are affected. In this work the effect of nanometric TiO2 on photosynthesis efficiency in green pea (Pisum sativum L.) was studied. Hydroponic cultivations with three Ti levels (10; 50 and 100 mg L-1) were applied. At all concentrations nanoparticles penetrated into plant tissues and were detected by the single particle ICP-MS/MS method. Nanoparticles altered the CO2 assimilation rate and gas exchange parameters (i.e. transpiration, stomatal conductance, sub-stomatal CO2 concentration). The most pronounced effects were observed for Ti 50 mg L-1 cultivation where photosynthesis efficiency, transpiration and stomatal conductance were increased by 14.69%, 4.58% and 8.92%, respectively. They were further confirmed by high maximum ribulose 1,5-bisphosphate carboxylation rate (27.40% increase), maximum electron transport rate (21.51% increase) and the lowest CO2 compensation point (45.19% decrease). Furthermore, concentrations of Cu, Mn, Zn, Fe, Mg, Ca, K and P were examined with the most pronounced changes observed for elements directly involved in photosynthesis (Cu, Zn, Mn, and Fe). The Cu concentrations in roots, stems and leaves for Ti 50 mg L-1 cultivation were below the control by 33.15%, 38.28% and 10.76%, respectively. The Zn content in analogous treatment and organs decreased by 30.24%, 26.69% and 13.35%. The Mn and Fe levels in leaves were increased by 72.22% and 50.32%, respectively. Our results indicated that plant defence mechanisms which restrain the water uptake have been overcome in pea by photocatalytic activity of nanoparticulate TiO2 which stimulated photosynthesis. On the contrary to the substantial stomatal conductance, the transpiration has been reduced because exceptional part of water flow was already consumed in chloroplasts and could not have been freed to the atmosphere.

Recombinant inbred lines derived from wide crosses in Pisum

Genomic resources are becoming available for Pisum but to link these to phenotypic diversity requires well marked populations segregating for relevant traits. Here we describe two such resources. Two recombinant inbred populations, derived from wide crosses in Pisum are described. One high resolution mapping population involves cv Caméor, for which the first pea whole genome assembly was obtained, crossed to JI0281, a basally divergent P. sativum sativum landrace from Ethiopia. The other is an inter sub-specific cross between P. s. sativum and the independently domesticated P. s. abyssinicum. The corresponding genetic maps provide information on chromosome level sequence assemblies and identify structural differences between the genomes of these two Pisum subspecies. In order to visualise chromosomal translocations that distinguish the mapping parents, we created a simplified version of Threadmapper to optimise it for interactive 3-dimensional display of multiple linkage groups. The genetic mapping of traits affecting seed coat roughness and colour, plant height, axil ring pigmentation, leaflet number and leaflet indentation enabled the definition of their corresponding genomic regions. The consequence of structural rearrangement for trait analysis is illustrated by leaf serration. These analyses pave the way for identification of the underlying genes and illustrate the utility of these publicly available resources. Segregating inbred populations derived from wide crosses in Pisum, together with the associated marker data, are made publicly available for trait dissection. Genetic analysis of these populations is informative about chromosome scale assemblies, structural diversity in the pea genome and has been useful for the fine mapping of several discrete and quantitative traits.

Representing true plant genomes: haplotype-resolved hybrid pepper genome with trio-binning

As sequencing costs decrease and availability of high fidelity long-read sequencing increases, generating experiment specific de novo genome assemblies becomes feasible. In many crop species, obtaining the genome of a hybrid or heterozygous individual is necessary for systems that do not tolerate inbreeding or for investigating important biological questions, such as hybrid vigor. However, most genome assembly methods that have been used in plants result in a merged single sequence representation that is not a true biologically accurate representation of either haplotype within a diploid individual. The resulting genome assembly is often fragmented and exhibits a mosaic of the two haplotypes, referred to as haplotype-switching. Important haplotype level information, such as causal mutations and structural variation is therefore lost causing difficulties in interpreting downstream analyses. To overcome this challenge, we have applied a method developed for animal genome assembly called trio-binning to an intra-specific hybrid of chili pepper (Capsicum annuum L. cv. HDA149 x Capsicum annuum L. cv. HDA330). We tested all currently available softwares for performing trio-binning, combined with multiple scaffolding technologies including Bionano to determine the optimal method of producing the best haplotype-resolved assembly. Ultimately, we produced highly contiguous biologically true haplotype-resolved genome assemblies for each parent, with scaffold N50s of 266.0 Mb and 281.3 Mb, with 99.6% and 99.8% positioned into chromosomes respectively. The assemblies captured 3.10 Gb and 3.12 Gb of the estimated 3.5 Gb chili pepper genome size. These assemblies represent the complete genome structure of the intraspecific hybrid, as well as the two parental genomes, and show measurable improvements over the currently available reference genomes. Our manuscript provides a valuable guide on how to apply trio-binning to other plant genomes.

Improved high-quality reference genome of red drum facilitates the processes of resistance-related gene exploration

Sciaenops ocellatus is among the most important artificially introduced farmed fish across 11 countries and regions. However, the frequent occurrence of extreme weather events and breeding escapes have placed great pressure on local marine biodiversity and ecosystems. We reported the de novo assembly and annotation with a contig N50 of 28.30 Mb using PacBio HiFi sequencing and Hi-C technologies, which resulted in a 283-fold increase in contig N50 length and improvement in continuity and quality in complex repetitive region for S. ocellatus compared to the previous version. In total, 257.36 Mb of repetitive sequences accounted for 35.48% of the genome, and 22,845 protein-coding genes associated with a BUSCO value of 98.32%, were identified by genome annotation. Moreover, 54 hub genes rapidly responding to hypoosmotic stress were identified by WGCNA. The high-quality chromosome-scale S. ocellatus genome and candidate resistance-related gene sets will not only provide a genomic basis for genetic improvement via molecular breeding, but will also lay an important foundation for investigating the molecular regulation of rapid responses to stress.

PanKmer: k-mer-based and reference-free pangenome analysis

SUMMARY

Pangenomes are replacing single reference genomes as the definitive representation of DNA sequence within a species or clade. Pangenome analysis predominantly leverages graph-based methods that require computationally intensive multiple genome alignments, do not scale to highly complex eukaryotic genomes, limit their scope to identifying structural variants (SVs), or incur bias by relying on a reference genome. Here, we present PanKmer, a toolkit designed for reference-free analysis of pangenome datasets consisting of dozens to thousands of individual genomes. PanKmer decomposes a set of input genomes into a table of observed k-mers and their presence-absence values in each genome. These are stored in an efficient k-mer index data format that encodes SNPs, INDELs, and SVs. It also includes functions for downstream analysis of the k-mer index, such as calculating sequence similarity statistics between individuals at whole-genome or local scales. For example, k-mers can be "anchored" in any individual genome to quantify sequence variability or conservation at a specific locus. This facilitates workflows with various biological applications, e.g. identifying cases of hybridization between plant species. PanKmer provides researchers with a valuable and convenient means to explore the full scope of genetic variation in a population, without reference bias.

AVAILABILITY AND IMPLEMENTATION

PanKmer is implemented as a Python package with components written in Rust, released under a BSD license. The source code is available from the Python Package Index (PyPI) at https://pypi.org/project/pankmer/ as well as Gitlab at https://gitlab.com/salk-tm/pankmer. Full documentation is available at https://salk-tm.gitlab.io/pankmer/.

Telomere-to-telomere genome assembly of melon (Cucumis melo L. var. inodorus) provides a high-quality reference for meta-QTL analysis of important traits

Melon is an important horticultural crop with extensive diversity in many horticultural groups. To explore its genomic diversity, it is necessary to assemble more high-quality complete genomes from different melon accessions. Meanwhile, a large number of QTLs have been mapped in several studies. Integration of the published QTLs onto a complete genome can provide more accurate information for candidate gene cloning. To address these problems, a telomere-to-telomere (T2T) genome of the elite melon landrace Kuizilikjiz (Cucumis melo L. var. inodorus) was de novo assembled and all the published QTLs were projected onto it in this study. The results showed that a high-quality Kuizilikjiz genome with the size of 379.2 Mb and N50 of 31.7 Mb was de novo assembled using the combination of short reads, PacBio high-fidelity long reads, Hi-C data, and a high-density genetic map. Each chromosome contained the centromere and telomeres at both ends. A large number of structural variations were observed between Kuizilikjiz and the other published genomes. A total of 1294 QTLs published in 67 studies were collected and projected onto the T2T genome. Several clustered, co-localized, and overlapped QTLs were determined. Furthermore, 20 stable meta-QTLs were identified, which significantly reduced the mapping intervals of the initial QTLs and greatly facilitated identification of the candidate genes. Collectively, the T2T genome assembly together with the numerous projected QTLs will not only broaden the high-quality genome resources but also provide valuable and abundant QTL information for cloning the genes controlling important traits in melon.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

RGB image-based method for phenotyping rust disease progress in pea leaves using R

BACKGROUND

Rust is a damaging disease affecting vital crops, including pea, and identifying highly resistant genotypes remains a challenge. Accurate measurement of infection levels in large germplasm collections is crucial for finding new resistance sources. Current evaluation methods rely on visual estimation of disease severity and infection type under field or controlled conditions. While they identify some resistance sources, they are error-prone and time-consuming. An image analysis system proves useful, providing an easy-to-use and affordable way to quickly count and measure rust-induced pustules on pea samples. This study aimed to develop an automated image analysis pipeline for accurately calculating rust disease progression parameters under controlled conditions, ensuring reliable data collection.

RESULTS

A highly efficient and automatic image-based method for assessing rust disease in pea leaves was developed using R. The method's optimization and validation involved testing different segmentation indices and image resolutions on 600 pea leaflets with rust symptoms. The approach allows automatic estimation of parameters like pustule number, pustule size, leaf area, and percentage of pustule coverage. It reconstructs time series data for each leaf and integrates daily estimates into disease progression parameters, including latency period and area under the disease progression curve. Significant variation in disease responses was observed between genotypes using both visual ratings and image-based analysis. Among assessed segmentation indices, the Normalized Green Red Difference Index (NGRDI) proved fastest, analysing 600 leaflets at 60% resolution in 62 s with parallel processing. Lin's concordance correlation coefficient between image-based and visual pustule counting showed over 0.98 accuracy at full resolution. While lower resolution slightly reduced accuracy, differences were statistically insignificant for most disease progression parameters, significantly reducing processing time and storage space. NGRDI was optimal at all time points, providing highly accurate estimations with minimal accumulated error.

CONCLUSIONS

A new image-based method for monitoring pea rust disease in detached leaves, using RGB spectral indices segmentation and pixel value thresholding, improves resolution and precision. It rapidly analyses hundreds of images with accuracy comparable to visual methods and higher than other image-based approaches. This method evaluates rust progression in pea, eliminating rater-induced errors from traditional methods. Implementing this approach to evaluate large germplasm collections will improve our understanding of plant-pathogen interactions and aid future breeding for novel pea cultivars with increased rust resistance.

Application of Multi-Perspectives in Tea Breeding and the Main Directions

Tea plants are an economically important crop and conducting research on tea breeding contributes to enhancing the yield and quality of tea leaves as well as breeding traits that satisfy the requirements of the public. This study reviews the current status of tea plants germplasm resources and their utilization, which has provided genetic material for the application of multi-omics, including genomics and transcriptomics in breeding. Various molecular markers for breeding were designed based on multi-omics, and available approaches in the direction of high yield, quality and resistance in tea plants breeding are proposed. Additionally, future breeding of tea plants based on single-cellomics, pangenomics, plant-microbe interactions and epigenetics are proposed and provided as references. This study aims to provide inspiration and guidance for advancing the development of genetic breeding in tea plants, as well as providing implications for breeding research in other crops.

Development of a knowledge graph framework to ease and empower translational approaches in plant research: a use-case on grain legumes

While the continuing decline in genotyping and sequencing costs has largely benefited plant research, some key species for meeting the challenges of agriculture remain mostly understudied. As a result, heterogeneous datasets for different traits are available for a significant number of these species. As gene structures and functions are to some extent conserved through evolution, comparative genomics can be used to transfer available knowledge from one species to another. However, such a translational research approach is complex due to the multiplicity of data sources and the non-harmonized description of the data. Here, we provide two pipelines, referred to as structural and functional pipelines, to create a framework for a NoSQL graph-database (Neo4j) to integrate and query heterogeneous data from multiple species. We call this framework Orthology-driven knowledge base framework for translational research (Ortho_KB). The structural pipeline builds bridges across species based on orthology. The functional pipeline integrates biological information, including QTL, and RNA-sequencing datasets, and uses the backbone from the structural pipeline to connect orthologs in the database. Queries can be written using the Neo4j Cypher language and can, for instance, lead to identify genes controlling a common trait across species. To explore the possibilities offered by such a framework, we populated Ortho_KB to obtain OrthoLegKB, an instance dedicated to legumes. The proposed model was evaluated by studying the conservation of a flowering-promoting gene. Through a series of queries, we have demonstrated that our knowledge graph base provides an intuitive and powerful platform to support research and development programmes.

The wall-associated kinase gene family in pea (Pisum sativum) and its function in response to B deficiency and Al toxicity

Plant cell walls are embedded in a pectin matrix which is physically linked with the wall-associated kinases (WAKs), a subfamily of receptor-like kinases that participate in the cell wall integrity (CWI) sensing. Since cell walls are also the main binding sites for boron (B) and aluminum (Al), WAK may be potentially associated with the regulation of plant responses to Al toxicity and B deficiency. Using pea as a model species, we have identified a total of 28 WAK genes in the genome and named them according to its chromosomal location. All the PsWAKs were phylogenetically grouped into three clades. Phylogenetic relationship and synteny analysis showed that the PsWAKs in pea and Glycine max or Medicago truncatula shared a relatively conserved evolutionary history. Protein domain, motif, and transmembrane analysis indicated that all PsWAK proteins were predicted to be localized to the plasma membrane, and most PsWAKs shared a similar structure to their homologs. The RNA-seq data showed that the expression pattern of WAK genes in response to B deficiency was similar to that of Al toxicity, with most of PsWAKs being up-regulated. The qRT-PCR results further confirmed that PsWAK5, PsWAK9 and PsWAK14 were more specific for both B-deficiency and Al toxicity, and the expression levels of PsWAK5, PsWAK9 and PsWAK14 were significantly higher in the Al-sensitive cultivar Hyogo than in the Al-resistant cultivar Alaska under Al toxicity. This study provided an important basis for the functional and evolutionary analysis of PsWAKs and linked them to responses to cell wall damage induced by B-deficiency and Al toxicity, suggesting that PsWAKs may play a key role in the perception of cell wall integrity under Al toxicity or B-deficiency, as well as in the regulation of Al tolerance in pea.

High-density linkage mapping and genetic dissection of resistance to broomrape (Orobanche crenata Forsk.) in pea (Pisum sativum L.)

Pea (Pisum sativum L.) is a widely cultivated legume of major importance for global food security and agricultural sustainability. Crenate broomrape (Orobanche crenata Forsk.) (Oc) is a parasitic weed severely affecting legumes, including pea, in the Mediterranean Basin and the Middle East. Previously, the identification of the pea line "ROR12", displaying resistance to Oc, was reported. Two-year field trials on a segregant population of 148 F₇ recombinant inbred lines (RILs), originating from a cross between "ROR12" and the susceptible cultivar "Sprinter", revealed high heritability (0.84) of the "ROR12" resistance source. Genotyping-by-sequencing (GBS) on the same RIL population allowed the construction of a high-density pea linkage map, which was compared with the pea reference genome and used for quantitative trait locus (QTL) mapping. Three QTLs associated with the response to Oc infection, named PsOcr-1, PsOcr-2, and PsOcr-3, were identified, with PsOcr-1 explaining 69.3% of the genotypic variance. Evaluation of the effects of different genotypic combinations indicated additivity between PsOcr-1 and PsOcr-2, and between PsOcr-1 and PsOcr-3, and epistasis between PsOcr-2 and PsOcr-3. Finally, three Kompetitive Allele Specific PCR (KASP) marker assays were designed on the single-nucleotide polymorphisms (SNPs) associated with the QTL significance peaks. Besides contributing to the development of pea genomic resources, this work lays the foundation for the obtainment of pea cultivars resistant to Oc and the identification of genes involved in resistance to parasitic Orobanchaceae.