A pangenome analysis pipeline provides insights into functional gene identification in rice

The developments and prospects of plant super-pangenomes: Demands, approaches, and applications

By integrating genomes from different accessions, pangenomes provide a more comprehensive and reference-bias-free representation of genetic information within a population compared to a single reference genome. With the rapid accumulation of genomic sequencing data and the expanding scope of plant research, plant pangenomics has gradually evolved from single-species to multi-species studies. This shift has given rise to the concept of a super-pangenome that covers all genomic sequences within a genus-level taxonomic group. By incorporating both cultivated and wild species, the super-pangenome has greatly enhanced the resolution of research in various areas such as plant genetic diversity, evolution, domestication, and molecular breeding. In this review, we present a comprehensive overview of the plant super-pangenome, emphasizing its development requirements, construction strategies, potential applications, and notable achievements. We also highlight the distinctive advantages and promising prospects of super-pangenomes while addressing current challenges and future directions.

A barley pan-transcriptome reveals layers of genotype-dependent transcriptional complexity

A pan-transcriptome describes the transcriptional and post-transcriptional consequences of genome diversity from multiple individuals within a species. We developed a barley pan-transcriptome using 20 inbred genotypes representing domesticated barley diversity by generating and analyzing short- and long-read RNA-sequencing datasets from multiple tissues. To overcome single reference bias in transcript quantification, we constructed genotype-specific reference transcript datasets (RTDs) and integrated these into a linear pan-genome framework to create a pan-RTD, allowing transcript categorization as core, shell or cloud. Focusing on the core (expressed in all genotypes), we observed significant transcript abundance variation among tissues and between genotypes driven partly by RNA processing, gene copy number, structural rearrangements and conservation of promotor motifs. Network analyses revealed conserved co-expression module::tissue correlations and frequent functional diversification. To complement the pan-transcriptome, we constructed a comprehensive cultivar (cv.) Morex gene-expression atlas and illustrate how these combined datasets can be used to guide biological inquiry.

Natural variation in CTF1 conferring cold tolerance at the flowering stage in rice

Improving cold tolerance at the flowering stage (CTF) in rice is crucial for minimising yield loss, making the identification and application of cold-tolerant genes and QTLs imperative for effective molecular breeding. The long lead time, dependence on cold treatment conditions, and the inherent complexity of the trait make studying the genetic basis of CTF in rice challenging. To date, the fine-mapping or cloning of QTLs specific to CTF has not yet been achieved. In this study, single segment substitution lines (SSSLs) were constructed using HJX74 as the recipient and IR58025B, known for good CTF, as the donor. This approach led to the identification of two cold tolerance QTLs, qCTF3 and qCTF6, in rice. qCTF6 has promising breeding potential. Further, we identified the causal gene CTF1 underlying qCTF6 through map-based cloning. CTF1 which encodes a conserved putative protein, has two SNPs within its coding sequence that influence CTF in rice. Additionally, genetic variations in the promoter of CTF1 also contributes to CTF. Thirteen variant sites of CTF1 in the four cold tolerance SSSLs are consistent with the IR58025B. Moreover, we analysed 307 accessions to characterise haplotypes based on the 13 variation sites, identifying five distinct haplotypes. The selection and evolutionary analysis indicate that the cold-tolerant haplotype of CTF1 is a newly generated mutation that has undergone selection in japonica during domestication. This study not only provides a novel favourable gene for molecular breeding of CTF but also highlights the potential of CTF1 in advancing rice breeding.

The role of pangenomics in orphan crop improvement

Global food security depends heavily on a few staple crops, while orphan crops, despite being less studied, offer the potential benefits of environmental adaptation and enhanced nutritional traits, especially in a changing climate. Major crops have benefited from genomics-based breeding, initially using single genomes and later pangenomes. Recent advances in DNA sequencing have enabled pangenome construction for several orphan crops, offering a more comprehensive understanding of genetic diversity. Orphan crop research has now entered the pangenomics era and applying these pangenomes with advanced selection methods and genome editing technologies can transform these neglected species into crops of broader agricultural significance.

Genomic and cell-specific regulation of benzylisoquinoline alkaloid biosynthesis in opium poppy

Opium poppy is a crop of great commercial value as a source of several opium alkaloids for the pharmaceutical industries including morphine, codeine, thebaine, noscapine, and papaverine. Most enzymes involved in benzylisoquinoline alkaloid (BIA) biosynthesis in opium poppy have been functionally characterized, and opium poppy currently serves as a model system to study BIA metabolism in plants. BIA biosynthesis in opium poppy involves two biosynthetic gene clusters associated respectively with the morphine and noscapine branches. Recent reports have shown that genes in the same cluster are co-expressed, suggesting they might also be co-regulated. However, the transcriptional regulation of opium poppy BIA biosynthesis is not well studied. Opium poppy BIA biosynthesis involves three cell types associated with the phloem system: companion cells, sieve elements, and laticifers. The transcripts and enzymes associated with BIA biosynthesis are distributed across cell types, requiring the translocation of key enzymes and pathway intermediates between cell types. Together, these suggest that the regulation of BIA biosynthesis in opium poppy is multilayered and complex, involving biochemical, genomic, and physiological mechanisms. In this review, we highlight recent advances in genome sequencing and single cell and spatial transcriptomics with a focus on how these efforts can improve our understanding of the genomic and cell-specific regulation of BIA biosynthesis. Such knowledge is vital for opium poppy genetic improvement and metabolic engineering efforts targeting the modulation of alkaloid yield and composition.

Centromeres are hotspots for chromosomal inversions and breeding traits in mango

Chromosomal inversions can preserve combinations of favorable alleles by suppressing recombination. Simultaneously, they reduce the effectiveness of purifying selection enabling deleterious alleles to accumulate. This study explores how areas of low recombination, including centromeric regions and chromosomal inversions, contribute to the accumulation of deleterious and favorable loci in 225 Mangifera indica genomes from the Australian Mango Breeding Program. Here, we identify 17 chromosomal inversions that cover 7.7% (29.7 Mb) of the M. indica genome: eight pericentric (inversion includes the centromere) and nine paracentric (inversion is on one arm of the chromosome). Our results show that these large pericentric inversions are accumulating deleterious loci, while the paracentric inversions show deleterious levels above and below the genome wide average. We find that despite their deleterious load, chromosomal inversions contain small effect loci linked to variation in crucial breeding traits. These results indicate that chromosomal inversions have likely facilitated the evolution of key mango breeding traits. Our study has important implications for selective breeding of favorable combinations of alleles in regions of low recombination.

Regulatory networks of coresident subgenomes during rapid fiber cell elongation in upland cotton

Cotton, an intriguing plant species shaped by polyploidization, evolution, and domestication, holds particular interest due to the complex mechanisms governing fiber traits across its two subgenomes. However, the regulatory elements or transcriptional networks between subgenomes during fiber elongation remain to be fully clarified. Here, we analyzed 1462 cotton fiber samples to reconstruct the gene-expression regulatory networks that influence fiber cell elongation. Inter-subgenome expression quantitative trait loci (eQTLs) largely dictate gene transcription, with a notable tendency for the D subgenome to regulate A-subgenome eGenes. This regulation reveals synchronized homoeologous gene expression driven by co-localized eQTLs and divergent patterns that diminish genetic correlations, thus leading to preferential expression in the A and D subgenomes. Hotspot456 emerged as a key regulator of fiber initiation and elongation, and artificial selection of trans-eQTLs in hotspot456 that positively regulate KCS1 has facilitated cell elongation. Experiments designed to clarify the roles of trans-eQTLs in improved fiber breeding confirmed the inhibition of GhTOL9 by a specific trans-eQTL via GhWRKY28, which negatively affects fiber elongation. We propose a model in which the GhWRKY28-GhTOL9 module regulates this process through the ESCRT (endosomal sorting complex required for transport) pathway. This research significantly advances our understanding of cotton's evolutionary and domestication processes and the intricate regulatory mechanisms that underlie significant plant traits.

Geography, Antimicrobial Resistance, and Genomics of Salmonella enterica (Serotypes Newport and Anatum) from Meat in Mexico (2021-2023)

Salmonella enterica non-typhoidal is a major contributor to diarrheal diseases, with over 2600 serovars identified across diverse environments. In Mexico, serovars Newport and Anatum have shown a marked increase, especially in foodborne disease, posing a public health problem. We conducted a cross-sectional study from 2021 to 2023 using active epidemiological surveillance to assess contamination in ground beef and pork at butcher shops nationwide. It involved isolation, phenotypic antimicrobial resistance, comparative genomics, spatial distribution, antimicrobial-resistance genes, and pangenome analysis. A total of 402 non-typhoidal S. enterica strains were isolated, including 59 Newport and 50 Anatum. After curating for redundancy, 45 Newport and 32 Anatum strains remained. We found that 75% of Newport strains exhibited multidrug resistance (MDR), compared to 25% of Anatum strains. Salmonella Newport also showed a broader distribution and stronger antibiotic-resistance capacity, particularly due to genes such as mphA and ramA. Our pangenome analysis showed a predominance of cell maintenance and survival-process genes in the accessory genome of both serotypes. Considering unique genes, Salmonella Anatum and Newport showed a notorious abundance of genes with functions related to replication, recombination, and repair. The substantial rise of Anatum and Newport strains in meat samples for human consumption presents an epidemiological alert, highlighting the critical need for stringent surveillance programs to mitigate human and ecosystem health risks.

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.

Unveiling the whole genomic features and potential probiotic characteristics of novel Lactiplantibacillus plantarum HMX2

This study investigates the genomic features and probiotic potential of Lactiplantibacillus plantarum HMX2, isolated from Chinese Sauerkraut, using whole-genome sequencing (WGS) and bioinformatics for the first time. This study also aims to find genetic diversity, antibiotic resistance genes, and functional capabilities to help us better understand its food safety applications and potential as a probiotic. L. plantarum HMX2 was cultured, and DNA was extracted for WGS. Genomic analysis comprised average nucleotide identity (ANI) prediction, genome annotation, pangenome, and synteny analysis. Bioinformatics techniques were used to identify CoDing Sequences (CDSs), transfer RNA (tRNA) and ribosomal RNA (rRNA) genes, and antibiotic resistance genes, as well as to conduct phylogenetic analysis to establish genetic diversity and evolution. The study found a significant genetic similarity (99.17% ANI) between L. plantarum HMX2 and the reference strain. Genome annotation revealed 3,242 coding sequences, 65 tRNA genes, and 16 rRNA genes. Significant genetic variety was found, including 25 antibiotic resistance genes. A phylogenetic study placed L. plantarum HMX2 among closely related bacteria, emphasizing its potential for probiotic and food safety applications. The genomic investigation of L. plantarum showed essential genes, including plnJK and plnEF, which contribute to antibacterial action against foodborne pathogens. Furthermore, genes such as MurA, Alr, and MprF improve food safety and probiotic potential by promoting bacterial survival under stress conditions in food and the gastrointestinal tract. This study introduces the new genomic features of L. plantarum HMX2 about specific genetics and its possibility of relevant uses in food security and technologies. These findings of specific genes involved in antimicrobial activity provide fresh possibilities for exploiting this strain in forming probiotic preparations and food preservation methods. The future research should focus on the experimental validation of antibiotic resistance genes, comparative genomics to investigate functional diversity, and the development of novel antimicrobial therapies that take advantage of L. plantarum's capabilities.

Historic breeding practices contribute to germplasm divergence in leaf specialized metabolism and ecophysiology in cultivated sunflower (Helianthus annuus)

PREMISE

The use of hybrid breeding systems to increase crop yields has been the cornerstone of modern agriculture and is exemplified in the breeding and improvement of cultivated sunflower (Helianthus annuus). However, it is poorly understood what effect supporting separate breeding pools in such systems, combined with continued selection for yield, may have on leaf ecophysiology and specialized metabolite variation.

METHODS

We analyzed 288 lines of cultivated H. annuus to examine the genomic basis of several specialized metabolites and agronomically important traits across major heterotic groups.

RESULTS

Heterotic group identity supports phenotypic divergences between fertility restoring and cytoplasmic male-sterility maintainer lines in leaf ecophysiology and specialized metabolism. However, the divergence is not associated with physical linkage to nuclear genes that support current hybrid breeding practices in cultivated H. annuus. Additionally, we identified four genomic regions associated with leaf ecophysiology and specialized metabolism that colocalize with previously identified QTLs for quantitative self-compatibility traits and with S-protein homolog (SPH) proteins, a recently discovered family of proteins associated with self-incompatibility and self/nonself recognition in Papaver rhoeas (common poppy) with suggested conserved downstream mechanisms among eudicots.

CONCLUSIONS

Further work is necessary to confirm the self-incompatibility mechanisms in cultivated H. annuus and their relationship to the integrative and polygenic architecture of leaf ecophysiology and specialized metabolism in cultivated sunflower. However, because self-compatibility is a derived quantitative trait in cultivated H. annuus, trait linkage to divergent phenotypic traits may have partially arisen as a potential unintended consequence of historical breeding practices and selection for yield.

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.

Previously unmeasured genetic diversity explains part of Lewontin's paradox in a k -mer-based meta-analysis of 112 plant species

At the molecular level, most evolution is expected to be neutral. A key prediction of this expectation is that the level of genetic diversity in a population should scale with population size. However, as was noted by Richard Lewontin in 1974 and reaffirmed by later studies, the slope of the population size-diversity relationship in nature is much weaker than expected under neutral theory. We hypothesize that one contributor to this paradox is that current methods relying on single nucleotide polymorphisms (SNPs) called from aligning short reads to a reference genome underestimate levels of genetic diversity in many species. To test this idea, we calculated nucleotide diversity ( π ) and k -mer-based metrics of genetic diversity across 112 plant species, amounting to over 205 terabases of DNA sequencing data from 27,488 individual plants. We then compared how these different metrics correlated with proxies of population size that account for both range size and population density variation across species. We found that our population size proxies scaled anywhere from about 3 to over 20 times faster with k -mer diversity than nucleotide diversity after adjusting for evolutionary history, mating system, life cycle habit, cultivation status, and invasiveness. The relationship between k -mer diversity and population size proxies also remains significant after correcting for genome size, whereas the analogous relationship for nucleotide diversity does not. These results suggest that variation not captured by common SNP-based analyses explains part of Lewontin's paradox in plants.

Endoplasmic reticulum stress response modulator OsbZIP39 regulates cadmium accumulation via activating the expression of defensin-like gene OsCAL2 in rice

Accumulation of cadmium (Cd) in rice is not only harmful to the growth of plants but also poses a threat to human health. Exposure to Cd triggers unfolded protein response (UPR) within cells, a process that is still not completely understood. The study demonstrated that the lack of OsbZIP39, an essential endoplasmic reticulum (ER)-resident regulator of the UPR, resulted in decreased Cd intake and reduced Cd levels in the roots, stems, and grains of rice. Upon exposure to Cd stress, GFP-OsbZIP39 translocated from ER to nucleus, initiating UPR. Further investigation revealed that Cd treatment caused changes in sphingolipid levels in the membrane, influencing the localization and activation of OsbZIP39. Yeast one-hybrid and dual-LUC assays were conducted to validate the interaction between activated OsbZIP39 and the promoter of the defensin-like gene OsCAL2, resulting in an increase in its expression. Different variations were identified in the coding region of OsbZIP39, which may explain the varying levels of Cd accumulation observed in the indica and japonica subspecies. Under Cd treatment, OsbZIP39ind exhibited a more significant enhancement in the transcription of OsCAL2 compared to OsbZIP39jap. Our data suggest that OsbZIP39 positively regulates Cd uptake in rice, offering an encouraging objective for the cultivation of low-Cd rice.

Triumphs of genomic-assisted breeding in crop improvement

Conventional breeding approaches have played a significant role in meeting the food demand remarkably well until now. However, the increasing population, yield plateaus in certain crops, and limited recombination necessitate using genomic resources for genomics-assisted crop improvement programs. As a result of advancements in the next-generation sequence technology, GABs have developed dramatically to characterize allelic variants and facilitate their rapid and efficient incorporation in crop improvement programs. Genomics-assisted breeding (GAB) has played an important role in harnessing the potential of modern genomic tools, exploiting allelic variation from genetic resources and developing cultivars over the past decade. The availability of pangenomes for major crops has been a significant development, albeit with varying degrees of completeness. Even though adopting these technologies is essentially determined on economic grounds and cost-effective assays, which create a wealth of information that can be successfully used to exploit the latent potential of crops. GAB has been instrumental in harnessing the potential of modern genomic resources and exploiting allelic variation for genetic enhancement and cultivar development. GAB strategies will be indispensable for designing future crops and are expected to play a crucial role in breeding climate-smart crop cultivars with higher nutritional value.

A Novel Function of GW5 on Controlling the Early Growth Vigor and its Haplotype Effect on Shoot Dry Weight and Grain Size in Rice (Oryza sativa L.)

Strong early growth vigor is an essential target in both direct seeded rice breeding and high-yielding rice breeding for rice varieties with relatively short growth duration in the double-cropping region. Shoot dry weight (SDW) is one of the important traits associated with early growth vigor, and breeders have been working to improve this trait. Finding stable QTLs or functional genes for SDW is crucial for improving the early growth vigor by implementing molecular breeding in rice. Here, a genome-wide association analysis revealed that the QTL for SDW, qSDW-5, was stably detected in the three cultivation methods commonly used in production practice. Through gene-based haplotype analysis of the annotated genes within the putative region of qSDW-5, and validated by gene expression and knockout transgenic experiments, LOC_Os05g09520, which is identical to the reported GW5/GSE5 controlling grain width (GW) and thousand grain weight (TGW) was identified as the causal gene for qSDW-5. Five main haplotypes of LOC_Os05g09520 were identified in the diverse international rice collection used in this study and their effects on SDW, GW and TGW were analyzed. Phenotypic comparisons of the major haplotypes of LOC_Os05g09520 in the three subpopulations (indica, japonica and aus) revealed the same patterns of wider GW and higher TGW along with higher SDW. Furtherly, the haplotype analysis of 138 rice varieties/lines widely used in southern China showed that 97.8% of the cultivars/lines carry Hap2LOC_Os05g09520. These results not only provide a promising gene source for the molecular breeding of rice varieties with strong early growth vigor, but also elucidate the effect of the LOC_Os05g09520 haplotypes on SDW, GW, and TGW in rice. Importantly, this study provides direct genetic evidence that these three traits are significantly correlated, and suggests a breeding strategy for developing high-yielding and slender grain-shaped indica cultivars with strong early growth vigor.

Genome-wide association study provided insights into the polled phenotype and polled intersex syndrome (PIS) in goats

BACKGROUND

Breeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals.

METHODS

The present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats.

RESULTS

The analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats.

CONCLUSION

We have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits.

GWAS and Transcriptomic Analysis Identify OsRING315 as a New Candidate Gene Controlling Amylose Content and Gel Consistency in Rice

Cooking quality is the main factor determining the market value of rice. Although several major genes and a certain number of QTLs controlling cooking quality have been identified, the genetic complexity and environmental susceptibility limit the further improvement for cooking quality by molecular breeding. This research conducted a genome-wide association study to elucidate the QTLs related to cooking quality including amylose content (AC), gel consistency (GC) and alkali spreading value (ASV) by using 450 rice accessions consisting of 300 indica and 150 japonica accessions in two distinct environments. A total of 54 QTLs were identified, including 25 QTLs for AC, 12 QTLs for GC and 17 QTLs for ASV. Among them, 10 QTLs were consistently observed by the same population in both environments. Six QTLs were co-localized with the reported QTLs or cloned genes. The Wx gene for AC and GC, and the ALK gene for ASV were identified in every population across the two environments. The qAC9-2 for AC and the qGC9-2 for GC were defined to the same interval. The OsRING315 gene, encoding an E3 ubiquitin ligase, was considered as the candidate gene for both qAC9-2 and qGC9-2. The higher expression of OsRING315 corresponded to the lower AC and higher GC. Three haplotypes of OsRING315 were identified. The Hap 1 mainly existed in the japonica accessions and had lower AC. The Hap 2 and Hap 3 were predominantly present in the indica accessions, associated with higher AC. Meanwhile, the GC of accessions harboring Hap 1 was higher than that of accessions harboring Hap 3. In addition, the distribution of the three haplotypes in several rice-growing regions was unbalanced. The three traits of cooking quality are controlled by both major and minor genes and susceptible to environmental factors. The expression level of OsRING315 is related to both AC and GC, and this gene can be a promising target in quality improvement by using the gene editing method. Moreover, the haplotypes of OsRING315 differentiate between indica and japonica, and reveal the differences in GC and AC between indica and japonica rice.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Progress in Rice Breeding Based on Genomic Research

The role of rice genomics in breeding progress is becoming increasingly important. Deeper research into the rice genome will contribute to the identification and utilization of outstanding functional genes, enriching the diversity and genetic basis of breeding materials and meeting the diverse demands for various improvements. Here, we review the significant contributions of rice genomics research to breeding progress over the last 25 years, discussing the profound impact of genomics on rice genome sequencing, functional gene exploration, and novel breeding methods, and we provide valuable insights for future research and breeding practices.

PPanG: a precision pangenome browser enabling nucleotide-level analysis of genomic variations in individual genomes and their graph-based pangenome

Graph-based pangenome is gaining more popularity than linear pangenome because it stores more comprehensive information of variations. However, traditional linear genome browser has its own advantages, especially the tremendous resources accumulated historically. With the fast-growing number of individual genomes and their annotations available, the demand for a genome browser to visualize genome annotation for many individuals together with a graph-based pangenome is getting higher and higher. Here we report a new pangenome browser PPanG, a precise pangenome browser enabling nucleotide-level comparison of individual genome annotations together with a graph-based pangenome. Nine rice genomes with annotations were provided by default as potential references, and any individual genome can be selected as the reference. Our pangenome browser provides unprecedented insights on genome variations at different levels from base to gene, and reveals how the structures of a gene could differ for individuals. PPanG can be applied to any species with multiple individual genomes available and it is available at https://cgm.sjtu.edu.cn/PPanG .

Pangenome-Wide Association Study and Transcriptome Analysis Reveal a Novel QTL and Candidate Genes Controlling both Panicle and Leaf Blast Resistance in Rice

Cultivating rice varieties with robust blast resistance is the most effective and economical way to manage the rice blast disease. However, rice blast disease comprises leaf and panicle blast, which are different in terms of resistance mechanisms. While many blast resistant rice cultivars were bred using genes conferring resistance to only leaf or panicle blast, mining durable and effective quantitative trait loci (QTLs) for both panicle and leaf blast resistance is of paramount importance. In this study, we conducted a pangenome-wide association study (panGWAS) on 9 blast resistance related phenotypes using 414 international diverse rice accessions from an international rice panel. This approach led to the identification of 74 QTLs associated with rice blast resistance. One notable locus, qPBR1, validated in a F4:5 population and fine-mapped in a Heterogeneous Inbred Family (HIF), exhibited broad-spectrum, major and durable blast resistance throughout the growth period. Furthermore, we performed transcriptomic analysis of 3 resistant and 3 sensitive accessions at different time points after infection, revealing 3,311 differentially expressed genes (DEGs) potentially involved in blast resistance. Integration of the above results identified 6 candidate genes within the qPBR1 locus, with no significant negative effect on yield. The results of this study provide valuable germplasm resources, QTLs, blast response genes and candidate functional genes for developing rice varieties with enduring and broad-spectrum blast resistance. The qPBR1, in particular, holds significant potential for breeding new rice varieties with comprehensive and durable resistance throughout their growth period.

Technological Development and Advances for Constructing and Analyzing Plant Pangenomes

A pangenome captures the genomic diversity for a species, derived from a collection of genetic sequences of diverse populations. Advances in sequencing technologies have given rise to three primary methods for pangenome construction and analysis: de novo assembly and comparison, reference genome-based iterative assembly, and graph-based pangenome construction. Each method presents advantages and challenges in processing varying amounts and structures of DNA sequencing data. With the emergence of high-quality genome assemblies and advanced bioinformatic tools, the graph-based pangenome is emerging as an advanced reference for exploring the biological and functional implications of genetic variations.

Haplotype-resolved genome of Prunus zhengheensis provides insight into its evolution and low temperature adaptation in apricot

Prunus zhengheensis, an extremely rare population of apricots, originated in warm South-East China and is an excellent material for genetic breeding. However, most apricots and two related species (P. sibirica, P. mandshurica) are found in the cold northern regions in China and the mechanism of their distribution is still unclear. In addition, the classification status of P. zhengheensis is controversial. Thus, we generated a high-quality haplotype-resolved genome for P. zhengheensis, exploring key genetic variations in its adaptation and the causes of phylogenetic incongruence. We found extensive phylogenetic discordances between the nuclear and organelle phylogenies of P. zhengheensis, which could be explained by incomplete lineage sorting. A 242.22-Mb pan-genome of the Armeniaca section was developed with 13 chromosomal genomes. Importantly, we identified a 566-bp insertion in the promoter of the HSFA1d gene in apricot and showed that the activity of the HSFA1d promoter increased under low temperatures. In addition, HSFA1d overexpression in Arabidopsis thaliana indicated that HSFA1d positively regulated plant growth under chilling. Therefore, we hypothesized that the insertion in the promoter of HSFA1d in apricot improved its low-temperature adaptation, allowing it to thrive in relatively cold locations. The findings help explain the weather adaptability of Armeniaca plants.

The Function of SD1 on Shoot Length and its Pyramiding Effect on Shoot Length and Plant Height in Rice (Oryza sativa L.)

Strong seedling vigor is imperative to achieve stable seedling establishment and enhance the competitiveness against weeds in rice direct seeding. Shoot length (SL) is one of the important traits associated with seedling vigor in rice, but few genes for SL have been cloned so far. In the previous study, we identified two tightly linked and stably expressed QTLs for SL, qSL-1f and qSL-1d by genome-wide association study, and cloned the causal gene (LOC_Os01g68500) underlying qSL-1f. In the present study, we identify LOC_Os01g66100 (i.e. the semidwarf gene SD1), a well-known gene controlling plant height (PH) at the adult-plant stage, as the causal gene underlying qSL-1d through gene-based haplotype analysis and knockout transgenic verification. By measuring the phenotypes (SL and PH) of various haplotypes of the two genes and their knockout lines, we found SD1 and LOC_ Os01g68500 controlled both SL and PH, and worked in the same direction, which provided the directly genetic evidence for a positive correlation between SL and PH combined with the analysis of SL and PH in the diverse natural population. Moreover, the knockout transgenic experiments suggested that SD1 had a greater effect on PH compared with LOC_ Os01g68500, but no significant difference in the effect on SL. Further investigation of the pyramiding effects of SD1 and LOC_Os01g68500 based on their haplotype combinations suggested that SD1 may play a dominant role in controlling SL and PH when the two genes coexist. In this study, the effect of SD1 on SL at the seedling stage is validated. In total, two causal genes, SD1 and LOC_ Os01g68500, for SL are cloned in our studies, which controlled both SL and PH, and the suitable haplotypes of SD1 and LOC_ Os01g68500 are beneficial to achieve the desired SL and PH in different rice breeding objectives. These results provide a new clue to develop rice varieties for direct seeding and provide new genetic resources for molecular breeding of rice with suitable PH and strong seedling vigor.

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.

Genome-Wide Association Analysis Unravels New Quantitative Trait Loci (QTLs) for Eight Lodging Resistance Constituent Traits in Rice (Oryza sativa L.)

Lodging poses a significant challenge to rice yield, prompting the need to identify elite alleles for lodging resistance traits to improve cultivated rice varieties. In this study, a natural population of 518 rice accessions was examined to identify elite alleles associated with plant height (PH), stem diameter (SD), stem anti-thrust (AT/S), and various internode lengths (first (FirINL), second (SecINL), third (ThirINL), fourth (ForINL), and fifth (FifINL) internode lengths). A total of 262 SSR markers linked to these traits were uncovered through association mapping in two environmental conditions. Phenotypic evaluations revealed striking differences among cultivars, and genetic diversity assessments showed polymorphisms across the accessions. Favorable alleles were identified for PH, SD, AT/S, and one to five internode lengths, with specific alleles displaying considerable effects. Noteworthy alleles include RM6811-160 bp on chromosome 6 (which reduces PH) and RM161-145 bp on chromosome 5 (which increases SD). The study identified a total of 42 novel QTLs. Specifically, seven QTLs were identified for PH, four for SD, five for AT/S, five for FirINL, six for SecINL, five for ThirINL, six for ForINL, and four for FifINL. QTLs qAT/S-2, qPH2.1, qForINL2.1, and qFifINL exhibited the most significant phenotypic variance (PVE) of 3.99% for the stem lodging trait. AT/S, PH, ForINL, and FifINL had additive effects of 5.31 kPa, 5.42 cm, 4.27 cm, and 4.27 cm, respectively, offering insights into eight distinct cross-combinations for enhancing each trait. This research suggests the potential for crossbreeding superior parents based on stacked alleles, promising improved rice cultivars with enhanced lodging resistance to meet market demands.

Association of a Specific OsCULLIN3c Haplotype with Salt Stress Responses in Local Thai Rice

We previously found that OsCUL3c is involved in the salt stress response. However, there are no definitive reports on the diversity of OsCUL3c in local Thai rice. In this study, we showed that the CUL3 group was clearly separated from the other CUL groups; next, we focused on OsCUL3c, the third CUL3 of the CUL3 family in rice, which is absent in Arabidopsis. A total of 111 SNPs and 28 indels over the OsCUL3c region, representing 79 haplotypes (haps), were found. Haplotyping revealed that group I (hap A and hap C) and group II (hap B1 and hap D) were different mutated variants, which showed their association with phenotypes under salt stress. These results were supported by cis-regulatory elements (CREs) and transcription factor binding sites (TFBSs) analyses. We found that LTR, MYC, [AP2; ERF], and NF-YB, which are related to salt stress, drought stress, and the response to abscisic acid (ABA), have distinct positions and numbers in the haplotypes of group I and group II. An RNA Seq analysis of the two predominant haplotypes from each group showed that the OsCUL3c expression of the group I representative was upregulated and that of group II was downregulated, which was confirmed by RT-qPCR. Promoter changes might affect the transcriptional responses to salt stress, leading to different regulatory mechanisms for the expression of different haplotypes. We speculate that OsCUL3c influences the regulation of salt-related responses, and haplotype variations play a role in this regulation.

Genome-wide association mapping combined with gene-based haplotype analysis identify a novel gene for shoot length in rice (Oryza sativa L.)

KEY MESSAGE

Genome-wide association mapping revealed a novel QTL for shoot length across multiple environments. Its causal gene, LOC_Os01g68500, was identified firstly through gene-based haplotype analysis, gene expression and knockout transgenic verification. Strong seedling vigor is an important breeding target for rice varieties used in direct seeding. Shoot length (SL) is one of the important traits associated with seedling vigor characterized by rapid growth of seedling, which enhance seedling emergence. Therefore, mining genes for SL and conducting molecular breeding help to develop varieties for direct seeding. However, few QTLs for SL have been fine mapped or cloned so far. In this study, a genome-wide association study of SL was performed in a diverse rice collection consisting of 391 accessions in two years, using phenotypes generated by different cultivation methods according to the production practice, and a total of twenty-four QTLs for SL were identified. Among them, the novel QTL qSL-1f on chromosome 1 could be stably detected across all three cultivation methods in the whole population and indica subpopulation. Through gene-based haplotype analysis of the annotated genes within the putative region of qSL-1f, and validated by gene expression and knockout transgenic experiments, LOC_Os01g68500 (i.e., Os01g0913100 in RAP-DB) was identified as the causal gene for SL, which has a single-base variation (C-to-A transversion) in its CDS region, resulting in the significant difference in SL of rice. LOC_Os01g68500 encodes a DUF538 (Domain of unknown function) containing protein, and the function of DUF538 protein gene on rice seedling growth is firstly reported in this study. These results provide a new clue for exploring the molecular mechanism regulating SL, and promising gene source for the molecular breeding in rice.

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.

Integrated Genomic Selection for Accelerating Breeding Programs of Climate-Smart Cereals

Rapidly rising population and climate changes are two critical issues that require immediate action to achieve sustainable development goals. The rising population is posing increased demand for food, thereby pushing for an acceleration in agricultural production. Furthermore, increased anthropogenic activities have resulted in environmental pollution such as water pollution and soil degradation as well as alterations in the composition and concentration of environmental gases. These changes are affecting not only biodiversity loss but also affecting the physio-biochemical processes of crop plants, resulting in a stress-induced decline in crop yield. To overcome such problems and ensure the supply of food material, consistent efforts are being made to develop strategies and techniques to increase crop yield and to enhance tolerance toward climate-induced stress. Plant breeding evolved after domestication and initially remained dependent on phenotype-based selection for crop improvement. But it has grown through cytological and biochemical methods, and the newer contemporary methods are based on DNA-marker-based strategies that help in the selection of agronomically useful traits. These are now supported by high-end molecular biology tools like PCR, high-throughput genotyping and phenotyping, data from crop morpho-physiology, statistical tools, bioinformatics, and machine learning. After establishing its worth in animal breeding, genomic selection (GS), an improved variant of marker-assisted selection (MAS), has made its way into crop-breeding programs as a powerful selection tool. To develop novel breeding programs as well as innovative marker-based models for genetic evaluation, GS makes use of molecular genetic markers. GS can amend complex traits like yield as well as shorten the breeding period, making it advantageous over pedigree breeding and marker-assisted selection (MAS). It reduces the time and resources that are required for plant breeding while allowing for an increased genetic gain of complex attributes. It has been taken to new heights by integrating innovative and advanced technologies such as speed breeding, machine learning, and environmental/weather data to further harness the GS potential, an approach known as integrated genomic selection (IGS). This review highlights the IGS strategies, procedures, integrated approaches, and associated emerging issues, with a special emphasis on cereal crops. In this domain, efforts have been taken to highlight the potential of this cutting-edge innovation to develop climate-smart crops that can endure abiotic stresses with the motive of keeping production and quality at par with the global food demand.

Genome-wide association mapping and gene expression analysis reveal candidate genes for grain chalkiness in rice

Grain chalkiness is the main factor determining the market value of rice. Reducing chalkiness is an important breeding goal for genetic improvement of high quality rice. Identification of QTLs or genes controlling chalkiness is the prerequisite for molecular breeding in rice. Here, we conducted a genome-wide association study to identify QTLs associated with grain chalkiness including percentage of grains with chalkiness (PGWC) and degree of endosperm chalkiness (DEC) in 450 rice accessions consisting of 300 indica and 150 japonica rice in two environments. A total of 34 QTLs were identified, including 14 QTLs for PGWC and 20 QTLs for DEC. Among them, seven QTLs were commonly identified in two environments, and eight QTLs were simultaneously related to two traits. Based on the haplotype analysis, LD decay analysis, RNA-sequencing, qRT-PCR confirmation and haplotype comparisons, four genes (LOC_Os10g36170, LOC_Os10g36260, LOC_Os10g36340 and LOC_Os10g36610) were considered as the candidate genes for qDEC-10c^1w,2wj , which could be identified in both environments and had the most significant p-value among the newly identified QTLs. These results provided new insight into the genetic basis of grain chalkiness and gene resources for improving quality by molecular breeding in rice.