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

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.

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.

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.

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.

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.