Graph-based pan-genome reveals structural and sequence variations related to agronomic traits and domestication in cucumber

Reference-free identification and pangenome analysis of accessory chromosomes in a major fungal plant pathogen

Accessory chromosomes, found in some but not all individuals of a species, play an important role in pathogenicity and host specificity in fungal plant pathogens. However, their variability complicates reference-based analysis, especially when these chromosomes are missing in the reference genome. Pangenome variation graphs offer a reference-free alternative for studying these chromosomes. Here, we constructed a pangenome variation graph for 73 diverse Fusarium oxysporum genomes, a major fungal plant pathogen with a compartmentalized genome that includes conserved core as well as variable accessory chromosomes. To obtain insights into accessory chromosome dynamics, we first constructed a chromosome similarity network using all-vs-all similarity mapping. We identified eleven core chromosomes conserved across all strains and a substantial number of highly variable accessory chromosomes. Some of these accessory chromosomes are host-specific and likely play a role in determining host range. Using a k-mer based approach, we further identified the presence of these accessory chromosomes in all available (581) F. oxysporum assemblies and corroborated the occurrence of host-specific accessory chromosomes. To further analyze the evolution of chromosomes in F. oxysporum, we constructed a pangenome variation graph per group of homologous chromosomes. This reveals that accessory chromosomes are composed of different stretches of accessory regions, and possibly rearrangements between accessory regions gave rise to these mosaic accessory chromosomes. Furthermore, we show that accessory chromosomes are likely horizontally transferred in natural populations. Our findings demonstrate that a pangenome variation graph is a powerful approach to elucidate the evolutionary dynamics of accessory chromosomes in F. oxysporum, which is not only a useful resource for Fusarium but also provides a framework for similar analyses in other species containing accessory chromosomes.

Identification of regulatory loci and candidate genes related to body weight traits in broilers based on different models

BACKGROUND

Growth traits are crucial for the economic viability in broiler production, as they significantly contribute to the cost of rearing. Maximizing body weight (BW) while minimizing feed intake is key to enhancing the efficiency of broiler breeding. Identifying the genetic architecture associated with BW trait is therefore a critical step in enhancing breeding strategies.

RESULTS

We conducted a genome-wide association study (GWAS) using two statistical approaches: single-trait GWAS and longitudinal GWAS. The study was performed on the BW trait at five developmental stages (72, 81, 89, 113, and 120 days) and mid-test metabolic weight (MWT) across four growth cycles. Transcriptome sequencing analysis was also included to investigate the differential expression of candidate genes identified through the GWAS models, particularly linked to BW and MWT traits. Utilizing the chicken 55K single nucleotide polymorphism (SNP) array, we identified 52,060 SNPs in the genomic data of 4,493 Wenchang chickens. The single-trait GWAS model revealed 42 BW-associated SNPs, corresponding to 18 potential genes. For MWT, 47 SNPs were associated, mapping to 31 candidate genes. The longitudinal GWAS model identified 34 BW-linked SNPs, annotated with 22 candidate genes, and 21 MWT-linked SNPs, annotated with 10 candidate genes. Notably, 16 SNPs on chromosome 4 were associated with both BW and MWT, located within the 73.08Mb-76.82Mb region. Nine genes were annotated from this region, including STIM2, SEL1L3, SEPSECS, LGI2, SOD3, KCNIP4, NCAPG, FAM184B, LDB2. Notably, there are 32 overlapping SNPs identified in both the single-trait and longitudinal GWAS models, suggesting consistent associations for both BW and MWT. These overlapping SNPs represent robust loci that may influence both traits across different statistical approaches. Transcriptome sequencing indicated differential expression of LDB2 and SEL1L3 between high and low BW groups.

CONCLUSION

Our study has uncovered novel candidate genes that are potentially involved in growth traits, providing valuable insights for broiler breeding. The identified SNPs and genes could serve as genetic markers for selecting broilers with improved growth efficiency, which may lead to more cost-effective and productive broiler farming.

Understanding the epidemiology and pathogenesis of Mycobacterium tuberculosis with non-redundant pangenome of epidemic strains in China

Tuberculosis is a major public health threat resulting in more than one million lives lost every year. Many challenges exist to defeat this deadly infectious disease which address the importance of a thorough understanding of the biology of the causative agent Mycobacterium tuberculosis (MTB). We generated a non-redundant pangenome of 420 epidemic MTB strains from China including 344 Lineage 2 strains, 69 Lineage 4 strains, six Lineage 3 strains, and one Lineage 1 strain. We estimate that MTB strains have a pangenome of 4,278 genes encoding 4,183 proteins, of which 3,438 are core genes. However, due to 99,694 interruptions in 2,447 coding genes, we can only confidently confirm 1,651 of these genes are translated in all samples. Of these interruptions, 67,315 (67.52%) could be classified by various genetic variations detected by currently available tools, and more than half of them are due to structural variations, mostly small indels. Assuming a proportion of these interruptions are artifacts, the number of active core genes would still be much lower than 3,438. We further described differential evolutionary patterns of genes under the influences of selective pressure, population structure and purifying selection. While selective pressure is ubiquitous among these coding genes, evolutionary adaptations are concentrated in 1,310 genes. Genes involved in cell wall biogenesis are under the strongest selective pressure, while the biological process of disruption of host organelles indicates the direction of the most intensive positive selection. This study provides a comprehensive view on the genetic diversity and evolutionary patterns of coding genes in MTB which may deepen our understanding of its epidemiology and pathogenicity.

Lake Malawi cichlid pangenome graph reveals extensive structural variation driven by transposable elements

Pangenome methods have the potential to uncover hitherto undiscovered sequences missing from established reference genomes, making them useful to study evolutionary and speciation processes in diverse organisms. The cichlid fishes of the East African Rift Lakes represent one of nature's most phenotypically diverse vertebrate radiations, but single-nucleotide polymorphism (SNP)-based studies have revealed little sequence difference, with 0.1%-0.25% pairwise divergence between Lake Malawi species. These were based on aligning short reads to a single linear reference genome and ignored the contribution of larger-scale structural variants (SVs). We constructed a pangenome graph that integrates six new and two existing long-read genome assemblies of Lake Malawi haplochromine cichlids. This graph intuitively represents complex and nested variation between the genomes and reveals that the SV landscape is dominated by large insertions, many exclusive to individual assemblies. The graph incorporates a substantial amount of extra sequence across seven species, the total size of which is 33.1% longer than that of a single cichlid genome. Approximately 4.73% to 9.86% of the assembly lengths are estimated as interspecies structural variation between cichlids, suggesting substantial genomic diversity underappreciated in SNP studies. Although coding regions remain highly conserved, our analysis uncovers a significant proportion of SV sequences as transposable element (TE) insertions, especially DNA, LINE, and LTR TEs. These findings underscore that the cichlid genome is shaped both by small-nucleotide mutations and large, TE-derived sequence alterations, both of which merit study to understand their interplay in cichlid evolution.

What Are We Learning from Plant Pangenomes?

A single reference genome does not fully capture species diversity. By contrast, a pangenome incorporates multiple genomes to capture the entire set of nonredundant genes in a given species, along with its genome diversity. New sequencing technologies enable researchers to produce multiple high-quality genome sequences and catalog diverse genetic variations with better precision. Pangenomic studies have detected structural variants in plant genomes, dissected the genetic architecture of agronomic traits, and helped unravel molecular underpinnings and evolutionary origins of plant phenotypes. The pangenome concept has further evolved into a so-called super-pangenome that includes wild relatives within a genus or clade and shifted to graph-based reference systems. Nevertheless, building pangenomes and representing complex structural variants remain challenging in many crops. Standardized computing pipelines and common data structures are needed to compare and interpret pangenomes. The growing body of plant pangenomics data requires new algorithms, huge data storage capacity, and training to help researchers and breeders take advantage of newly discovered genes and genetic variants.

Research Advances and Perspectives on Early Flowering Traits in Cucumber

Early flowering refers to the phenomenon in which the first flower appears in fewer days than normal, regardless of the sex of the flower. It is a significant feature impacting the early maturity and economic yield of cucumbers. The early flowering trait of cucumber is influenced by several factors. Considering its heritability, technologies such as whole-genome sequencing, genetic modification, bioinformatics analysis, quantitative trait locus (QTL) mapping, molecular marker-assisted selection, and gene editing are widely used to explore the regulatory genes and molecular mechanisms of the early flowering trait in cucumbers. This review aimed to summarize the factors, QTL mapping, molecular regulation mechanisms, and omics analysis related to early flowering traits in cucumbers. This review contributes theoretical insights to support both cucumber breeding for early flowering and fundamental research on early flowering traits.

Structural variation-based and gene-based pangenome construction reveals untapped diversity of hexaploid wheat

Increasing number of structural variations (SVs) have been identified as causative mutations for diverse agronomic traits. However, the systematic exploration of SVs quantity, distribution, and contribution in wheat was lacking. Here, we report high-quality gene-based and SV-based pangenomes comprising 22 hexaploid wheat assemblies showing a wide range of chromosome size, gene number, and TE component, which indicates their representativeness of wheat genetic diversity. Pan-gene analyses uncover 140,261 distinct gene families, of which only 23.2 % are shared in all accessions. Moreover, we build a ∼16.15 Gb graph pangenome containing 695,897 bubbles, intersecting 5132 genes and 230,307 cis-regulatory regions. Pairwise genome comparisons identify ∼1,978,221 non-redundant SVs and 497 SV hotspots. Notably, the density of bubbles as well as SVs show remarkable aggregation in centromeres, which probably play an important role in chromosome plasticity and stability. As for functional SVs exploration, we identify 2769 SVs with absolute relative frequency differences exceeding 0.7 between spring and winter growth habit groups. Additionally, several reported functional genes in wheat display complex structural graphs, for example, PPD-A1, VRT-A2, and TaNAAT2-A. These findings deepen our understanding of wheat genetic diversity, providing valuable graphical pangenome and variation resources to improve the efficiency of genome-wide association mapping in wheat.

Graphical pangenomics-enabled characterization of structural variant impact on gene expression in Brassica napus

KEY MESSAGE

Pangenome graphs enable population-scale genotyping and improve expression analysis, revealing that structural variations (SVs), particularly transposable elements (TEs), significantly contribute to gene expression variation in winter oilseed rape. Structural variations (SVs) impact important traits, from yield to flowering behaviour and stress responses. Pangenome graphs capture population-level diversity, including SVs, within a single data structure and provide a robust framework for downstream applications. They have the potential to serve as unbiased references for SV genotyping, pan-transcriptomic analyses, and association studies, offering significant advantages over single reference genomes. However, their full potential for expression quantitative trait locus (eQTL) analysis is yet to be explored. We combined long and short-read whole genome sequencing data with expression profiling of Brassica napus (oilseed rape) to assess the impact of SVs on gene expression regulation and explored the utility of pangenome graphs for eQTL analysis. Over 90,000 SVs were discovered from 57 long-read datasets. Pangenome graph as reference was evaluated and used for SV genotyping with short reads and transcript expression quantification. Using SVs genotyped from the graph and 100 expression datasets, we identified 267 gene proximal (cis) SV-eQTLs. Over 70% of eQTL-SVs had similarity to transposable elements (TEs), especially Helitrons. The highest proportion of cis-eQTL-SVs were found in promoter regions. About a third of transcripts whose expression was associated with SVs, had no associated SNPs, suggesting that including SVs allows capturing of relationship which would be missed in SNP-only analyses. This study demonstrated that pangenome graphs provide a unifying framework for eQTL analysis by allowing population-scale SV genotyping and gene expression quantification. We also showed that SVs make an appreciable contribution to gene expression variation in winter oilseed rape.

Gastric cancer genomics study using reference human pangenomes

A pangenome is the sum of the genetic information of all individuals in a species or a population. Genomics research has been gradually shifted to a paradigm using a pangenome as the reference. However, in disease genomics study, pangenome-based analysis is still in its infancy. In this study, we introduced a graph-based pangenome GGCPan from 185 patients with gastric cancer. We then systematically compared the cancer genomics study results using GGCPan, a linear pangenome GCPan, and the human reference genome as the reference. For small variant detection and microsatellite instability status identification, there is little difference in using three different genomes. Using GGCPan as the reference had a significant advantage in structural variant identification. A total of 24 candidate gastric cancer driver genes were detected using three different reference genomes, of which eight were common and five were detected only based on pangenomes. Our results showed that disease-specific pangenome as a reference is promising and a whole set of tools are still to be developed or improved for disease genomics study in the pangenome era.

Genomics Research on the Road of Studying Biology and Virulence of Cereal Rust Fungi

Rust fungi are highly destructive pathogens that pose a significant threat to crop production worldwide, especially cereals. Obligate biotrophy and, in many cases, complex life cycles make rust fungi particularly challenging to study. However, recent rapid advances in sequencing technologies and genomic analysis tools have revolutionised rust fungal research. It is anticipated that the increasing availability and ongoing substantial improvements in genome assemblies will propel the field of rust biology into the post-genomic era, instigating a cascade of research endeavours encompassing multi-omics and gene discoveries. This is especially the case for many cereal rust pathogens, for which continental-scale studies of virulence have been conducted over many years and historical collections of viable isolates have been sequenced and assembled. Genomic analysis plays a crucial role in uncovering the underlying causes of the high variability of virulence and the complexity of population dynamics in rust fungi. Here, we provide an overview of progress in rust genomics, discuss the strategies employed in genomic analysis, and elucidate the strides that will drive cereal rust biology into the post-genomic era.

The Marchantia polymorpha pangenome reveals ancient mechanisms of plant adaptation to the environment

Plant adaptation to terrestrial life started 450 million years ago and has played a major role in the evolution of life on Earth. The genetic mechanisms allowing this adaptation to a diversity of terrestrial constraints have been mostly studied by focusing on flowering plants. Here, we gathered a collection of 133 accessions of the model bryophyte Marchantia polymorpha and studied its intraspecific diversity using selection signature analyses, a genome-environment association study and a pangenome. We identified adaptive features, such as peroxidases or nucleotide-binding and leucine-rich repeats (NLRs), also observed in flowering plants, likely inherited from the first land plants. The M. polymorpha pangenome also harbors lineage-specific accessory genes absent from seed plants. We conclude that different land plant lineages still share many elements from the genetic toolkit evolved by their most recent common ancestor to adapt to the terrestrial habitat, refined by lineage-specific polymorphisms and gene family evolution.

Meiotic crossovers revealed by differential visualization of homologous chromosomes using enhanced haplotype oligo-painting in cucumber

The interaction dynamics of homologous chromosomes during meiosis, such as recognition, pairing, synapsis, recombination, and segregation are vital for species fertility and genetic diversity within populations. Meiotic crossover (CO), a prominent feature of meiosis, ensures the faithful segregation of homologous chromosomes and enriches genetic diversity within a population. Nevertheless, visually distinguishing homologous chromosomes and COs remains an intractable challenge in cytological studies, particularly in non-model or plants with small genomes, limiting insights into meiotic dynamics. In the present study, we developed a robust and reliable enhanced haplotype oligo-painting (EHOP) technique to image small amounts of oligos, enabling visual discrimination of homologous chromosomes. Using EHOP developed based on sequence polymorphisms and reconstructed oligonucleotides, we visually distinguished parental and most recombinant chromosomes in cucumber F1 hybrids and F2 populations. Results from EHOP revealed that meiotic CO events preferentially occur in the 30-60% intervals of chromosome arms with lower sequence polymorphisms and significant recombination bias exists between cultivated and ancestral chromosomes. Due to the occupation of extensive heterochromatin occupancy, it is not yet possible to precisely identify the meiotic COs present in the central portion of chr2 and chr4. Notably, CO accessibility was universally detected in the cytological centromere region in F2 populations, a feature rarely observed in crops with large genomes. EHOP demonstrated exceptional performance in distinguishing homologous chromosomes and holds significant potential for broad application in studying homologous chromosome interactions.

Comparative analysis using a chromosome-scale genome assembly for Funaria hygrometrica suggests greater collinearity in mosses than in seed plants

Mosses, the largest lineage of seed-free plants, have smaller and less variable genome sizes than flowering plants. Nevertheless, whether this difference results from divergent genome dynamics is poorly known. Here, we use newly generated chromosome-scale genome assemblies for Funaria hygrometrica and comparative analysis with other moss and seed plant genomes to investigate moss genome dynamics. Although some aspects of moss genome dynamics are seed plant-like, such as the mechanism of genome size change and de novo gain/loss of genes, moss genomes retain higher synteny, and collinearity over evolutionary time than seed plant genomes. Furthermore, transposable elements and genes are more evenly distributed along chromosomes in mosses than in seed plants, a feature shared with other sequenced seed-free plant genomes. Overall, our findings support the hypothesis that large-scale genome structure and dynamics of mosses and seed plants differ. In particular, our data suggest a lower rate of gene order reshuffling along chromosomes in mosses compared to seed plants. We speculate that such lower rate of structural genomic variation and unique chromosome structure in mosses may contribute to their relatively smaller and less variable genome sizes.

Dispensable genome and segmental duplications drive the genome plasticity in Fusarium solani

Fusarium solani is a species complex encompassing a large phylogenetic clade with diverse members occupying varied habitats. We recently reported a unique opportunistic F. solani associated with unusual dark galls in sugarbeet. We assembled the chromosome-level genome of the F. solani sugarbeet isolate strain SB1 using Oxford Nanopore and Hi-C sequencing. The average size of F. solani genomes is 54 Mb, whereas SB1 has a larger genome of 59.38 Mb, organized into 15 chromosomes. The genome expansion of strain SB1 is due to the high repeats and segmental duplications within its three potentially accessory chromosomes. These chromosomes are absent in the closest reference genome with chromosome-level assembly, F. vanettenii 77-13-4. Segmental duplications were found in three chromosomes but are most extensive between two specific SB1 chromosomes, suggesting that this isolate may have doubled its accessory genes. Further comparison of the F. solani strain SB1 genome demonstrates inversions and syntenic regions to an accessory chromosome of F. vanettenii 77-13-4. The pan-genome of 12 publicly available F. solani isolates nearly reached gene saturation, with few new genes discovered after the addition of the last genome. Based on orthogroups and average nucleotide identity, F. solani is not grouped by lifestyle or origin. The pan-genome analysis further revealed the enrichment of several enzymes-coding genes within the dispensable (accessory + unique genes) genome, such as hydrolases, transferases, oxidoreductases, lyases, ligases, isomerase, and dehydrogenase. The evidence presented here suggests that genome plasticity, genetic diversity, and adaptive traits in Fusarium solani are driven by the dispensable genome with significant contributions from segmental duplications.

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.

Improvement of crop production in controlled environment agriculture through breeding

Controlled environment agriculture (CEA) represents one of the fastest-growing sectors of horticulture. Production in controlled environments ranges from highly controlled indoor environments with 100% artificial lighting (vertical farms or plant factories) to high-tech greenhouses with or without supplemental lighting, to simpler greenhouses and high tunnels. Although food production occurs in the soil inside high tunnels, most CEA operations use various hydroponic systems to meet crop irrigation and fertility needs. The expansion of CEA offers promise as a tool for increasing food production in and near urban systems as these systems do not rely on arable agricultural land. In addition, CEA offers resilience to climate instability by growing inside protective structures. Products harvested from CEA systems tend to be of high quality, both internal and external, and are sought after by consumers. Currently, CEA producers rely on cultivars bred for production in open-field agriculture. Because of high energy and other production costs in CEA, only a limited number of food crops have proven themselves to be profitable to produce. One factor contributing to this situation may be a lack of optimized cultivars. Indoor growing operations offer opportunities for breeding cultivars that are ideal for these systems. To facilitate breeding these specialized cultivars, a wide range of tools are available for plant breeders to help speed this process and increase its efficiency. This review aims to cover breeding opportunities and needs for a wide range of horticultural crops either already being produced in CEA systems or with potential for CEA production. It also reviews many of the tools available to breeders including genomics-informed breeding, marker-assisted selection, precision breeding, high-throughput phenotyping, and potential sources of germplasm suitable for CEA breeding. The availability of published genomes and trait-linked molecular markers should enable rapid progress in the breeding of CEA-specific food crops that will help drive the growth of this industry.

Pan-genome analyses of 11 Fraxinus species provide insights into salt adaptation in ash trees

Ash trees (Fraxinus) exhibit rich genetic diversity and wide adaptation to various ecological environments, and several species are highly salt tolerant. Dissecting the genomic basis of salt adaptation in Fraxinus is vital for its resistance breeding. Here, we present 11 high-quality chromosome-level genome assemblies for Fraxinus species, which reveal two unequal subgenome compositions and two recent whole-genome triplication events in their evolutionary history. A Fraxinus pan-genome was constructed on the basis of structural variations and revealed that presence-absence variations (PAVs) of transmembrane transport genes have likely contributed to salt adaptation in Fraxinus. Through whole-genome resequencing of an F1 population from an interspecies cross of F. velutina 'Lula 3' (salt tolerant) with F. pennsylvanica 'Lula 5' (salt sensitive), we mapped salt-tolerance PAV-based quantitative trait loci (QTLs) and pinpointed two PAV-QTLs and candidate genes associated with Fraxinus salt tolerance. Mechanistically, FvbHLH85 enhances salt tolerance by mediating reactive oxygen species and Na+/K+ homeostasis, whereas FvSWEET5 enhances salt tolerance by mediating osmotic homeostasis. Collectively, these findings provide valuable genomic resources for Fraxinus salt-resistance breeding and the research community.

Frequent genetic exchanges revealed by a pan-mitogenome graph of a fungal plant pathogen

Mitochondria are present in almost all eukaryotic lineages. The mitochondrial genomes (mitogenomes) evolve separately from nuclear genomes, and they can therefore provide relevant insights into the evolution of their host species. Fusarium oxysporum is a major fungal plant pathogen that is assumed to reproduce clonally. However, horizontal chromosome transfer between strains can occur through heterokaryon formation, and recently, signs of sexual recombination have been observed. Similarly, signs of recombination in F. oxysporum mitogenomes challenged the prevailing assumption of clonal reproduction in this species. Here, we construct, to our knowledge, the first fungal pan-mitogenome graph of nearly 500 F. oxysporum mitogenome assemblies to uncover the variation and evolution. In general, the gene order of fungal mitogenomes is not well conserved, yet the mitogenome of F. oxysporum and related species are highly colinear. We observed two strikingly contrasting regions in the F. oxysporum pan-mitogenome, comprising a highly conserved core mitogenome and a long variable region (6-16 kb in size), of which we identified three distinct types. The pan-mitogenome graph reveals that only five intron insertions occurred in the core mitogenome and that the long variable regions drive the difference between mitogenomes. Moreover, we observed that their evolution is neither concurrent with the core mitogenome nor with the nuclear genome. Our large-scale analysis of long variable regions uncovers frequent recombination between mitogenomes, even between strains that belong to different taxonomic clades. This challenges the common assumption of incompatibility between genetically diverse F. oxysporum strains and provides new insights into the evolution of this fungal species.IMPORTANCEInsights into plant pathogen evolution is essential for the understanding and management of disease. Fusarium oxysporum is a major fungal pathogen that can infect many economically important crops. Pathogenicity can be transferred between strains by the horizontal transfer of pathogenicity chromosomes. The fungus has been thought to evolve clonally, yet recent evidence suggests active sexual recombination between related isolates, which could at least partially explain the horizontal transfer of pathogenicity chromosomes. By constructing a pan-genome graph of nearly 500 mitochondrial genomes, we describe the genetic variation of mitochondria in unprecedented detail and demonstrate frequent mitochondrial recombination. Importantly, recombination can occur between genetically diverse isolates from distinct taxonomic clades and thus can shed light on genetic exchange between fungal strains.

A genome-wide association study of the racing performance traits in Yili horses based on Blink and FarmCPU models

Racing performance traits are the main indicators for evaluating the performance and value of sport horses. The aim of this study was to identify the key genes for racing performance traits in Yili horses by performing a genome-wide association study (GWAS). Breeding values for racing performance traits were calculated for Yili horses (n = 827) using an animal model. Genome-wide association analysis of racing performance traits in horses (n = 236) was carried out using the Blink, and FarmCPU models in GAPIT software, and genes within the significant regions were functionally annotated. The results of GWAS showed that a total of 24 significant SNP markers (P < 6.05 × 10- 9) and 22 suggestive SNP markers (P < 1.21 × 10- 7) were identified. Among them, the Blink associated 16 significant SNP loci and FarmCPU associated 12 significant SNP loci. A total of 127 candidate genes (50 significant) were annotated. Among these, CNTN6 (motor coordination), NIPA1 (neuronal development), and DCC (dopamine pathway maturation) may be the main candidate genes affecting speed traits. SHANK2 (neuronal synaptic regulation), ISCA1 (mitochondrial protein assembly), and KCNIP4 (neuronal excitability) may be the main candidate genes affecting ranking score traits. A common locus (ECA1: 22698579) was significantly associated with racing performance traits, and the function of the genes at this locus needs to be studied in depth. These findings will provide new insights into the detection and selection of genetic variants for racing performance and will help to accelerate the genetic improvement of Yili horses.

Pan-genome analysis of 13 Spinacia accessions reveals structural variations associated with sex chromosome evolution and domestication traits in spinach

Structural variations (SVs) are major genetic variants that can be involved in the origin, adaptation and domestication of species. However, the identification and characterization of SVs in Spinacia species are rare due to the lack of a pan-genome. Here, we report eight chromosome-scale assemblies of cultivated spinach and its two wild species. After integration with five existing assemblies, we constructed a comprehensive Spinacia pan-genome and identified 193 661 pan-SVs, which were genotyped in 452 Spinacia accessions. Our pan-SVs enabled genome-wide association study identified signals associated with sex and clarified the evolutionary direction of spinach. Most sex-linked SVs (86%) were biased to occur on the Y chromosome during the evolution of the sex-linked region, resulting in reduced Y-linked gene expression. The frequency of pan-SVs among Spinacia accessions further illustrated the contribution of these SVs to domestication, such as bolting time and seed dormancy. Furthermore, compared with SNPs, pan-SVs act as efficient variants in genomic selection (GS) because of their ability to capture missing heritability information and higher prediction accuracy. Overall, this study provides a valuable resource for spinach genomics and highlights the potential utility of pan-SV in crop improvement and breeding programmes.

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.

Analyzes of pan-genome and resequencing atlas unveil the genetic basis of jujube domestication

Jujube (Ziziphus jujuba Mill.), belonging to the Rhamnaceae family, is gaining increasing prominence as a perennial fruit crop with significant economic and medicinal values. Here, we conduct de novo assembly of four reference-grade genomes, encompassing one wild and three cultivated jujube accessions. We present insights into the population structure, genetic diversity, and genomic variations within a diverse collection of 1059 jujube accessions. Analyzes of the jujube pan-genome, based on our four assemblies and four previously released genomes, reveal extensive genomic variations within domestication-associated regions, potentially leading to the discovery of a candidate gene that regulates flowering and fruit ripening. By leveraging the pan-genome and a large-scale resequencing population, we identify two candidate genes involved in domestication traits, including the seed-setting rate, the bearing-shoot length and the leaf size in jujube. These genomic resources will accelerate evolutionary and functional genomics studies of jujube.

Low-dose 60Co-γ-ray irradiation promotes the growth of cucumber seedlings by inducing CsSAUR37 expression

Cucumber (Cucumis sativus L.) is a major vegetable crop grown globally, with a cultivation history of more than 3000 years. The limited genetic diversity, low rate of intraspecific variation, and extended periods of traditional breeding have resulted in slow progress in their genetic research and the development of new varieties. Gamma (γ)-ray irradiation potentially accelerates the breeding progress; however, the biological and molecular effects of γ-ray irradiation on cucumbers are unknown. Exposing cucumber seeds to 0, 50, 100, 150, 200, and 250 Gy doses of 60Co-γ-ray irradiation, this study aimed to investigate the resulting phenotype and physiological characteristics of seedling treatment to determine the optimal irradiation dose. The results showed that low irradiation doses (50-100 Gy) enhanced root growth, hypocotyl elongation, and lateral root numbers, promoting seedling growth. However, high irradiation doses (150-250 Gy) significantly inhibited seed germination and growth, decreasing the survival rate of seedlings. More than 100 Gy irradiation significantly decreased the total chlorophyll content while increasing the malondialdehyde (MDA) and H2O2 content in cucumber. Transcriptome sequencing analysis at 0, 50, 100, 150, 200, and 250 Gy doses showed that gene expression significantly differed between low and high irradiation doses. Gene Ontology enrichment and functional pathway enrichment analyses revealed that the auxin response pathway played a crucial role in seedling growth under low irradiation doses. Further, gene function analysis revealed that small auxin up-regulated gene CsSAUR37 was a key gene that was overexpressed in response to low irradiation doses, promoting primary root elongation and enhancing lateral root numbers by regulating the expression of protein phosphatase 2Cs (PP2Cs) and auxin synthesis genes.

Genetic architecture of kernel-related traits in sweet and waxy maize revealed by genome-wide association analysis

Introduction

Maize (Zea mays L.) is one of the most important crops worldwide, the kernel size-related traits are the major components of maize grain yield.

Methods

To dissect the genetic architecture of four kernel-related traits of 100-kernel weight, kernel length, kernel width, and kernel diameter, a genome-wide association study (GWAS) was conducted in the waxy and sweet maize panel comprising of 447 maize inbred lines re-sequenced at the 5× coverage depth. GWAS analysis was carried out with the mixed linear model using 1,684,029 high-quality SNP markers.

Results

In total, 49 SNPs significantly associated with the four kernel-related traits were identified, including 46 SNPs on chromosome 3, two SNPs on chromosome 4, and one SNP on chromosome 7. Haplotype regression analysis identified 338 haplotypes that significantly affected these four kernel-related traits. Genomic selection (GS) results revealed that a set of 10,000 SNPs and a training population size of 30% are sufficient for the application of GS in waxy and sweet maize breeding for kernel weight and kernel size. Forty candidate genes associated with the four kernel-related traits were identified, including both Zm00001d000707 and Zm00001d044139 expressed in the kernel development tissues and stages with unknown functions.

Discussion

These significant SNPs and important haplotypes provide valuable information for developing functional markers for the implementation of marker-assisted selection in breeding. The molecular mechanism of Zm00001d000707 and Zm00001d044139 regulating these kernel-related traits needs to be investigated further.

The CsDof1.8-CsLIPOXYGENASE09 module regulates C9 aroma production in cucumber

Nine-carbon aldehydes and their relative alcohols (C9 aromas) are the main aroma compounds of cucumber (Cucumis sativus L.) fruits and provide a unique cucumber-like note. However, the key regulators of C9 aroma accumulation in cucumber fruit are poorly characterized. Based on C9 aroma dynamic analysis and transcriptome analysis during fruit development of two different cucumber inbred lines, Q16 and Q24, Lipoxygenase09 (CsLOX09) was identified as a candidate gene for C9 aroma accumulation. Additionally, Q24 with higher CsLOX09 expression accumulated more C9 aromas than Q16. To verify the function of CsLOX09, Cslox09 homozygote knockout lines were created. C9 aroma content decreased by 80.79% to 99.16% in these mutants compared to the wild type. To further explore the reasons for the difference in CsLOX09 expression between Q16 and Q24 fruits, a co-expression network was constructed by integrating the C9 aroma-associated metabolism and transcriptomic data. Eighteen candidate transcription factors were highly correlated with the expression of CsLOX09. DNA binding with One Finger 1.8 (CsDof1.8) was confirmed to bind directly to the A/TAAAG motif of the CsLOX09 promoter through dual-luciferase, yeast one-hybrid, chromatin immunoprecipitation-qPCR and electrophoretic mobility shift assays. Furthermore, C9 aroma content and CsLOX09 expression were significantly increased in the CsDof1.8 overexpression lines. Overall, these data elucidate the metabolic regulation of C9 aromas in cucumber and provide a foundation for facilitating the regulation of flavor in cucumber breeding.

The era of panomics-driven gene discovery in plants

Panomics is an approach to integrate multiple 'omics' datasets, generated using different individuals or natural variations. Considering their diverse phenotypic spectrum, the phenome is inherently associated with panomics-based science, which is further combined with genomics, transcriptomics, metabolomics, and other omics techniques, either independently or collectively. Panomics has been accelerated through recent technological advancements in the field of genomics that enable the detection of population-wide structural variations (SVs) and hence offer unprecedented insights into the genetic variations contributing to important agronomic traits. The present review provides the recent trends of panomics-driven gene discovery toward various traits related to plant development, stress tolerance, accumulation of specialized metabolites, and domestication/dedomestication. In addition, the success stories are highlighted in the broader context of enhancing crop productivity.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

In-depth exploration of the genomic diversity in tea varieties based on a newly constructed pangenome of Camellia sinensis

Tea, one of the most widely consumed beverages globally, exhibits remarkable genomic diversity in its underlying flavour and health-related compounds. In this study, we present the construction and analysis of a tea pangenome comprising a total of 11 genomes, with a focus on three newly sequenced genomes comprising the purple-leaved assamica cultivar "Zijuan", the temperature-sensitive sinensis cultivar "Anjibaicha" and the wild accession "L618" whose assemblies exhibited excellent quality scores as they profited from latest sequencing technologies. Our analysis incorporates a detailed investigation of transposon complement across the tea pangenome, revealing shared patterns of transposon distribution among the studied genomes and improved transposon resolution with long read technologies, as shown by long terminal repeat (LTR) Assembly Index analysis. Furthermore, our study encompasses a gene-centric exploration of the pangenome, exploring the genomic landscape of the catechin pathway with our study, providing insights on copy number alterations and gene-centric variants, especially for Anthocyanidin synthases. We constructed a gene-centric pangenome by structurally and functionally annotating all available genomes using an identical pipeline, which both increased gene completeness and allowed for a high functional annotation rate. This improved and consistently annotated gene set will allow for a better comparison between tea genomes. We used this improved pangenome to capture the core and dispensable gene repertoire, elucidating the functional diversity present within the tea species. This pangenome resource might serve as a valuable resource for understanding the fundamental genetic basis of traits such as flavour, stress tolerance, and disease resistance, with implications for tea breeding programmes.

Machine learning-enabled high-throughput industry screening of edible oils

Long-term consumption of mixed fraudulent edible oils increases the risk of developing of chronic diseases which has been a threat to the public health globally. The complicated global supply-chain is making the industry malpractices had often gone undetected. In order to restore the confidence of consumers, traceability (and accountability) of every level in the supply chain is vital. In this work, we shown that machine learning (ML) assisted windowed spectroscopy (e.g., visible-band, infra-red band) produces high-throughput, non-destructive, and label-free authentication of edible oils (e.g., olive oils, sunflower oils), offers the feasibility for rapid analysis of large-scale industrial screening. We report achieving high-level of discriminant (AUC > 0.96) in the large-scale (n ≈ 11,500) of adulteration in olive oils. Notably, high clustering fidelity of 'spectral fingerprints' achieved created opportunity for (hypothesis-free) self-sustaining large database compilation which was never possible without machine learning. (137 words).

Genome Assembly and Structural Variation Analysis of Luffa acutangula Provide Insights on Flowering Time and Ridge Development

Luffa spp. is an important worldwide cultivated vegetable and medicinal plant from the Cucurbitaceae family. In this study, we report a high-quality chromosome-level genome of the high-generation inbred line SG261 of Luffa acutangula. The genomic sequence was determined by PacBio long reads, Hi-C sequencing reads, and 10× Genomics sequencing, with an assembly size of 739.82 Mb, contig N50 of 18.38 Mb, and scaffold N50 of 56.08 Mb. The genome of L. acutangula SG261 was predicted to contain 27,312 protein-coding genes and 72.56% repetitive sequences, of which long terminal repeats (LTRs) were an important form of repetitive sequences, accounting for 67.84% of the genome. Phylogenetic analysis reveals that L. acutangula evolved later than Luffa cylindrica, and Luffa is closely related to Momodica charantia. Comparing the genome of L. acutangula SG261 and L. cylindrica with PacBio data, 67,128 high-quality structural variations (SVs) and 55,978 presence-absence variations (PAVs) were identified in SG261, resulting in 2424 and 1094 genes with variation in the CDS region, respectively, and there are 287 identical genes affected by two different structural variation analyses. In addition, we found that the transcription factor FY (FLOWERING LOCUS Y) families had a large expansion in L. acutangula SG261 (flowering in the morning) compared to L. cylindrica (flowering in the afternoon), which may result in the early flowering time in L. acutangula SG261. This study provides valuable reference for the breeding of and pan-genome research into Luffa species.

Genetic characterization of cucumber genetic resources in the NARO Genebank indicates their multiple dispersal trajectories to the East

KEY MESSAGE

Genotyping-by-sequencing of 723 worldwide cucumber genetic resources revealed that cucumbers were dispersed eastward via at least three distinct routes, one to Southeast Asia and two from different directions to East Asia. The cucumber (Cucumis sativus) is an economically important vegetable crop cultivated and consumed worldwide. Despite its popularity, the manner in which cucumbers were dispersed from their origin in South Asia to the rest of the world, particularly to the east, remains a mystery due to the lack of written records. In this study, we performed genotyping-by-sequencing (GBS) on 723 worldwide cucumber accessions, mainly deposited in the Japanese National Agriculture and Food Research Organization (NARO) Genebank, to characterize their genetic diversity, relationships, and population structure. Analyses based on over 60,000 genome-wide single-nucleotide polymorphisms identified by GBS revealed clear genetic differentiation between Southeast and East Asian populations, suggesting that they reached their respective region independently, not progressively. A deeper investigation of the East Asian population identified two subpopulations with different fruit characteristics, supporting the traditional classification of East Asian cucumbers into two types thought to have been introduced by independent routes. Finally, we developed a core collection of 100 accessions representing at least 93.2% of the genetic diversity present in the entire collection. The genetic relationships and population structure, their associations with geographic distribution and phenotypic traits, and the core collection presented in this study are valuable resources for elucidating the dispersal history and promoting the efficient use and management of genetic resources for research and breeding in cucumber.

Graph-Based Pangenome of Actinidia chinensis Reveals Structural Variations Mediating Fruit Degreening

Fruit ripening is associated with the degreening process (loss of chlorophyll) that occurs in most fruit species. Kiwifruit is one of the special species whose fruits may maintain green flesh by accumulating a large amount of chlorophyll even after ripening. However, little is known about the genetic variations related to the fruit degreening process. Here, a graph-based kiwifruit pangenome by analyzing 14 chromosome-scale haplotype-resolved genome assemblies from seven representative cultivars or lines in Actinidia chinensis is built. A total of 49,770 non-redundant gene families are identified, with core genes constituting 46.6%, and dispensable genes constituting 53.4%. A total of 84,591 non-redundant structural variations (SVs) are identified. The pangenome graph integrating both reference genome sequences and variant information facilitates the identification of SVs related to fruit color. The SV in the promoter of the AcBCM gene determines its high expression in the late developmental stage of fruits, which causes chlorophyll accumulation in the green-flesh fruits by post-translationally regulating AcSGR2, a key enzyme of chlorophyll catabolism. Taken together, a high-quality pangenome is constructed, unraveled numerous genetic variations, and identified a novel SV mediating fruit coloration and fruit quality, providing valuable information for further investigating genome evolution and domestication, QTL genes function, and genomics-assisted breeding.

An alignment- and reference-free strategy using k-mer present pattern for population genomic analyses

Pangenomes are replacing single reference genomes to capture all variants within a species or clade, but their analysis predominantly leverages graph-based methods that require multiple high-quality genomes and computationally intensive multiple-genome alignments. K-mer decomposition is an alternative to graph-based pangenomes. However, how to directly use k-mers for the population genetic analyses is unknown. Here, we developed a novel strategy that uses the variants of k-mer count in the genome for population analyses. To test the effectivity of this method, we compared it directly to the SNP-based method on the analysis of population structure and genetic diversity of 267 Saccharomyces cerevisiae strains within two simulated datasets and a real sequence dataset. The population structure identified with k-mers recapitulates that obtained using SNPs, indicating the effectiveness of k-mer-based approach, and higher genetic diversity within real dataset supported k-mers contained more genetic variants. Based on k-mer frequency, we found not only SNP but also some insertion/deletion and horizontal gene transfer (HGT) fragments related to the adaptive evolution of S. cerevisiae. Our study creates a framework for the alignment- and reference-free (ARF) method in population genetic analyses, which will be more pronounced in the species with no complete genome or highly diverged species.

A nanopore-based cucumber genome assembly reveals structural variations at two QTLs controlling hypocotyl elongation

The Xishuangbanna (XIS) cucumber (Cucumis sativus var. xishuangbannanesis) is a semiwild variety that has many distinct agronomic traits. Here, long reads generated by Nanopore sequencing technology helped assembling a high-quality genome (contig N50 = 8.7 Mb) of landrace XIS49. A total of 10,036 structural/sequence variations (SVs) were identified when comparing with Chinese Long (CL), and known SVs controlling spines, tubercles, and carpel number were confirmed in XIS49 genome. Two QTLs of hypocotyl elongation under low light, SH3.1 and SH6.1, were fine-mapped using introgression lines (donor parent, XIS49; recurrent parent, CL). SH3.1 encodes a red-light receptor Phytochrome B (PhyB, CsaV3_3G015190). A ∼4 kb region with large deletion and highly divergent regions (HDRs) were identified in the promoter of the PhyB gene in XIS49. Loss of function of this PhyB caused a super-long hypocotyl phenotype. SH6.1 encodes a CCCH-type zinc finger protein FRIGIDA-ESSENTIAL LIKE (FEL, CsaV3_6G050300). FEL negatively regulated hypocotyl elongation but it was transcriptionally suppressed by long terminal repeats retrotransposon insertion in CL cucumber. Mechanistically, FEL physically binds to the promoter of CONSTITUTIVE PHOTOMORPHOGENIC 1a (COP1a), regulating the expression of COP1a and the downstream hypocotyl elongation. These above results demonstrate the genetic mechanism of cucumber hypocotyl elongation under low light.

A wild melon reference genome provides novel insights into the domestication of a key gene responsible for melon fruit acidity

KEY MESSAGE

A wild melon reference genome elucidates the genomic basis of fruit acidity domestication. Structural variants (SVs) have been reported to impose major effects on agronomic traits, representing a significant contributor to crop domestication. However, the landscape of SVs between wild and cultivated melons is elusive and how SVs have contributed to melon domestication remains largely unexplored. Here, we report a 379-Mb chromosome-scale genome of a wild progenitor melon accession "P84", with a contig N50 of 14.9 Mb. Genome comparison identifies 10,589 SVs between P84 and four cultivated melons with 6937 not characterized in previously analysis of 25 melon genome sequences. Furthermore, the population-scale genotyping of these SVs was determined in 1175 accessions, and 18 GWAS signals including fruit acidity, fruit length, fruit weight, fruit color and sex determination were detected. Based on these genotyped SVs, we identified 3317 highly diverged SVs between wild and cultivated melons, which could be the potential SVs associated with domestication-related traits. Furthermore, we identify novel SVs affecting fruit acidity and proposed the diverged evolutionary trajectories of CmPH, a key regulator of melon fruit acidity, during domestication and selection of different populations. These results will offer valuable resources for genomic studies and genetic improvement in melon.

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

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.

Multi-Omics Analysis Reveals the Distinct Features of Metabolism Pathways Supporting the Fruit Size and Color Variation of Giant Pumpkin

Pumpkin (Cucurbita maxima) is an important vegetable crop of the Cucurbitaceae plant family. The fruits of pumpkin are often used as directly edible food or raw material for a number of processed foods. In nature, mature pumpkin fruits differ in size, shape, and color. The Atlantic Giant (AG) cultivar has the world's largest fruits and is described as the giant pumpkin. AG is well-known for its large and bright-colored fruits with high ornamental and economic value. At present, there are insufficient studies that have focused on the formation factors of the AG cultivar. To address these knowledge gaps, we performed comparative transcriptome, proteome, and metabolome analysis of fruits from the AG cultivar and a pumpkin with relatively small fruit (Hubbard). The results indicate that up-regulation of gene-encoded expansins contributed to fruit cell expansion, and the increased presence of photoassimilates (stachyose and D-glucose) and jasmonic acid (JA) accumulation worked together in terms of the formation of large fruit in the AG cultivar. Notably, perhaps due to the rapid transport of photoassimilates, abundant stachyose that was not converted into glucose in time was detected in giant pumpkin fruits, implying that a unique mode of assimilate unloading is in existence in the AG cultivar. The potential molecular regulatory network of photoassimilate metabolism closely related to pumpkin fruit expansion was also investigated, finding that three MYB transcription factors, namely CmaCh02G015900, CmaCh01G018100, and CmaCh06G011110, may be involved in metabolic regulation. In addition, neoxanthin (a type of carotenoid) exhibited decreased accumulation that was attributed to the down-regulation of carotenoid biosynthesis genes in AG fruits, which may lead to pigmentation differences between the two pumpkin cultivars. Our current work will provide new insights into the potential formation factors of giant pumpkins for further systematic elucidation.

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

Pan-Genome Analysis of TIFY Gene Family and Functional Analysis of CsTIFY Genes in Cucumber

Cucumbers are frequently affected by gray mold pathogen Botrytis cinerea, a pathogen that causes inhibited growth and reduced yield. Jasmonic acid (JA) plays a primary role in plant responses to biotic stresses, and the jasmonate-ZIM-Domain (JAZ) proteins are key regulators of the JA signaling pathway. In this study, we used the pan-genome of twelve cucumber varieties to identify cucumber TIFY genes. Our findings revealed that two CsTIFY genes were present in all twelve cucumber varieties and showed no differences in protein sequence, gene structure, and motif composition. This suggests their evolutionary conservation across different cucumber varieties and implies that they may play a crucial role in cucumber growth. On the other hand, the other fourteen CsTIFY genes exhibited variations in protein sequence and gene structure or conserved motifs, which could be the result of divergent evolution, as these genes adapt to different cultivation and environmental conditions. Analysis of the expression profiles of the CsTIFY genes showed differential regulation by B. cinerea. Transient transfection plants overexpressing CsJAZ2, CsJAZ6, or CsZML2 were found to be more susceptible to B. cinerea infection compared to control plants. Furthermore, these plants infected by the pathogen showed lower levels of the enzymatic activities of POD, SOD and CAT. Importantly, after B. cinerea infection, the content of JA was upregulated in the plants, and cucumber cotyledons pretreated with exogenous MeJA displayed increased resistance to B. cinerea infection compared to those pretreated with water. Therefore, this study explored key TIFY genes in the regulation of cucumber growth and adaptability to different cultivation environments based on bioinformatics analysis and demonstrated that CsJAZs negatively regulate cucumber disease resistance to gray mold via multiple signaling pathways.

Pan-genome Analysis of WOX Gene Family and Function Exploration of CsWOX9 in Cucumber

Cucumber is an economically important vegetable crop, and the warts (composed of spines and Tubercules) of cucumber fruit are an important quality trait that influences its commercial value. WOX transcription factors are known to have pivotal roles in regulating various aspects of plant growth and development, but their studies in cucumber are limited. Here, genome-wide identification of cucumber WOX genes was performed using the pan-genome analysis of 12 cucumber varieties. Our findings revealed diverse CsWOX genes in different cucumber varieties, with variations observed in protein sequences and lengths, gene structure, and conserved protein domains, possibly resulting from the divergent evolution of CsWOX genes as they adapt to diverse cultivation and environmental conditions. Expression profiles of the CsWOX genes demonstrated that CsWOX9 was significantly expressed in unexpanded ovaries, especially in the epidermis. Additionally, analysis of the CsWOX9 promoter revealed two binding sites for the C2H2 zinc finger protein. We successfully executed a yeast one-hybrid assay (Y1H) and a dual-luciferase (LUC) transaction assay to demonstrate that CsWOX9 can be transcriptionally activated by the C2H2 zinc finger protein Tu, which is crucial for fruit Tubercule formation in cucumber. Overall, our results indicated that CsWOX9 is a key component of the molecular network that regulates wart formation in cucumber fruits, and provide further insight into the function of CsWOX genes in cucumber.

Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis

Tea is one of the world's oldest crops and is cultivated to produce beverages with various flavours. Despite advances in sequencing technologies, the genetic mechanisms underlying key agronomic traits of tea remain unclear. In this study, we present a high-quality pangenome of 22 elite cultivars, representing broad genetic diversity in the species. Our analysis reveals that a recent long terminal repeat burst contributed nearly 20% of gene copies, introducing functional genetic variants that affect phenotypes such as leaf colour. Our graphical pangenome improves the efficiency of genome-wide association studies and allows the identification of key genes controlling bud flush timing. We also identified strong correlations between allelic variants and flavour-related chemistries. These findings deepen our understanding of the genetic basis of tea quality and provide valuable genomic resources to facilitate its genomics-assisted breeding.

PanPA: generation and alignment of panproteome graphs

Motivation

Compared to eukaryotes, prokaryote genomes are more diverse through different mechanisms, including a higher mutation rate and horizontal gene transfer. Therefore, using a linear representative reference can cause a reference bias. Graph-based pangenome methods have been developed to tackle this problem. However, comparisons in DNA space are still challenging due to this high diversity. In contrast, amino acid sequences have higher similarity due to evolutionary constraints, whereby a single amino acid may be encoded by several synonymous codons. Coding regions cover the majority of the genome in prokaryotes. Thus, panproteomes present an attractive alternative leveraging the higher sequence similarity while not losing much of the genome in non-coding regions.

Results

We present PanPA, a method that takes a set of multiple sequence alignments of protein sequences, indexes them, and builds a graph for each multiple sequence alignment. In the querying step, it can align DNA or amino acid sequences back to these graphs. We first showcase that PanPA generates correct alignments on a panproteome from 1350 Escherichia coli. To demonstrate that panproteomes allow comparisons at longer phylogenetic distances, we compare DNA and protein alignments from 1073 Salmonella enterica assemblies against E.coli reference genome, pangenome, and panproteome using BWA, GraphAligner, and PanPA, respectively; with PanPA aligning around 22% more sequences. We also aligned a DNA short-reads whole genome sequencing (WGS) sample from S.enterica against the E.coli reference with BWA and the panproteome with PanPA, where PanPA was able to find alignment for 68% of the reads compared to 5% with BWA.

Availalability and implementation

PanPA is available at https://github.com/fawaz-dabbaghieh/PanPA.

Comprehensive dissection of meiotic DNA double-strand breaks and crossovers in cucumber

Meiotic recombination drives genetic diversity and crop genome optimization. In plant breeding, parents with favorable traits are crossed to create elite varieties. Different hybridizations produce diverse types of segment reshuffling between homologous chromosomes. However, little is known about the factors that cause hybrid-specific changes in crossovers (COs). Here, we constructed 2 F2 populations from crosses between a semiwild and 2 domesticated cucumber (Cucumis sativus) accessions and examined CO events. COs mainly occurred around genes and differed unevenly along chromosomes between the 2 hybrids. Fine-scale CO distributions were suppressed in regions of heterozygous structural variations (SVs) and were accelerated by high sequence polymorphism. C. sativus RADiation sensitive 51A (CsRAD51A) binding, histone H3 lysine 4 trimethylation (H3K4me3) modification, chromatin accessibility, and hypomethylation were positively associated with global CO landscapes and in local DNA double-strand break (DSB) hotspots and genes. The frequency and suppression of COs could be roughly predicted based on multiomic information. Differences in CO events between hybrids could be partially traced to distinct genetic and epigenetic features and were significantly associated with specific DSB hotspots and heterozygous SVs. Our findings identify the genomic and epigenetic features that contribute to CO formation and hybrid-specific divergence in cucumber and provide theoretical support for selecting parental combinations and manipulating recombination events at target genomic regions during plant breeding.

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/.

Pan-genome analysis sheds light on structural variation-based dissection of agronomic traits in melon crops

Sweetness and appearance of fresh fruits are key palatable and preference attributes for consumers and are often controlled by multiple genes. However, fine-mapping the key loci or genes of interest by single genome-based genetic analysis is challenging. Herein, we present the chromosome-level genome assembly of 1 landrace melon accession (Cucumis melo ssp. agrestis) with wild morphologic features and thus construct a melon pan-genome atlas via integrating sequenced melon genome datasets. Our comparative genomic analysis reveals a total of 3.4 million genetic variations, of which the presence/absence variations (PAVs) are mainly involved in regulating the function of genes for sucrose metabolism during melon domestication and improvement. We further resolved several loci that are accountable for sucrose contents, flesh color, rind stripe, and suture using a structural variation (SV)-based genome-wide association study. Furthermore, via bulked segregation analysis (BSA)-seq and map-based cloning, we uncovered that a single gene, (CmPIRL6), determines the edible or inedible characteristics of melon fruit exocarp. These findings provide important melon pan-genome information and provide a powerful toolkit for future pan-genome-informed cultivar breeding of melon.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.