Structural variations in plant genomes

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.

Pangenome characterization and analysis of the NAC gene family reveals genes for Sclerotinia sclerotiorum resistance in sunflower (Helianthus annuus)

BACKGROUND

Sunflower (Helianthus annuus) is one of the most important economic crops in oilseed production worldwide. The different cultivars exhibit variability in their resistance genes. The NAC transcription factor (TF) family plays diverse roles in plant development and stress responses. With the completion of the H. annuus genome sequence, the entire complement of genes coding for NACs has been identified. However, the reference genome of a single individual cannot cover all the genetic information of the species.

RESULTS

Considering only a single reference genome to study gene families will miss many meaningful genes. A pangenome-wide survey and characterization of the NAC genes in sunflower species were conducted. In total, 139 HaNAC genes are identified, of which 114 are core and 25 are variable. Phylogenetic analysis of sunflower NAC proteins categorizes these proteins into 16 subgroups. 138 HaNACs are randomly distributed on 17 chromosomes. SNP-based haplotype analysis shows haplotype diversity of the HaNAC genes in wild accessions is richer than in landraces and modern cultivars. Ten HaNAC genes in the basal stalk rot (BSR) resistance quantitative trait loci (QTL) are found. A total of 26 HaNAC genes are differentially expressed in response to Sclerotinia head rot (SHR). A total of 137 HaNAC genes are annotated in Gene Ontology (GO) and are classified into 24 functional groups. GO functional enrichment analysis reveals that HaNAC genes are involved in various functions of the biological process.

CONCLUSIONS

We identified NAC genes in H. annuus (HaNAC) on a pangenome-wide scale and analyzed S. sclerotiorum resistance-related NACs. This study provided a theoretical basis for further genomic improvement targeting resistance-related NAC genes in sunflowers.

Recent advances in plant translational genomics for crop improvement

The growing population, climate change, and limited agricultural resources put enormous pressure on agricultural systems. A plateau in crop yields is occurring and extreme weather events and urbanization threaten the livelihood of farmers. It is imperative that immediate attention is paid to addressing the increasing food demand, ensuring resilience against emerging threats, and meeting the demand for more nutritious, safer food. Under uncertain conditions, it is essential to expand genetic diversity and discover novel crop varieties or variations to develop higher and more stable yields. Genomics plays a significant role in developing abundant and nutrient-dense food crops. An alternative to traditional breeding approach, translational genomics is able to improve breeding programs in a more efficient and precise manner by translating genomic concepts into practical tools. Crop breeding based on genomics offers potential solutions to overcome the limitations of conventional breeding methods, including improved crop varieties that provide more nutritional value and are protected from biotic and abiotic stresses. Genetic markers, such as SNPs and ESTs, contribute to the discovery of QTLs controlling agronomic traits and stress tolerance. In order to meet the growing demand for food, there is a need to incorporate QTLs into breeding programs using marker-assisted selection/breeding and transgenic technologies. This chapter primarily focuses on the recent advances that are made in translational genomics for crop improvement and various omics techniques including transcriptomics, metagenomics, pangenomics, single cell omics etc. Numerous genome editing techniques including CRISPR Cas technology and their applications in crop improvement had been discussed.

Effect of structural variation in the promoter region of RsMYB1.1 on the skin color of radish taproot

Accumulation of anthocyanins in the taproot of radish is an agronomic trait beneficial for human health. Several genetic loci are related to a red skin or flesh color of radish, however, the functional divergence of candidate genes between non-red and red radishes has not been investigated. Here, we report that a novel genetic locus on the R2 chromosome, where RsMYB1.1 is located, is associated with the red color of the skin of radish taproot. A genome-wide association study (GWAS) of 66 non-red-skinned (nR) and 34 red-skinned (R) radish accessions identified three nonsynonymous single nucleotide polymorphisms (SNPs) in the third exon of RsMYB1.1. Although the genotypes of SNP loci differed between the nR and R radishes, no functional difference in the RsMYB1.1 proteins of nR and R radishes in their physical interaction with RsTT8 was detected by yeast-two hybrid assay or in anthocyanin accumulation in tobacco and radish leaves coexpressing RsMYB1.1 and RsTT8. By contrast, insertion- or deletion-based GWAS revealed that one large AT-rich low-complexity sequence of 1.3-2 kb was inserted in the promoter region of RsMYB1.1 in the nR radishes (RsMYB1.1nR), whereas the R radishes had no such insertion; this represents a presence/absence variation (PAV). This insertion sequence (RsIS) was radish specific and distributed among the nine chromosomes of Raphanus genomes. Despite the extremely low transcription level of RsMYB1.1nR in the nR radishes, the inactive RsMYB1.1nR promoter could be functionally restored by deletion of the RsIS. The results of a transient expression assay using radish root sections suggested that the RsIS negatively regulates the expression of RsMYB1.1nR, resulting in the downregulation of anthocyanin biosynthesis genes, including RsCHS, RsDFR, and RsANS, in the nR radishes. This work provides the first evidence of the involvement of PAV in an agronomic trait of radish.

Pan-genome analysis of 13 Malus accessions reveals structural and sequence variations associated with fruit traits

Structural variations (SVs) and copy number variations (CNVs) contribute to trait variations in fleshy-fruited species. Here, we assemble 10 genomes of genetically diverse Malus accessions, including the ever-green cultivar 'Granny Smith' and the widely cultivated cultivar 'Red Fuji'. Combining with three previously reported genomes, we assemble the pan-genome of Malus species and identify 20,220 CNVs and 317,393 SVs. We also observe CNVs that are positively correlated with expression levels of the genes they are associated with. Furthermore, we show that the noncoding RNA generated from a 209 bp insertion in the intron of mitogen-activated protein kinase homology encoding gene, MMK2, regulates the gene expression and affects fruit coloration. Moreover, we identify overlapping SVs associated with fruit quality and biotic resistance. This pan-genome uncovers possible contributions of CNVs to gene expression and highlights the role of SVs in apple domestication and economically important traits.

A comprehensive lettuce variation map reveals the impact of structural variations in agronomic traits

BACKGROUND

As an important vegetable crop, cultivated lettuce is grown worldwide and a great variety of agronomic traits have been preserved within germplasm collections. The mechanisms underlying these phenotypic variations remain to be elucidated in association with sequence variations. Compared with single nucleotide polymorphisms, structural variations (SVs) that have more impacts on gene functions remain largely uncharacterized in the lettuce genome.

RESULTS

Here, we produced a comprehensive SV set for 333 wild and cultivated lettuce accessions. Comparison of SV frequencies showed that the SVs prevalent in L. sativa affected the genes enriched in carbohydrate derivative catabolic and secondary metabolic processes. Genome-wide association analysis of seven agronomic traits uncovered potentially causal SVs associated with seed coat color and leaf anthocyanin content.

CONCLUSION

Our work characterized a great abundance of SVs in the lettuce genome, and provides a valuable genomic resource for future lettuce breeding.

A fast and globally optimal solution for RNA-seq quantification

Alignment-based RNA-seq quantification methods typically involve a time-consuming alignment process prior to estimating transcript abundances. In contrast, alignment-free RNA-seq quantification methods bypass this step, resulting in significant speed improvements. Existing alignment-free methods rely on the Expectation-Maximization (EM) algorithm for estimating transcript abundances. However, EM algorithms only guarantee locally optimal solutions, leaving room for further accuracy improvement by finding a globally optimal solution. In this study, we present TQSLE, the first alignment-free RNA-seq quantification method that provides a globally optimal solution for transcript abundances estimation. TQSLE adopts a two-step approach: first, it constructs a k-mer frequency matrix A for the reference transcriptome and a k-mer frequency vector b for the RNA-seq reads; then, it directly estimates transcript abundances by solving the linear equation ATAx = ATb. We evaluated the performance of TQSLE using simulated and real RNA-seq data sets and observed that, despite comparable speed to other alignment-free methods, TQSLE outperforms them in terms of accuracy. TQSLE is freely available at https://github.com/yhg926/TQSLE.

Evolutionary roles of polyploidization-derived structural variations in the phenotypic diversification of Panax species

Genomic structural variations (SVs) are widespread in plant and animal genomes and play important roles in phenotypic novelty and species adaptation. Frequent whole genome duplications followed by (re)diploidizations have resulted in high diversity of genome architecture among extant species. In this study, we identified abundant genomic SVs in the Panax genus that are hypothesized to have occurred through during the repeated polyploidizations/(re)diploidizations. Our genome-wide comparisons demonstrated that although these polyploidization-derived SVs have evolved at distinct evolutionary stages, a large number of SV-intersecting genes showed enrichment in functionally important pathways related to secondary metabolites, photosynthesis and basic cellular activities. In line with these observations, our metabolic analyses of these Panax species revealed high diversity of primary and secondary metabolites both at the tissue and interspecific levels. In particular, genomic SVs identified at ginsenoside biosynthesis genes, including copy number variation and large fragment deletion, appear to have played important roles in the evolution and diversification of ginsenosides. A further herbivore deterrence experiment demonstrated that, as major triterpenoidal saponins found exclusively in Panax, ginsenosides provide protection against insect herbivores. Our study provides new insights on how polyploidization-derived SVs have contributed to phenotypic novelty and plant adaptation.

Two complementary genes in a presence-absence variation contribute to indica-japonica reproductive isolation in rice

Understanding the evolutionary forces in speciation is a central goal in evolutionary biology. Asian cultivated rice has two subspecies, indica and japonica, but the underlying mechanism of the partial reproductive isolation between them remains obscure. Here we show a presence-absence variation (PAV) at the Se locus functions as an indica-japonica reproductive barrier by causing hybrid sterility (HS) in indica-japonica crosses. The locus comprises two adjacent genes: ORF3 encodes a sporophytic pollen killer, whereas ORF4 protects pollen in a gametophytic manner. In F₁ of indica-japonica crosses, pollen with the japonica haplotype, which lacks the sequence containing the protective ORF4, is aborted due to the pollen-killing effect of ORF3 from indica. Evolutionary analysis suggests ORF3 is a gene associated with the Asian cultivated rice species complex, and the PAV has contributed to the reproductive isolation between the two subspecies of Asian cultivated rice. Our analyses provide perspectives on rice inter-subspecies post-zygotic isolation, and will promote efforts to overcome reproductive barriers in indica-japonica hybrid rice breeding.

Recent advancements in the breeding of sorghum crop: current status and future strategies for marker-assisted breeding

Sorghum is emerging as a model crop for functional genetics and genomics of tropical grasses with abundant uses, including food, feed, and fuel, among others. It is currently the fifth most significant primary cereal crop. Crops are subjected to various biotic and abiotic stresses, which negatively impact on agricultural production. Developing high-yielding, disease-resistant, and climate-resilient cultivars can be achieved through marker-assisted breeding. Such selection has considerably reduced the time to market new crop varieties adapted to challenging conditions. In the recent years, extensive knowledge was gained about genetic markers. We are providing an overview of current advances in sorghum breeding initiatives, with a special focus on early breeders who may not be familiar with DNA markers. Advancements in molecular plant breeding, genetics, genomics selection, and genome editing have contributed to a thorough understanding of DNA markers, provided various proofs of the genetic variety accessible in crop plants, and have substantially enhanced plant breeding technologies. Marker-assisted selection has accelerated and precised the plant breeding process, empowering plant breeders all around the world.

A comprehensive analysis of copy number variations in diverse apple populations

BACKGROUND

As an important source of genetic variation, copy number variation (CNV) can alter the dosage of DNA segments, which in turn may affect gene expression level and phenotype. However, our knowledge of CNV in apple is still limited. Here, we obtained high-confidence CNVs and investigated their functional impact based on genome resequencing data of two apple populations, cultivars and wild relatives.

RESULTS

In this study, we identified 914,610 CNVs comprising 14,839 CNV regions (CNVRs) from 346 apple accessions, including 289 cultivars and 57 wild relatives. CNVRs summed to 71.19 Mb, accounting for 10.03% of the apple genome. Under the low linkage disequilibrium (LD) with nearby SNPs, they could also accurately reflect the population structure of apple independent of SNPs. Furthermore, A total of 3,621 genes were covered by CNVRs and functionally involved in biological processes such as defense response, reproduction and metabolic processes. In addition, the population differentiation index ([Formula: see text]) analysis between cultivars and wild relatives revealed 127 CN-differentiated genes, which may contribute to trait differences in these two populations.

CONCLUSIONS

This study was based on identification of CNVs from 346 diverse apple accessions, which to our knowledge was the largest dataset for CNV analysis in apple. Our work presented the first comprehensive CNV map and provided valuable resources for understanding genomic variations in apple.

Exploring the sorghum race level diversity utilizing 272 sorghum accessions genomic resources

Due to evolutionary divergence, sorghum race populations exhibit significant genetic and morphological variation. A k-mer-based sorghum race sequence comparison identified the conserved k-mers of all 272 accessions from sorghum and the race-specific genetic signatures identified the gene variability in 10,321 genes (PAVs). To understand sorghum race structure, diversity and domestication, a deep learning-based variant calling approach was employed in a set of genotypic data derived from a diverse panel of 272 sorghum accessions. The data resulted in 1.7 million high-quality genome-wide SNPs and identified selective signature (both positive and negative) regions through a genome-wide scan with different (iHS and XP-EHH) statistical methods. We discovered 2,370 genes associated with selection signatures including 179 selective sweep regions distributed over 10 chromosomes. Co-localization of these regions undergoing selective pressure with previously reported QTLs and genes revealed that the signatures of selection could be related to the domestication of important agronomic traits such as biomass and plant height. The developed k-mer signatures will be useful in the future to identify the sorghum race and for trait and SNP markers for assisting in plant breeding programs.

Identification of structural variations related to drought tolerance in wheat (Triticum aestivum L.)

Structural variations are common in plant genomes, affecting meiotic recombination and distorted segregation in wheat. And presence/absence variations can significantly affect drought tolerance in wheat. Drought is a major abiotic stress limiting wheat production. Common wheat has a complex genome with three sub-genomes, which host large numbers of structural variations (SVs). SVs play critical roles in understanding the genetic contributions of plant domestication and phenotypic plasticity, but little is known about their genomic characteristics and their effects on drought tolerance. In the present study, high-resolution karyotypes of 180 doubled haploids (DHs) were developed. Signal polymorphisms between the parents involved with 8 presence-absence variations (PAVs) of tandem repeats (TR) distributed on the 7 (2A, 4A, 5A, 7A, 3B, 7B, and 2D) of 21 chromosomes. Among them, PAV on chromosome 2D showed distorted segregation, others transmit normal conforming to a 1:1 segregation ration in the population; and a PAVs recombination occurred on chromosome 2A. Association analysis of PAV and phenotypic traits under different water regimes, we found PAVs on chromosomes 4A, 5A, and 7B showed negative effect on grain length (GL) and grain width (GW); PAV.7A had opposite effect on grain thickness (GT) and spike length (SL), with the effect on traits differing under different water regimes. PAVs on linkage group 2A, 4A, 7A, 2D, and 7B associated with the drought tolerance coefficients (DTCs), and significant negative effect on drought resistance values (D values) were detected in PAV.7B. Additionally, quantitative trait loci (QTL) associated with phenotypic traits using the 90 K SNP array showed QTL for DTCs and grain-related traits in chromosomes 4A, and 5A, 3B were co-localized in differential regions of PAVs. These PAVs can cause the differentiation of the target region of SNP and could be used for genetic improvement of agronomic traits under drought stress via marker-assisted selection (MAS) breeding.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

Application of crop wild relatives in modern breeding: An overview of resources, experimental and computational methodologies

Global agricultural industries are under pressure to meet the future food demand; however, the existing crop genetic diversity might not be sufficient to meet this expectation. Advances in genome sequencing technologies and availability of reference genomes for over 300 plant species reveals the hidden genetic diversity in crop wild relatives (CWRs), which could have significant impacts in crop improvement. There are many ex-situ and in-situ resources around the world holding rare and valuable wild species, of which many carry agronomically important traits and it is crucial for users to be aware of their availability. Here we aim to explore the available ex-/in- situ resources such as genebanks, botanical gardens, national parks, conservation hotspots and inventories holding CWR accessions. In addition we highlight the advances in availability and use of CWR genomic resources, such as their contribution in pangenome construction and introducing novel genes into crops. We also discuss the potential and challenges of modern breeding experimental approaches (e.g. de novo domestication, genome editing and speed breeding) used in CWRs and the use of computational (e.g. machine learning) approaches that could speed up utilization of CWR species in breeding programs towards crop adaptability and yield improvement.

Chromosome-level genome assembly of Dongxiang wild rice (Oryza rufipogon) provides insights into resistance to disease and freezing

Dongxiang wild rice (DXWR, Oryza rufipogon Griff.) belongs to common wild rice O. rufipogon, which is the well-known ancestral progenitor of cultivated rice, possessing important gene resources for rice breeding. However, the distribution of DXWR is decreasing rapidly, and no reference genome has been published to date. In this study, we constructed a chromosome-level reference genome of DXWR by Oxford Nanopore Technology (ONT) and High-through chromosome conformation capture (Hi-C). A total of 58.41 Gb clean data from ONT were de novo assembled into 231 contigs with the total length of 413.46 Mb and N50 length of 5.18 Mb. These contigs were clustered and ordered into 12 pseudo-chromosomes covering about 97.39% assembly with Hi-C data, with a scaffold N50 length of 33.47 Mb. Moreover, 54.10% of the genome sequences were identified as repeat sequences. 33,862 (94.21%) genes were functionally annotated from a total of predicted 35,942 protein-coding sequences. Compared with other species of Oryza genus, the genes related to disease and cold resistance in DXWR had undergone a large-scale expansion, which may be one of the reasons for the stronger disease resistance and cold resistance of DXWR. Comparative transcriptome analysis also determined a list of differentially expressed genes under normal and cold treatment, which supported DXWR as a cold-tolerant variety. The collinearity between DXWR and cultivated rice was high, but there were still some significant structural variations, including a specific inversion on chromosome 11, which may be related to the differentiation of DXWR. The high-quality chromosome-level reference genome of DXWR assembled in this study will become a valuable resource for rice molecular breeding and genetic research in the future.

Legume Pangenome: Status and Scope for Crop Improvement

In the last decade, legume genomics research has seen a paradigm shift due to advances in genome sequencing technologies, assembly algorithms, and computational genomics that enabled the construction of high-quality reference genome assemblies of major legume crops. These advances have certainly facilitated the identification of novel genetic variants underlying the traits of agronomic importance in many legume crops. Furthermore, these robust sequencing technologies have allowed us to study structural variations across the whole genome in multiple individuals and at the species level using 'pangenome analysis.' This review updates the progress of constructing pangenome assemblies for various legume crops and discusses the prospects for these pangenomes and how to harness the information to improve various traits of economic importance through molecular breeding to increase genetic gain in legumes and tackle the increasing global food crisis.

Location of Tandem Repeats on Wheat Chromosome 5B and the Breakpoint on the 5BS Arm in Wheat Translocation T7BS.7BL-5BS Using Single-Copy FISH Analysis

Wheat (Triticum aestivum L.) is rich in tandem repeats, and this is helpful in studying its karyotypic evolution. Some tandem repeats have not been assembled into the wheat genome sequence. Alignment using the blastn tool in the B2DSC web server indicated that the genomic sequence of 5B chromosome (IWGSC RefSeq v2.1) does not contain the tandem repeat pTa-275, and the tandem repeat (GA)₂₆ distributed throughout the whole 5B chromosome. The nondenaturing fluorescence in situ hybridization (ND-FISH) using the oligonucleotide (oligo) probes derived from pTa-275 and (GA)₂₆ indicated that one signal band of pTa-275 and two signal bands of (GA)₂₆ appeared on the 5B chromosome of Chinese Spring wheat, indicating the aggregative distribution patterns of the two kinds of tandem repeats. Single-copy FISH indicated that the clustering region of pTa-275 and the two clustering regions of (GA)₂₆ were located in ~160-201 Mb, ~153-157 Mb, and ~201-234 Mb intervals, respectively. Using ND-FISH and single-copy FISH technologies, the translocation breakpoint on the 5BS portion of the translocation T7BS.7BL-5BS, which exists widely in north-western European wheat cultivars, was located in the region from 157,749,421 bp to 158,555,080 bp (~0.8 Mb), and this region mainly contains retrotransposons, and no gene was found. The clustering regions of two kinds of tandem repeats on wheat chromosome 5B were determined and this will be helpful to improve the future sequence assembly of this chromosome. The sequence characteristics of the translocation breakpoint on the translocation T7BS.7BL-5BS obtained in this study are helpful to understand the mechanism of wheat chromosome translocation.

Chromosome-Level Genome Assembly of a Fragrant Japonica Rice Cultivar 'Changxianggeng 1813' Provides Insights into Genomic Variations between Fragrant and Non-Fragrant Japonica Rice

East Asia has an abundant resource of fragrant japonica rice that is gaining increasing interest among both consumers and producers. However, genomic resources and in particular complete genome sequences currently available for the breeding of fragrant japonica rice are still scarce. Here, integrating Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods, we presented a high-quality chromosome-level genome assembly (~378.78 Mb) for a new fragrant japonica cultivar ‘Changxianggeng 1813’, with 31,671 predicated protein-coding genes. Based on the annotated genome sequence, we demonstrated that it was the badh2-E2 type of deletion (a 7-bp deletion in the second exon) that caused fragrance in ‘Changxianggeng 1813’. Comparative genomic analyses revealed that multiple gene families involved in the abiotic stress response were expanded in the ‘Changxianggeng 1813’ genome, which further supported the previous finding that no generalized loss of abiotic stress tolerance associated with the fragrance phenotype. Although the ‘Changxianggeng 1813’ genome showed high genomic synteny with the genome of the non-fragrant japonica rice cultivar Nipponbare, a total of 289,970 single nucleotide polymorphisms (SNPs), 96,093 small insertion-deletion polymorphisms (InDels), and 8690 large structure variants (SVs, >1000 bp) were identified between them. Together, these genomic resources will be valuable for elucidating the mechanisms underlying economically important traits and have wide-ranging implications for genomics-assisted breeding in fragrant japonica rice.

A telomere-to-telomere gap-free reference genome of watermelon and its mutation library provide important resources for gene discovery and breeding

Watermelon, Citrullus lanatus, is the world's third largest fruit crop. Reference genomes with gaps and a narrow genetic base hinder functional genomics and genetic improvement of watermelon. Here, we report the assembly of a telomere-to-telomere gap-free genome of the elite watermelon inbred line G42 by incorporating high-coverage and accurate long-read sequencing data with multiple assembly strategies. All 11 chromosomes have been assembled into single-contig pseudomolecules without gaps, representing the highest completeness and assembly quality to date. The G42 reference genome is 369 321 829 bp in length and contains 24 205 predicted protein-coding genes, with all 22 telomeres and 11 centromeres characterized. Furthermore, we established a pollen-EMS mutagenesis protocol and obtained over 200 000 M1 seeds from G42 . In a sampling pool, 48 monogenic phenotypic mutations, selected from 223 M1 and 78 M2 mutants with morphological changes, were confirmed. The average mutation density was 1 SNP/1.69 Mb and 1 indel/4.55 Mb per M1 plant and 1 SNP/1.08 Mb and 1 indel/6.25 Mb per M2 plant. Taking advantage of the gap-free G42 genome, 8039 mutations from 32 plants sampled from M1 and M2 families were identified with 100% accuracy, whereas only 25% of the randomly selected mutations identified using the 97103v2 reference genome could be confirmed. Using this library and the gap-free genome, two genes responsible for elongated fruit shape and male sterility (ClMS1) were identified, both caused by a single base change from G to A. The validated gap-free genome and its EMS mutation library provide invaluable resources for functional genomics and genetic improvement of watermelon.

Identification of Heat-Tolerant Genes in Non-Reference Sequences in Rice by Integrating Pan-Genome, Transcriptomics, and QTLs

The availability of large-scale genomic data resources makes it very convenient to mine and analyze genes that are related to important agricultural traits in rice. Pan-genomes have been constructed to provide insight into the genome diversity and functionality of different plants, which can be used in genome-assisted crop improvement. Thus, a pan-genome comprising all genetic elements is crucial for comprehensive variation study among the heat-resistant and -susceptible rice varieties. In this study, a rice pan-genome was firstly constructed by using 45 heat-tolerant and 15 heat-sensitive rice varieties. A total of 38,998 pan-genome genes were identified, including 37,859 genes in the reference and 1141 in the non-reference contigs. Genomic variation analysis demonstrated that a total of 76,435 SNPs were detected and identified as the heat-tolerance-related SNPs, which were specifically present in the highly heat-resistant rice cultivars and located in the genic regions or within 2 kbp upstream and downstream of the genes. Meanwhile, 3214 upregulated and 2212 downregulated genes with heat stress tolerance-related SNPs were detected in one or multiple RNA-seq datasets of rice under heat stress, among which 24 were located in the non-reference contigs of the rice pan-genome. We then mapped the DEGs with heat stress tolerance-related SNPs to the heat stress-resistant QTL regions. A total of 1677 DEGs, including 990 upregulated and 687 downregulated genes, were mapped to the 46 heat stress-resistant QTL regions, in which 2 upregulated genes with heat stress tolerance-related SNPs were identified in the non-reference sequences. This pan-genome resource is an important step towards the effective and efficient genetic improvement of heat stress resistance in rice to help meet the rapidly growing needs for improved rice productivity under different environmental stresses. These findings provide further insight into the functional validation of a number of non-reference genes and, especially, the two genes identified in the heat stress-resistant QTLs in rice.

Fine mapping of Pm58 from Aegilops tauschii conferring powdery mildew resistance

The powdery mildew resistance gene Pm58 was traced to a 141.3-kb interval with the co-segregating marker Xkasp68500 in wheat breeding. Pm58 is a powdery mildew resistance gene identified in Aegilops tauschii accession TA1662 and effective in a common wheat background. To finely map Pm58, an F₂ population of 676 plants derived from the cross T093 × TA1662 was used for recombinant screening. We obtained 13 recombinants that occurred between the flanking markers Xhnu670 and Xhnu186. Genotyping and phenotyping these recombinant F_2:3 families delimited Pm58 to a 0.22-cM interval (Xsts20220-Xkasp61553) on chromosome arm 2DS. The region carrying the Pm58 locus was approximately 141.3-kb, which contained eight annotated genes according to the reference genome sequence of Ae. tauschii AL8/78. Haplotype analysis of 178 Ae. tauschii accessions using the candidate gene-specific markers identified a disease resistance gene AET2Gv20068500 as a candidate for Pm58. Comparative mapping of the Pm58-containing interval revealed two presence/absence variations (PAVs) between AL8/78 and common wheat Chinese Spring. PAV-1 resides in the 3'-end of AET2Gv20068500. The majority of 158 common wheat cultivars (84.8%) displayed the absence of a 14.1-kb fragment in the PAV-1 region, which was confirmed by aligning the targeted genome sequences of the other sequenced Ae. tauschii accessions and common wheat cultivars. A co-segregating marker Xkasp68500 developed from AET2Gv20068500 can distinguish TA1662 from all randomly selected common wheat cultivars and will be instrumental for tracking Pm58 in breeding programs.

Oxford Nanopore and Bionano Genomics technologies evaluation for plant structural variation detection

Background

Structural Variations (SVs) are genomic rearrangements derived from duplication, deletion, insertion, inversion, and translocation events. In the past, SVs detection was limited to cytological approaches, then to Next-Generation Sequencing (NGS) short reads and partitioned assemblies. Nowadays, technologies such as DNA long read sequencing and optical mapping have revolutionized the understanding of SVs in genomes, due to the enhancement of the power of SVs detection.

This study aims to investigate performance of two techniques, 1) long-read sequencing obtained with the MinION device (Oxford Nanopore Technologies) and 2) optical mapping obtained with Saphyr device (Bionano Genomics) to detect and characterize SVs in the genomes of the two ecotypes of Arabidopsis thaliana, Columbia-0 (Col-0) and Landsberg erecta 1 (Ler-1).

Results

We described the SVs detected from the alignment of the best ONT assembly and DLE-1 optical maps of A. thaliana Ler-1 against the public reference genome Col-0 TAIR10.1. After filtering (SV > 1 kb), 1184 and 591 Ler-1 SVs were retained from ONT and Bionano technologies respectively. A total of 948 Ler-1 ONT SVs (80.1%) corresponded to 563 Bionano SVs (95.3%) leading to 563 common locations. The specific locations were scrutinized to assess improvement in SV detection by either technology. The ONT SVs were mostly detected near TE and gene features, and resistance genes seemed particularly impacted.

Conclusions

Structural variations linked to ONT sequencing error were removed and false positives limited, with high quality Bionano SVs being conserved. When compared with the Col-0 TAIR10.1 reference genome, most of the detected SVs discovered by both technologies were found in the same locations. ONT assembly sequence leads to more specific SVs than Bionano one, the latter being more efficient to characterize large SVs. Even if both technologies are complementary approaches, ONT data appears to be more adapted to large scale populations studies, while Bionano performs better in improving assembly and describing specificity of a genome compared to a reference.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08499-4.

Genomic prediction of the recombination rate variation in barley - A route to highly recombinogenic genotypes

Meiotic recombination is not only fundamental to the adaptation of sexually reproducing eukaryotes in nature but increased recombination rates facilitate the combination of favourable alleles into a single haplotype in breeding programmes. The main objectives of this study were to (i) assess the extent and distribution of the recombination rate variation in cultivated barley (Hordeum vulgare L.), (ii) quantify the importance of the general and specific recombination effects, and (iii) evaluate a genomic selection approach's ability to predict the recombination rate variation. Genetic maps were created for the 45 segregating populations that were derived from crosses among 23 spring barley inbreds with origins across the world. The genome-wide recombination rate among populations ranged from 0.31 to 0.73 cM/Mbp. The crossing design used in this study allowed to separate the general recombination effects (GRE) of individual parental inbreds from the specific recombination effects (SRE) caused by the combinations of parental inbreds. The variance of the genome-wide GRE was found to be about eight times the variance of the SRE. This finding indicated that parental inbreds differ in the efficiency of their recombination machinery. The ability to predict the chromosome or genome-wide recombination rate of an inbred ranged from 0.80 to 0.85. These results suggest that a reliable screening of large genetic materials for their potential to cause a high extent of genetic recombination in their progeny is possible, allowing to systematically manipulate the recombination rate using natural variation.

Pangenomics in Microbial and Crop Research: Progress, Applications, and Perspectives

Advances in sequencing technologies and bioinformatics tools have fueled a renewed interest in whole genome sequencing efforts in many organisms. The growing availability of multiple genome sequences has advanced our understanding of the within-species diversity, in the form of a pangenome. Pangenomics has opened new avenues for future research such as allowing dissection of complex molecular mechanisms and increased confidence in genome mapping. To comprehensively capture the genetic diversity for improving plant performance, the pangenome concept is further extended from species to genus level by the inclusion of wild species, constituting a super-pangenome. Characterization of pangenome has implications for both basic and applied research. The concept of pangenome has transformed the way biological questions are addressed. From understanding evolution and adaptation to elucidating host–pathogen interactions, finding novel genes or breeding targets to aid crop improvement to design effective vaccines for human prophylaxis, the increasing availability of the pangenome has revolutionized several aspects of biological research. The future availability of high-resolution pangenomes based on reference-level near-complete genome assemblies would greatly improve our ability to address complex biological problems.

Pangenomics in crop improvement-from coding structural variations to finding regulatory variants with pangenome graphs

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

Combined use of Oxford Nanopore and Illumina sequencing yields insights into soybean structural variation biology

BACKGROUND

Structural variants (SVs), including deletions, insertions, duplications, and inversions, are relatively long genomic variations implicated in a diverse range of processes from human disease to ecology and evolution. Given their complex signatures, tendency to occur in repeated regions, and large size, discovering SVs based on short reads is challenging compared to single-nucleotide variants. The increasing availability of long-read technologies has greatly facilitated SV discovery; however, these technologies remain too costly to apply routinely to population-level studies. Here, we combined short-read and long-read sequencing technologies to provide a comprehensive population-scale assessment of structural variation in a panel of Canadian soybean cultivars.

RESULTS

We used Oxford Nanopore long-read sequencing data (~12× mean coverage) for 17 samples to both benchmark SV calls made from Illumina short-read data and predict SVs that were subsequently genotyped in a population of 102 samples using Illumina data. Benchmarking results show that variants discovered using Oxford Nanopore can be accurately genotyped from the Illumina data. We first use the genotyped deletions and insertions for population genetics analyses and show that results are comparable to those based on single-nucleotide variants. We observe that the population frequency and distribution within the genome of deletions and insertions are constrained by the location of genes. Gene Ontology and PFAM domain enrichment analyses also confirm previous reports that genes harboring high-frequency deletions and insertions are enriched for functions in defense response. Finally, we discover polymorphic transposable elements from the deletions and insertions and report evidence of the recent activity of a Stowaway MITE.

CONCLUSIONS

We show that structural variants discovered using Oxford Nanopore data can be genotyped with high accuracy from Illumina data. Our results demonstrate that long-read and short-read sequencing technologies can be efficiently combined to enhance SV analysis in large populations, providing a reusable framework for their study in a wider range of samples and non-model species.

Comparison of structural variants in the whole genome sequences of two Medicago truncatula ecotypes: Jemalong A17 and R108

BACKGROUND

Structural variants (SVs) constitute a large proportion of the genomic variation that results in phenotypic variation in plants. However, they are still a largely unexplored feature in most plant genomes. Here, we present the whole-genome landscape of SVs between two model legume Medicago truncatula ecotypes-Jemalong A17 and R108- that have been extensively used in various legume biology studies.

RESULTS

To catalogue SVs, we first resolved the previously published R108 genome assembly (R108 v1.0) to chromosome-scale using 124 × Hi-C data, resulting in a high-quality genome assembly. The inter-chromosomal reciprocal translocations between chromosomes 4 and 8 were confirmed by performing syntenic analysis between the two genomes. Combined with the Hi-C data, it appears that these translocation events had a significant effect on chromatin organization. Using both whole-genome and short-read alignments, we identified the genomic landscape of SVs between the two genomes, some of which may account for several phenotypic differences, including their differential responses to aluminum toxicity and iron deficiency, and the development of different anthocyanin leaf markings. We also found extensive SVs within the nodule-specific cysteine-rich gene family which encodes antimicrobial peptides essential for terminal bacteroid differentiation during nitrogen-fixing symbiosis.

CONCLUSIONS

Our results provide a near-complete R108 genome assembly and the first genomic landscape of SVs obtained by comparing two M. truncatula ecotypes. This may provide valuable genomic resources for the functional and molecular research of legume biology in the future.

Extensive allele mining discovers novel genetic diversity in the loci controlling frost tolerance in barley

Exome sequencing-based allele mining for frost tolerance suggests HvCBF14 rather than CNV at Fr-H2 locus is the main responsible of frost tolerance in barley. Wild relatives, landraces and old cultivars of barley represent a reservoir of untapped and potentially important genes for crop improvement, and the recent sequencing technologies provide the opportunity to mine the existing genetic diversity and to identify new genes/alleles for the traits of interest. In the present study, we use frost tolerance and vernalization requirement as case studies to demonstrate the power of allele mining carried out on exome sequencing data generated from > 400 barley accessions. New deletions in the first intron of VRN-H1 were identified and linked to a reduced vernalization requirement, while the allelic diversity of HvCBF2a, HvCBF4b and HvCBF14 was investigated by combining the analysis of SNPs and read counts. This approach has proven very effective to identify gene paralogs and copy number variants of HvCBF2 and the HvCBF4b-HvCBF2a segment. A multiple linear regression model which considers allelic variation at these genes suggests a major involvement of HvCBF14, rather than copy number variation of HvCBF4b-HvCBF2a, in controlling frost tolerance in barley. Overall, the present study provides powerful resource and tools to discover novel alleles at relevant genes in barley.

Genome-wide exploration of oil biosynthesis genes in cultivated olive tree varieties (Olea europaea): insights into regulation of oil biosynthesis

Genome-wide oil biosynthesis was explored by de novo sequencing two cultivated olive tree (Olea europaea) varieties (cv. Ayvalik and Picual). This is the first report of the former variety sequencing. As outgroups, raw reads of cv. Leccino and scaffold-level assembly of cv. Farga were also retrieved. Each of these four cultivars was chromosome-scale assembled into 23 pseudochromosomes, with 1.31 Gbp (Farga), 0.93 Gbp (Ayvalik), 0.7 Gbp (Picual), and 0.54 Gbp (Leccino) in size. Ab initio gene finding was performed on these assemblies, using wild olive tree (oleaster)-trained parameters. High numbers of gene models were predicted and anchored to the pseudochromosomes: 69,028 (Ayvalik), 55,073 (Picual), 63,785 (Farga), and 40,449 (Leccino). Using previously reported oil biosynthesis genes from wild olive tree genome project, the following homologous sequences were identified: 1,355 (Ayvalik), 1,269 (Farga), 812 (Leccino), and 774 (Picual). Of these, 358 sequences were commonly shared by all cultivars. Besides, some sequences were cultivar unique: Ayvalik (126), Farga (118), Leccino (46), and Picual (52). These putative sequences were assigned to various GO terms, ranging from lipid metabolism to stress tolerance, from signal transactions to development, and to many others, implicating that oil biosynthesis is synergistically regulated with involvement of various other pathways.

Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families

Long-read DNA sequencing technologies require high molecular weight (HMW) DNA of adequate purity and integrity, which can be difficult to isolate from plant material. Plant leaves usually contain high levels of carbohydrates and secondary metabolites that can impact DNA purity, affecting downstream applications. Several protocols and kits are available for HMW DNA extraction, but they usually require a high amount of input material and often lead to substantial DNA fragmentation, making sequencing suboptimal in terms of read length and data yield. We here describe a protocol for plant HMW DNA extraction from low input material (0.1 g) which is easy to follow and quick (2.5 h). This method successfully enabled us to extract HMW from four species from different families (Orchidaceae, Poaceae, Brassicaceae, Asteraceae). In the case of recalcitrant species, we show that an additional purification step is sufficient to deliver a clean DNA sample. We demonstrate the suitability of our protocol for long-read sequencing on the Oxford Nanopore Technologies PromethION^® platform, with and without the use of a short fragment depletion kit.

An SGSGeneloss-Based Method for Constructing a Gene Presence-Absence Table Using Mosdepth

Presence-absence variants (PAV) are genomic regions present in some individuals of a species, but not others. PAVs have been shown to contribute to genomic diversity, especially in bacteria and plants. These structural variations have been linked to traits and can be used to track a species' evolutionary history. PAVs are usually called by aligning short read sequence data from one or more individuals to a reference genome or pangenome assembly, and then comparing coverage. Regions where reads do not align define absence in that individual, and the regions are classified as PAVs. The method below details how to align sequence reads to a reference and how to use the sequencing-coverage calculator Mosdepth to identify PAVs and construct a PAV table for use in downstream comparative genome analysis.

Making a Pangenome Using the Iterative Mapping Approach

Pangenomes have replaced single reference genomes as genetic references, as they contain a better scope of the diversity found in a single species. This protocol outlines the iterative mapping approach in constructing a pangenome, including how to check the raw data, align the data to a reference, how to assemble the data, and how to remove potential contaminants from the final assembly.

Omics-Facilitated Crop Improvement for Climate Resilience and Superior Nutritive Value

Novel crop improvement approaches, including those that facilitate for the exploitation of crop wild relatives and underutilized species harboring the much-needed natural allelic variation are indispensable if we are to develop climate-smart crops with enhanced abiotic and biotic stress tolerance, higher nutritive value, and superior traits of agronomic importance. Top among these approaches are the "omics" technologies, including genomics, transcriptomics, proteomics, metabolomics, phenomics, and their integration, whose deployment has been vital in revealing several key genes, proteins and metabolic pathways underlying numerous traits of agronomic importance, and aiding marker-assisted breeding in major crop species. Here, citing several relevant examples, we appraise our understanding on the recent developments in omics technologies and how they are driving our quest to breed climate resilient crops. Large-scale genome resequencing, pan-genomes and genome-wide association studies are aiding the identification and analysis of species-level genome variations, whilst RNA-sequencing driven transcriptomics has provided unprecedented opportunities for conducting crop abiotic and biotic stress response studies. Meanwhile, single cell transcriptomics is slowly becoming an indispensable tool for decoding cell-specific stress responses, although several technical and experimental design challenges still need to be resolved. Additionally, the refinement of the conventional techniques and advent of modern, high-resolution proteomics technologies necessitated a gradual shift from the general descriptive studies of plant protein abundances to large scale analysis of protein-metabolite interactions. Especially, metabolomics is currently receiving special attention, owing to the role metabolites play as metabolic intermediates and close links to the phenotypic expression. Further, high throughput phenomics applications are driving the targeting of new research domains such as root system architecture analysis, and exploration of plant root-associated microbes for improved crop health and climate resilience. Overall, coupling these multi-omics technologies to modern plant breeding and genetic engineering methods ensures an all-encompassing approach to developing nutritionally-rich and climate-smart crops whose productivity can sustainably and sufficiently meet the current and future food, nutrition and energy demands.

Calling large indels in 1047 Arabidopsis with IndelEnsembler

Large indels greatly impact the observable phenotypes in different organisms including plants and human. Hence, extracting large indels with high precision and sensitivity is important. Here, we developed IndelEnsembler to detect large indels in 1047 Arabidopsis whole-genome sequencing data. IndelEnsembler identified 34 093 deletions, 12 913 tandem duplications and 9773 insertions. Our large indel dataset was more comprehensive and accurate compared with the previous dataset of AthCNV (1). We captured nearly twice of the ground truth deletions and on average 27% more ground truth duplications compared with AthCNV, though our dataset has less number of large indels compared with AthCNV. Our large indels were positively correlated with transposon elements across the Arabidopsis genome. The non-homologous recombination events were the major formation mechanism of deletions in Arabidopsis genome. The Neighbor joining (NJ) tree constructed based on IndelEnsembler's deletions clearly divided the geographic subgroups of 1047 Arabidopsis. More importantly, our large indels represent a previously unassessed source of genetic variation. Approximately 49% of the deletions have low linkage disequilibrium (LD) with surrounding single nucleotide polymorphisms. Some of them could affect trait performance. For instance, using deletion-based genome-wide association study (DEL-GWAS), the accessions containing a 182-bp deletion in AT1G11520 had delayed flowering time and all accessions in north Sweden had the 182-bp deletion. We also found the accessions with 65-bp deletion in the first exon of AT4G00650 (FRI) flowered earlier than those without it. These two deletions cannot be detected in AthCNV and, interestingly, they do not co-occur in any Arabidopsis thaliana accession. By SNP-GWAS, surrounding SNPs of these two deletions do not correlate with flowering time. This example demonstrated that existing large indel datasets miss phenotypic variations and our large indel dataset filled in the gap.

Omics for the Improvement of Abiotic, Biotic, and Agronomic Traits in Major Cereal Crops: Applications, Challenges, and Prospects

Omics technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal crop breeding program, as they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in omics assure time efficiency and cost benefits to improve cereal crops. This review provides a comprehensive overview of the established omics methods in five major cereals, namely rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants; (2) high-density transcriptomics data to study gene expression patterns; (3) global and targeted proteome profiling to study protein structure and interaction; (4) metabolomic profiling to quantify organ-level, small-density metabolites, and their composition; and (5) high-resolution, high-throughput, image-based phenomics approaches are surveyed in this review.

Next-Generation Breeding Strategies for Climate-Ready Crops

Climate change is a threat to global food security due to the reduction of crop productivity around the globe. Food security is a matter of concern for stakeholders and policymakers as the global population is predicted to bypass 10 billion in the coming years. Crop improvement via modern breeding techniques along with efficient agronomic practices innovations in microbiome applications, and exploiting the natural variations in underutilized crops is an excellent way forward to fulfill future food requirements. In this review, we describe the next-generation breeding tools that can be used to increase crop production by developing climate-resilient superior genotypes to cope with the future challenges of global food security. Recent innovations in genomic-assisted breeding (GAB) strategies allow the construction of highly annotated crop pan-genomes to give a snapshot of the full landscape of genetic diversity (GD) and recapture the lost gene repertoire of a species. Pan-genomes provide new platforms to exploit these unique genes or genetic variation for optimizing breeding programs. The advent of next-generation clustered regularly interspaced short palindromic repeat/CRISPR-associated (CRISPR/Cas) systems, such as prime editing, base editing, and de nova domestication, has institutionalized the idea that genome editing is revamped for crop improvement. Also, the availability of versatile Cas orthologs, including Cas9, Cas12, Cas13, and Cas14, improved the editing efficiency. Now, the CRISPR/Cas systems have numerous applications in crop research and successfully edit the major crop to develop resistance against abiotic and biotic stress. By adopting high-throughput phenotyping approaches and big data analytics tools like artificial intelligence (AI) and machine learning (ML), agriculture is heading toward automation or digitalization. The integration of speed breeding with genomic and phenomic tools can allow rapid gene identifications and ultimately accelerate crop improvement programs. In addition, the integration of next-generation multidisciplinary breeding platforms can open exciting avenues to develop climate-ready crops toward global food security.

Chasing breeding footprints through structural variations in Cucumis melo and wild relatives

Cucumis melo (melon or muskmelon) is an important crop in the family of the Cucurbitaceae. Melon is cross pollinated and domesticated at several locations throughout the breeding history, resulting in highly diverse genetic structure in the germplasm. Yet, the relations among the groups and cultivars are still incomplete. We shed light on the melonbreeding history, analyzing structural variations ranging from 50 bp up to 100 kb, identified from whole genome sequences of 100 selected melon accessions and wild relatives. Phylogenetic trees based on SV types completely resolve cultivars and wild accessions into two monophyletic groups and clustering of cultivars largely correlates with their geographic origin. Taking into account morphology, we found six mis-categorized cultivars. Unique inversions are more often shared between cultivars, carrying advantageous genes and do not directly originate from wild species. Approximately 60% of the inversion breaks carry a long poly A/T motif, and following observations in other plant species, suggest that inversions in melon likely resulted from meiotic recombination events. We show that resistance genes in the linkage V region are expanded in the cultivar genomes compared to wild relatives. Furthermore, particular agronomic traits such as fruit ripening, fragrance, and stress response are specifically selected for in the melon subspecies. These results represent distinctive footprints of selective breeding that shaped today's melon. The sequences and genomic relations between land races, wild relatives, and cultivars will serve the community to identify genetic diversity, optimize experimental designs, and enhance crop development.

Building pan-genome infrastructures for crop plants and their use in association genetics

Pan-genomic studies aim at representing the entire sequence diversity within a species to provide useful resources for evolutionary studies, functional genomics and breeding of cultivated plants. Cost reductions in high-throughput sequencing and advances in sequence assembly algorithms have made it possible to create multiple reference genomes along with a catalogue of all forms of genetic variations in plant species with large and complex or polyploid genomes. In this review, we summarize the current approaches to building pan-genomes as an in silico representation of plant sequence diversity and outline relevant methods for their effective utilization in linking structural with phenotypic variation. We propose as future research avenues (i) transcriptomic and epigenomic studies across multiple reference genomes and (ii) the development of user-friendly and feature-rich pan-genome browsers.

Positional-based cloning 'fail-safe' approach is overpowered by wheat chromosome structural variation

Positional-based cloning is a foundational method for understanding the genes and gene networks that control valuable agronomic traits such as grain yield components. In this study, we sought to positionally clone the causal genetic variant of a 1000-grain weight (TGW) quantitative trait loci (QTL) on wheat (Triticum aestivum L.) chromosome arm 5AL. We developed heterogenous inbred families (HIFs) (>5,000 plants) for enhanced genotypic resolution and fine-mapped the QTL to a 10-Mbp region. The transcriptome of developing grains from positive and negative control HIF haplotypes revealed presence-absence chromosome arm 5AS structural variation and unexpectedly no differential expression of genes within the chromosome arm 5AL candidate region. Evaluation of genomic, transcriptomic, and phenotypic data, and predicted function of genes, identified that the 5AL QTL was the result of strong linkage disequilibrium (LD) with chromosome arm 5AS presence or absence (HIF r² = 0.91). Structural variation is common in wheat, and our results highlight that the redundant polyploid genome's masking of such variation is a significant barrier to positional cloning. We propose recommendations for more efficient and robust detection of structural variation, including transitioning from a single nucleotide polymorphism (SNP) to a haplotype-based approach to identify positional cloning targets. We also present nine candidate genes for grain yield components based on chromosome arm 5AS presence or absence, which may unveil hidden variation of homoeolog dosage-dependent genes across the group five chromosome short arms. Taken together, our discovery demonstrates the phenotypic resiliency of polyploid genomic structural variation and highlights a considerable challenge to routine positional cloning in wheat.

Myb10-D confers PHS-3D resistance to pre-harvest sprouting by regulating NCED in ABA biosynthesis pathway of wheat

Pre-harvest sprouting (PHS), the germination of grain before harvest, is a serious problem resulting in wheat yield and quality losses. Here, we mapped the PHS resistance gene PHS-3D from synthetic hexaploid wheat to a 2.4 Mb presence-absence variation (PAV) region and found that its resistance effect was attributed to the pleiotropic Myb10-D by integrated omics and functional analyses. Three haplotypes were detected in this PAV region among 262 worldwide wheat lines and 16 Aegilops tauschii, and the germination percentages of wheat lines containing Myb10-D was approximately 40% lower than that of the other lines. Transcriptome and metabolome profiling indicated that Myb10-D affected the transcription of genes in both the flavonoid and abscisic acid (ABA) biosynthesis pathways, which resulted in increases in flavonoids and ABA in transgenic wheat lines. Myb10-D activates 9-cis-epoxycarotenoid dioxygenase (NCED) by biding the secondary wall MYB-responsive element (SMRE) to promote ABA biosynthesis in early wheat seed development stages. We revealed that the newly discovered function of Myb10-D confers PHS resistance by enhancing ABA biosynthesis to delay germination in wheat. The PAV harboring Myb10-D associated with grain color and PHS will be useful for understanding and selecting white grained PHS resistant wheat cultivars.

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Advances in Cereal Crop Genomics for Resilience under Climate Change

Adapting to climate change, providing sufficient human food and nutritional needs, and securing sufficient energy supplies will call for a radical transformation from the current conventional adaptation approaches to more broad-based and transformative alternatives. This entails diversifying the agricultural system and boosting productivity of major cereal crops through development of climate-resilient cultivars that can sustainably maintain higher yields under climate change conditions, expanding our focus to crop wild relatives, and better exploitation of underutilized crop species. This is facilitated by the recent developments in plant genomics, such as advances in genome sequencing, assembly, and annotation, as well as gene editing technologies, which have increased the availability of high-quality reference genomes for various model and non-model plant species. This has necessitated genomics-assisted breeding of crops, including underutilized species, consequently broadening genetic variation of the available germplasm; improving the discovery of novel alleles controlling important agronomic traits; and enhancing creation of new crop cultivars with improved tolerance to biotic and abiotic stresses and superior nutritive quality. Here, therefore, we summarize these recent developments in plant genomics and their application, with particular reference to cereal crops (including underutilized species). Particularly, we discuss genome sequencing approaches, quantitative trait loci (QTL) mapping and genome-wide association (GWAS) studies, directed mutagenesis, plant non-coding RNAs, precise gene editing technologies such as CRISPR-Cas9, and complementation of crop genotyping by crop phenotyping. We then conclude by providing an outlook that, as we step into the future, high-throughput phenotyping, pan-genomics, transposable elements analysis, and machine learning hold much promise for crop improvements related to climate resilience and nutritional superiority.

SVLR: Genome Structural Variant Detection Using Long-Read Sequencing Data

Genome structural variants (SVs) have great impacts on human phenotype and diversity, and have been linked to numerous diseases. Long-read sequencing technologies arise to make it possible to find SVs of as long as 10,000 nucleotides. Thus, long read-based SV detection has been drawing attention of many recent research projects, and many tools have been developed for long reads to detect SVs recently. In this article, we present a new method, called SVLR, to detect SVs based on long-read sequencing data. Comparing with existing methods, SVLR can detect three new kinds of SVs: block replacements, block interchanges, and translocations. Although these new SVs are structurally more complicated, SVLR achieves accuracies that are comparable with those of the classic SVs. Moreover, for the classic SVs that can be detected by state-of-the-art methods (e.g., SVIM and Sniffles), our experiments demonstrate recall improvements of up to 38% without harming the precisions (i.e., >78%). We also point out three directions to further improve SV detection in the future. Source codes: https://github.com/GWYSDU/SVLR.

Structural variations in papaya genomes

BACKGROUND

Structural variations (SVs) are a type of mutations that have not been widely detected in plant genomes and studies in animals have shown their role in the process of domestication. An in-depth study of SVs will help us to further understand the impact of SVs on the phenotype and environmental adaptability during papaya domestication and provide genomic resources for the development of molecular markers.

RESULTS

We detected a total of 8083 SVs, including 5260 deletions, 552 tandem duplications and 2271 insertions with deletion being the predominant, indicating the universality of deletion in the evolution of papaya genome. The distribution of these SVs is non-random in each chromosome. A total of 1794 genes overlaps with SV, of which 1350 genes are expressed in at least one tissue. The weighted correlation network analysis (WGCNA) of these expressed genes reveals co-expression relationship between SVs-genes and different tissues, and functional enrichment analysis shows their role in biological growth and environmental responses. We also identified some domesticated SVs genes related to environmental adaptability, sexual reproduction, and important agronomic traits during the domestication of papaya. Analysis of artificially selected copy number variant genes (CNV-genes) also revealed genes associated with plant growth and environmental stress.

CONCLUSIONS

SVs played an indispensable role in the process of papaya domestication, especially in the reproduction traits of hermaphrodite plants. The detection of genome-wide SVs and CNV-genes between cultivated gynodioecious populations and wild dioecious populations provides a reference for further understanding of the evolution process from male to hermaphrodite in papaya.

Genomic patterns of structural variation among diverse genotypes of Sorghum bicolor and a potential role for deletions in local adaptation

Genomic structural mutations, especially deletions, are an important source of variation in many species and can play key roles in phenotypic diversification and evolution. Previous work in many plant species has identified multiple instances of structural variations (SVs) occurring in or near genes related to stress response and disease resistance, suggesting a possible role for SVs in local adaptation. Sorghum [Sorghum bicolor (L.) Moench] is one of the most widely grown cereal crops in the world. It has been adapted to an array of different climates as well as bred for multiple purposes, resulting in a striking phenotypic diversity. In this study, we identified genome-wide SVs in the Biomass Association Panel, a collection of 347 diverse sorghum genotypes collected from multiple countries and continents. Using Illumina-based, short-read whole-genome resequencing data from every genotype, we found a total of 24,648 SVs, including 22,359 deletions. The global site frequency spectrum of deletions and other types of SVs fit a model of neutral evolution, suggesting that the majority of these mutations were not under any types of selection. Clustering results based on single nucleotide polymorphisms separated the genotypes into eight clusters which largely corresponded with geographic origins, with many of the large deletions we uncovered being unique to a single cluster. Even though most deletions appeared to be neutral, a handful of cluster-specific deletions were found in genes related to biotic and abiotic stress responses, supporting the possibility that at least some of these deletions contribute to local adaptation in sorghum.

A novel deletion in FLOWERING LOCUS T modulates flowering time in winter oilseed rape

A novel structural variant was discovered in the FLOWERING LOCUS T orthologue BnaFT.A02 by long-read sequencing. Nested association mapping in an elite winter oilseed rape population revealed that this 288 bp deletion associates with early flowering, putatively by modification of binding-sites for important flowering regulation genes. Perfect timing of flowering is crucial for optimal pollination and high seed yield. Extensive previous studies of flowering behavior in Brassica napus (canola, rapeseed) identified mutations in key flowering regulators which differentiate winter, semi-winter and spring ecotypes. However, because these are generally fixed in locally adapted genotypes, they have only limited relevance for fine adjustment of flowering time in elite cultivar gene pools. In crosses between ecotypes, the ecotype-specific major-effect mutations mask minor-effect loci of interest for breeding. Here, we investigated flowering time in a multiparental mapping population derived from seven elite winter oilseed rape cultivars which are fixed for major-effect mutations separating winter-type rapeseed from other ecotypes. Association mapping revealed eight genomic regions on chromosomes A02, C02 and C03 associating with fine modulation of flowering time. Long-read genomic resequencing of the seven parental lines identified seven structural variants coinciding with candidate genes for flowering time within chromosome regions associated with flowering time. Segregation patterns for these variants in the elite multiparental population and a diversity set of winter types using locus-specific assays revealed significant associations with flowering time for three deletions on chromosome A02. One of these was a previously undescribed 288 bp deletion within the second intron of FLOWERING LOCUS T on chromosome A02, emphasizing the advantage of long-read sequencing for detection of structural variants in this size range. Detailed analysis revealed the impact of this specific deletion on flowering-time modulation under extreme environments and varying day lengths in elite, winter-type oilseed rape.

Recent advances in Cannabis sativa genomics research

Cannabis (Cannabis sativa L.) is one of the oldest cultivated plants purported to have unique medicinal properties. However, scientific research of cannabis has been restricted by the Single Convention on Narcotic Drugs of 1961, an international treaty that prohibits the production and supply of narcotic drugs except under license. Legislation governing cannabis cultivation for research, medicinal and even recreational purposes has been relaxed recently in certain jurisdictions. As a result, there is now potential to accelerate cultivar development of this multi-use and potentially medically useful plant species by application of modern genomics technologies. Whilst genomics has been pivotal to our understanding of the basic biology and molecular mechanisms controlling key traits in several crop species, much work is needed for cannabis. In this review we provide a comprehensive summary of key cannabis genomics resources and their applications. We also discuss prospective applications of existing and emerging genomics technologies for accelerating the genetic improvement of cannabis.