Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica

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.

Development of Omni InDel and supporting database for maize

Insertions-deletions (InDels) are the second most abundant molecular marker in the genome and have been widely used in molecular biology research along with simple sequence repeats (SSR) and single-nucleotide polymorphisms (SNP). However, InDel variant mining and marker development usually focuses on a single type of dimorphic InDel, which does not reflect the overall InDel diversity across the genome. Here, we developed Omni InDels for maize, soybean, and rice based on sequencing data and genome assembly that included InDel variants with base lengths from 1 bp to several Mb, and we conducted a detailed classification of Omni InDels. Moreover, we screened a set of InDels that are easily detected and typed (Perfect InDels) from the Omni InDels, verified the site authenticity using 3,587 germplasm resources from 11 groups, and analyzed the germplasm resources. Furthermore, we developed a Multi-InDel set based on the Omni InDels; each Multi-InDel contains multiple InDels, which greatly increases site polymorphism, they can be detected in multiple platforms such as fluorescent capillary electrophoresis and sequencing. Finally, we developed an online database website to make Omni InDels easy to use and share and developed a visual browsing function called "Variant viewer" for all Omni InDel sites to better display the variant distribution.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

FrangiPANe, a tool for creating a panreference using left behind reads

We present here FrangiPANe, a pipeline developed to build panreference using short reads through a map-then-assemble strategy. Applying it to 248 African rice genomes using an improved CG14 reference genome, we identified an average of 8 Mb of new sequences and 5290 new contigs per individual. In total, 1.4 G of new sequences, consisting of 1 306 676 contigs, were assembled. We validated 97.7% of the contigs of the TOG5681 cultivar individual assembly from short reads on a newly long reads genome assembly of the same TOG5681 cultivar. FrangiPANe also allowed the anchoring of 31.5% of the new contigs within the CG14 reference genome, with a 92.5% accuracy at 2 kb span. We annotated in addition 3252 new genes absent from the reference. FrangiPANe was developed as a modular and interactive application to simplify the construction of a panreference using the map-then-assemble approach. It is available as a Docker image containing (i) a Jupyter notebook centralizing codes, documentation and interactive visualization of results, (ii) python scripts and (iii) all the software and libraries requested for each step of the analysis. We foreseen our approach will help leverage large-scale illumina dataset for pangenome studies in GWAS or detection of selection.

ANTHRACNOSE RESISTANCE GENE2 confers fungal resistance in sorghum

Sorghum is an important food and feed crop globally; its production is hampered by anthracnose disease caused by the fungal pathogen Colletotrichum sublineola (Cs). Here, we report identification and characterization of ANTHRACNOSE RESISTANCE GENE 2 (ARG2) encoding a nucleotide-binding leucine-rich repeat (NLR) protein that confers race-specific resistance to Cs strains. ARG2 is one of a cluster of several NLR genes initially identified in the sorghum differential line SC328C that is resistant to some Cs strains. This cluster shows structural and copy number variations in different sorghum genotypes. Different sorghum lines carrying independent ARG2 alleles provided the genetic validation for the identity of the ARG2 gene. ARG2 expression is induced by Cs, and chitin induces ARG2 expression in resistant but not in susceptible lines. ARG2-mediated resistance is accompanied by higher expression of defense and secondary metabolite genes at early stages of infection, and anthocyanin and zeatin metabolisms are upregulated in resistant plants. Interestingly, ARG2 localizes to the plasma membrane when transiently expressed in Nicotiana benthamiana. Importantly, ARG2 plants produced higher shoot dry matter than near-isogenic lines carrying the susceptible allele suggesting an absence of an ARG2 associated growth trade-off. Furthermore, ARG2-mediated resistance is stable at a wide range of temperatures. Our observations open avenues for resistance breeding and for dissecting mechanisms of resistance.

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.

Development of an inclusive 580K SNP array and its application for genomic selection and genome-wide association studies in rice

Rice is a globally cultivated crop and is primarily a staple food source for more than half of the world's population. Various single-nucleotide polymorphism (SNP) arrays have been developed and utilized as standard genotyping methods for rice breeding research. Considering the importance of SNP arrays with more inclusive genetic information for GWAS and genomic selection, we integrated SNPs from eight different data resources: resequencing data from the Korean World Rice Collection (KRICE) of 475 accessions, 3,000 rice genome project (3 K-RGP) data, 700 K high-density rice array, Affymetrix 44 K SNP array, QTARO, Reactome, and plastid and GMO information. The collected SNPs were filtered and selected based on the breeder's interest, covering all key traits or research areas to develop an integrated array system representing inclusive genomic polymorphisms. A total of 581,006 high-quality SNPs were synthesized with an average distance of 200 bp between adjacent SNPs, generating a 580 K Axiom Rice Genotyping Chip (580 K _ KNU chip). Further validation of this array on 4,720 genotypes revealed robust and highly efficient genotyping. This has also been demonstrated in genome-wide association studies (GWAS) and genomic selection (GS) of three traits: clum length, heading date, and panicle length. Several SNPs significantly associated with cut-off, -log10 p-value >7.0, were detected in GWAS, and the GS predictabilities for the three traits were more than 0.5, in both rrBLUP and convolutional neural network (CNN) models. The Axiom 580 K Genotyping array will provide a cost-effective genotyping platform and accelerate rice GWAS and GS studies.

High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence-absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar 'Dempsey' and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.

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.

Pan-genome analysis of three main Chinese chestnut varieties

Chinese chestnut (Castanea mollissima Blume) is one of the earliest domesticated plants of high nutritional and ecological value, yet mechanisms of C. mollissima underlying its growth and development are poorly understood. Although individual chestnut species differ greatly, the molecular basis of the formation of their characteristic traits remains unknown. Though the draft genomes of chestnut have been previously released, the pan-genome of different variety needs to be studied. We report the genome sequence of three cultivated varieties of chestnut herein, namely Hei-Shan-Zhai-7 (H7, drought-resistant variety), Yan-Hong (YH, easy-pruning variety), and Yan-Shan-Zao-Sheng (ZS, early-maturing variety), to expedite convenience and efficiency in its genetics-based breeding. We obtained three chromosome-level chestnut genome assemblies through a combination of Oxford Nanopore technology, Illumina HiSeq X, and Hi-C mapping. The final genome assemblies are 671.99 Mb (YH), 790.99 Mb (ZS), and 678.90 Mb (H7), across 12 chromosomes, with scaffold N50 sizes of 50.50 Mb (YH), 65.05 Mb (ZS), and 52.16 Mb (H7). Through the identification of homologous genes and the cluster analysis of gene families, we found that H7, YH and ZS had 159, 131, and 91 unique gene families, respectively, and there were 13,248 single-copy direct homologous genes in the three chestnut varieties. For the convenience of research, the chestnut genome database was constructed. Based on the results of gene family identification, the presence/absence variations (PAVs) information of the three sample genes was calculated, and a total of 2,364, 2,232, and 1,475 unique genes were identified in H7, YH and ZS, respectively. Our results suggest that the GBSS II-b gene family underwent expansion in chestnut (relative to nearest source species). Overall, we developed high-quality and well-annotated genome sequences of three C. mollissima varieties, which will facilitate clarifying the molecular mechanisms underlying important traits, and shortening the breeding process.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

The Fits and Starts of Indian Rice Domestication: How the Movement of Rice Across Northwest India Impacted Domestication Pathways and Agricultural Stories

Rice is currently the staple food for over 3.5 billion people and is arguably the most important crop exploited by humans. Understanding how we came to the point where a single crop dominates the lives of almost half of the Earth’s population has major significance for our future, even more so given the climatic instability we face today, as rice is a cereal that is dependent on water to an extreme degree. In this study, the nature of early rice agriculture in South Asia is explored, looking at how this critical crop may have begun to be exploited, cultivated, and then brought under agricultural regimes during the long span between c.6500 and 1500 BC. There is now clear evidence for early Holocene cultivation of rice in the Middle Gangetic plains of northern India, but there is still considerable debate about the timing of when this cultivation began and whether it involved domestication of rice. By 3200 BC, however, rice agriculture was present outside the Ganges in the Indus Civilization. The data show accelerated domestication in the Indus environment and agricultural systems that played a part in later hybridization with the arrival of Chinese rice. Understanding how this move from its place of origin to a new environment may have become entangled in the domestication pathways of South Asia rice prior to the arrival of Chinese rice c.1500 BC are important to the overall rice story, as they play into modern concerns relating to biodiversity and different ways of growing and watering rice.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

Expanding Gene-Editing Potential in Crop Improvement with Pangenomes

Pangenomes aim to represent the complete repertoire of the genome diversity present within a species or cohort of species, capturing the genomic structural variance between individuals. This genomic information coupled with phenotypic data can be applied to identify genes and alleles involved with abiotic stress tolerance, disease resistance, and other desirable traits. The characterisation of novel structural variants from pangenomes can support genome editing approaches such as Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR associated protein Cas (CRISPR-Cas), providing functional information on gene sequences and new target sites in variant-specific genes with increased efficiency. This review discusses the application of pangenomes in genome editing and crop improvement, focusing on the potential of pangenomes to accurately identify target genes for CRISPR-Cas editing of plant genomes while avoiding adverse off-target effects. We consider the limitations of applying CRISPR-Cas editing with pangenome references and potential solutions to overcome these limitations.

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.

Current status of structural variation studies in plants

Structural variations (SVs) including gene presence/absence variations and copy number variations are a common feature of genomes in plants and, together with single nucleotide polymorphisms and epigenetic differences, are responsible for the heritable phenotypic diversity observed within and between species. Understanding the contribution of SVs to plant phenotypic variation is important for plant breeders to assist in producing improved varieties. The low resolution of early genetic technologies and inefficient methods have previously limited our understanding of SVs in plants. However, with the rapid expansion in genomic technologies, it is possible to assess SVs with an ever-greater resolution and accuracy. Here, we review the current status of SV studies in plants, examine the roles that SVs play in phenotypic traits, compare current technologies and assess future challenges for SV studies.

The pan-genome of the cultivated soybean (PanSoy) reveals an extraordinarily conserved gene content

Studies on structural variation in plants have revealed the inadequacy of a single reference genome for an entire species and suggest that it is necessary to build a species-representative genome called a pan-genome to better capture the extent of both structural and nucleotide variation. Here, we present a pan-genome of cultivated soybean (Glycine max), termed PanSoy, constructed using the de novo genome assembly of 204 phylogenetically and geographically representative improved accessions selected from the larger GmHapMap collection. PanSoy uncovers 108 Mb (˜11%) of novel nonreference sequences encompassing 3621 protein-coding genes (including 1659 novel genes) absent from the soybean 'Williams 82' reference genome. Nonetheless, the core genome represents an exceptionally large proportion of the genome, with >90.6% of genes being shared by >99% of the accessions. A majority of PAVs encompassing genes could be confirmed with long-read sequencing on a subset of accessions. The PanSoy is a major step towards capturing the extent of genetic variation in cultivated soybean and provides a resource for soybean genomics research and breeding.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

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.

The Genetic Basis of Tomato Aroma

Tomato (Solanum lycopersicum L.) aroma is determined by the interaction of volatile compounds (VOCs) released by the tomato fruits with receptors in the nose, leading to a sensorial impression, such as "sweet", "smoky", or "fruity" aroma. Of the more than 400 VOCs released by tomato fruits, 21 have been reported as main contributors to the perceived tomato aroma. These VOCs can be grouped in five clusters, according to their biosynthetic origins. In the last decades, a vast array of scientific studies has investigated the genetic component of tomato aroma in modern tomato cultivars and their relatives. In this paper we aim to collect, compare, integrate and summarize the available literature on flavour-related QTLs in tomato. Three hundred and 5ifty nine (359) QTLs associated with tomato fruit VOCs were physically mapped on the genome and investigated for the presence of potential candidate genes. This review makes it possible to (i) pinpoint potential donors described in literature for specific traits, (ii) highlight important QTL regions by combining information from different populations, and (iii) pinpoint potential candidate genes. This overview aims to be a valuable resource for researchers aiming to elucidate the genetics underlying tomato flavour and for breeders who aim to improve tomato aroma.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

How the pan-genome is changing crop genomics and improvement

Crop genomics has seen dramatic advances in recent years due to improvements in sequencing technology, assembly methods, and computational resources. These advances have led to the development of new tools to facilitate crop improvement. The study of structural variation within species and the characterization of the pan-genome has revealed extensive genome content variation among individuals within a species that is paradigm shifting to crop genomics and improvement. Here, we review advances in crop genomics and how utilization of these tools is shifting in light of pan-genomes that are becoming available for many crop species.

Whole-genome de novo assemblies reveal extensive structural variations and dynamic organelle-to-nucleus DNA transfers in African and Asian rice

Asian cultivated rice (Oryza sativa) and African cultivated rice (Oryza glaberrima) originated from the wild rice species Oryza rufipogon and Oryza barthii, respectively. The genomes of both cultivated species have undergone profound changes during domestication. Whole-genome de novo assemblies of O. barthii, O. glaberrima, O. rufipogon and Oryza nivara, produced using PacBio single-molecule real-time (SMRT) and next-generation sequencing (NGS) technologies, showed that Gypsy-like retrotransposons are the major contributors to genome size variation in African and Asian rice. Through the detection of genome-wide structural variations (SVs), we observed that besides 28 shared SV hot spots, another 67 hot spots existed in either the Asian or African rice genomes. Based on gene annotation information of the SVs, we established that organelle-to-nucleus DNA transfers resulted in numerous SVs that participated in the nuclear genome divergence of rice species and subspecies. We detected 52 giant nuclear integrants of organelle DNA (NORGs, defined as >10 kb) in six Oryza AA genomes. In addition, we developed an effective method to genotype giant NORGs, based on genome assembly, and first showed the dynamic change in the distribution of giant NORGs in rice natural population. Interestingly, 16 highly differentiated giant NORGs tended to accumulate in natural populations of Asian rice from higher latitude regions, grown at lower temperatures and light intensities. Our study provides new insight into the genome divergence of African and Asian rice, and establishes that organelle-to-nucleus DNA transfers, as potentially powerful contributors to environmental adaptation during rice evolution, play a major role in producing SVs in rice genomes.

Prospects of next generation sequencing in lentil breeding

Lentil is an important food legume crop that has large and complex genome. During past years, considerable attention has been given on the use of next generation sequencing for enriching the genomic resources including identification of SSR and SNP markers, development of unigenes, transcripts, and identification of candidate genes for biotic and abiotic stresses, analysis of genetic diversity and identification of genes/ QTLs for agronomically important traits. However, in other crops including pulses, next generation sequencing has revolutionized the genomic research and helped in genomic assisted breeding rapidly and cost effectively. The present review discuss current status and future prospects of the use NGS based breeding in lentil.

Breeding and biotechnological interventions for trait improvement: status and prospects

Present review describes the molecular tools and strategies deployed in the trait discovery and improvement of major crops. The prospects and challenges associated with these approaches are discussed. Crop improvement relies on modulating the genes and genomic regions underlying key traits, either directly or indirectly. Direct approaches include overexpression, RNA interference, genome editing, etc., while breeding majorly constitutes the indirect approach. With the advent of latest tools and technologies, these strategies could hasten the improvement of crop species. Next-generation sequencing, high-throughput genotyping, precision editing, use of space technology for accelerated growth, etc. had provided a new dimension to crop improvement programmes that work towards delivering better varieties to cope up with the challenges. Also, studies have widened from understanding the response of plants to single stress to combined stress, which provides insights into the molecular mechanisms regulating tolerance to more than one stress at a given point of time. Altogether, next-generation genetics and genomics had made tremendous progress in delivering improved varieties; however, the scope still exists to expand its horizon to other species that remain underutilized. In this context, the present review systematically analyses the different genomics approaches that are deployed for trait discovery and improvement in major species that could serve as a roadmap for executing similar strategies in other crop species. The application, pros, and cons, and scope for improvement of each approach have been discussed with examples, and altogether, the review provides comprehensive coverage on the advances in genomics to meet the ever-growing demands for agricultural produce.

Characterization of the Common Japonica-Originated Genomic Regions in the High-Yielding Varieties Developed from Inter-Subspecific Crosses in Temperate Rice (Oryza sativa L.)

The inter-subspecific crossing between indica and japonica subspecies in rice have been utilized to improve the yield potential of temperate rice. In this study, a comparative study of the genomic regions in the eight high-yielding varieties (HYVs) was conducted with those of the four non-HYVs. The Next-Generation Sequencing (NGS) mapping on the Nipponbare reference genome identified a total of 14 common genomic regions of japonica-originated alleles. Interestingly, the HYVs shared japonica-originated genomic regions on nine chromosomes, although they were developed through different breeding programs. A panel of 94 varieties was classified into four varietal groups with 38 single nucleotide polymorphism (SNP) markers from 38 genes residing in the japonica-originated genomic regions and 16 additional trait-specific SNPs. As expected, the japonica-originated genomic regions were only present in the japonica (JAP) and HYV groups, except for Chr4-1 and Chr4-2. The Wx gene, located within Chr6-1, was present in the HYV and JAP variety groups, while the yield-related genes were conserved as indica alleles in HYVs. The japonica-originated genomic regions and alleles shared by HYVs can be employed in molecular breeding programs to further develop the HYVs in temperate rice.

The broad-spectrum rice blast resistance (R) gene Pita2 encodes a novel R protein unique from Pita

BACKGROUND

Rice blast is generally considered the most devastating rice disease worldwide. The development of resistant varieties has been proven to be the most economical strategy to control the disease. A cluster of resistant (R) genes on rice chromosome 12 including Pita, Pita2 and Ptr has been studies for decades. However, the relationship between these R genes has not been well established.

RESULTS

In this study, we compared the resistance spectra controlled by Pita2 and Pita by testing their monogenic lines (MLs) in four hotspots found in the Philippines and Burundi from 2014 to 2018. The reaction patterns were distinct in two countries and that Pita2-mediated field resistance was relatively prevalent. Pathogenicity tests using 328 single-spore isolates in greenhouse further verified that IRBLta2-Re for Pita2 conferred a relatively broader spectrum resistance than those of Pita. Rough and fine mapping of Pita2 were conducted using F₂ and F₃ populations derived from IRBLta2-Re [CO] and CO 39 consisting of 4344 progeny to delimit Pita2 in a genomic interval flanked by two markers 12 g18530 and 12 g18920 proximal to the centromere of chromosome 12. Alignment of the markers to the genomic sequence of IR64, which harbors Pita2 verified by genetic analysis, approximately delimited the candidate gene(s) within 313-kb genomic fragment. The two Pita2 suppressive mutants that contain mutations within Pita2 were verified and identified. Comparative sequence analysis in these two mutants further identified that each individual allele contains a single nucleotide substitution at a different position resulting in nonsense and missense mutations in the protein product of LOC_Os12g18729. On the contrary, no sequence mutation was detected in other candidate genes, indicating that mutations in LOC_Os12g18729 were responsible for the loss of function of Pita2. Pita2 encodes a novel R protein unique from Pita, which is exactly identical to the previously cloned Ptr. Moreover, based on the resistance gene analysis of rice varieties and mutants containing Pita, it was found that Pita2 rather than Pita was responsible for the specificity to some differential isolates with AvrPita. The diagnosis and survey of Pita2 in IRRI released varieties showed relatively low frequency, implying a high value of its application for breeding resistant varieties against rice blast via marker assisted selection.

CONCLUSION

Our study clarified the relationship between Pita, Pita2 and Ptr. Pita2 is identical to Ptr and distinct from Pita in both sequence and chromosomal location although Pita2 and Pita are genetically linked to each other. The loss of function of Pita2 but not Pita eliminate the specificity to some AvrPita containing isolates, however, the mechanism underlying the recognition between Pita2/Pita and AvrPita remains elusive.

Super-Pangenome by Integrating the Wild Side of a Species for Accelerated Crop Improvement

The pangenome provides genomic variations in the cultivated gene pool for a given species. However, as the crop's gene pool comprises many species, especially wild relatives with diverse genetic stock, here we suggest using accessions from all available species of a given genus for the development of a more comprehensive and complete pangenome, which we refer to as a super-pangenome. The super-pangenome provides a complete genomic variation repertoire of a genus and offers unprecedented opportunities for crop improvement. This opinion article focuses on recent developments in crop pangenomics, the need for a super-pangenome that should include wild species, and its application for crop improvement.

Agricultural systems in Bangladesh: the first archaeobotanical results from Early Historic Wari-Bateshwar and Early Medieval Vikrampura

The present paper reports the first systematic archaeobotanical evidence from Bangladesh together with direct AMS radiocarbon dates on crop remains. Macro-botanical remains were collected by flotation from two sites, Wari-Bateshwar (WB), an Early Historic archaeological site, dating mainly between 400 and 100 BC, with a later seventh century AD temple complex, and Raghurampura Vikrampura (RV), a Buddhist Monastery (vihara) located within the Vikrampura city site complex and dating to the eleventh and sixteenth centuries AD. Despite being a tropical country, with high rainfall and intensive soil processes, our work demonstrates that conventional archaeobotany, the collection of macro-remains through flotation, has much potential towards putting together a history of crops and agricultural systems in Bangladesh. The archaeobotanical assemblage collected from both sites indicates the predominance of rice agriculture, which would have been practiced in summer. Spikelet bases are of domesticated type rice, while grain metrics suggest the majority of rice was probably subspecies japonica. The presence of some wetland weeds suggests at least some of the rice was grown in wet (flooded) systems, but much of it may have been rainfed as inferred from the Southeast Asian weed Acmella paniculata. Other crops include winter cereals, barley and possible oat, and small numbers of summer millets (Pennisetum glaucum, Sorghum bicolor, Setaria italica), a wide diversity of summer and winter pulses (14 spp.), cotton, sesame and mustard seed. Pulse crops included many known from India. Thus, while most crops indicate diffusion of crops from India eastwards, the absence of indica rice could also indicate some diffusion from Southeast Asia. The later site RV also produced evidence of the rice bean (Vigna umbellata), a domesticate of mainland Southeast Asia. These data provide the first empirical evidence for reconstructing past agriculture in Bangladesh and for the role of connections to both India and mainland Southeast Asia in the development of crop diversity in the Ganges delta region.

Identification of Loci Through Genome-Wide Association Studies to Improve Tolerance to Sulfur Deficiency in Rice

Sulfur (S) is an essential nutrient for plant growth and development; however, S supply for crop production is decreasing due to reduced inputs from atmospheric deposition and reduced application of S-containing fertilizers. Sulfur deficiency in soil is therefore becoming a widespread cause of reduced grain yield and quality in rice (Oryza sativa L). We therefore assessed the genotypic variation for tolerance to S deficiency in rice and identified loci associated with improved tolerance. Plants were grown in nutrient solution with either low (0.01 mM) or high (1.0 mM) supply of S. Plants grown under low-S treatment showed a reduction in total biomass, mainly due to a marked reduction in shoot biomass, while root biomass and root-to-shoot ratio increased, relative to plants under high-S treatment. Genome-wide association studies (GWAS) identified loci associated with root length (qSUE2-3, qSUE4, and qSUE9), and root (qSUE1, qSUE2-1, and qSUE3-1 and qSUE3-2) or total dry matter (qSUE2, qSUE3-1, and qSUE11). Candidate genes identified at associated loci coded for enzymes involved in secondary S metabolic pathways (sulfotransferases), wherein the sulfated compounds play several roles in plant responses to abiotic stress; cell wall metabolism including wall loosening and modification (carbohydrate hydrolases: beta-glucosidase and beta-gluconase) important for root growth; and cell detoxification (glutathione S-transferase). This study confirmed the existence of genetic variation conferring tolerance to S deficiency among traditional aus rice varieties. The advantageous haplotypes identified could be exploited through marker assisted breeding to improve tolerance to S-deficiency in modern cultivars in order to achieve sustainable crop production and food security.

Genome Size Variation and Comparative Genomics Reveal Intraspecific Diversity in Brassica rapa

Traditionally, reference genomes in crop species rely on the assembly of one accession, thus occulting most of intraspecific diversity. However, rearrangements, gene duplications, and transposable element content may have a large impact on the genomic structure, which could generate new phenotypic traits. Comparing two Brassica rapa genomes recently sequenced and assembled using long-read technology and optical mapping, we investigated structural variants and repetitive content between the two accessions and genome size variation among a core collection. We explored the structural consequences of the presence of large repeated sequences in B. rapa 'Z1' genome vs. the B. rapa 'Chiifu' genome, using comparative genomics and cytogenetic approaches. First, we showed that large genomic variants on chromosomes A05, A06, A09, and A10 are due to large insertions and inversions when comparing B. rapa 'Z1' and B. rapa 'Chiifu' at the origin of important length differences in some chromosomes. For instance, lengths of 'Z1' and 'Chiifu' A06 chromosomes were estimated in silico to be 55 and 29 Mb, respectively. To validate these observations, we compared using fluorescent in situ hybridization (FISH) the two A06 chromosomes present in an F1 hybrid produced by crossing these two varieties. We confirmed a length difference of 17.6% between the A06 chromosomes of 'Z1' compared to 'Chiifu.' Alternatively, using a copy number variation approach, we were able to quantify the presence of a higher number of rDNA and gypsy elements in 'Z1' genome compared to 'Chiifu' on different chromosomes including A06. Using flow cytometry, the total genome size of 12 Brassica accessions corresponding to a B. rapa available core collection was estimated and revealed a genome size variation of up to 16% between these accessions as well as some shared inversions. This study revealed the contribution of long-read sequencing of new accessions belonging to different cultigroups of B. rapa and highlighted the potential impact of differential insertion of repeat elements and inversions of large genomic regions in genome size intraspecific variability.

1k-RiCA (1K-Rice Custom Amplicon) a novel genotyping amplicon-based SNP assay for genetics and breeding applications in rice

BACKGROUND

While a multitude of genotyping platforms have been developed for rice, the majority of them have not been optimized for breeding where cost, turnaround time, throughput and ease of use, relative to density and informativeness are critical parameters of their utility. With that in mind we report the development of the 1K-Rice Custom Amplicon, or 1k-RiCA, a robust custom sequencing-based amplicon panel of ~ 1000-SNPs that are uniformly distributed across the rice genome, designed to be highly informative within indica rice breeding pools, and tailored for genomic prediction in elite indica rice breeding programs.

RESULTS

Empirical validation tests performed on the 1k-RiCA showed average marker call rates of 95% with marker repeatability and concordance rates of 99%. These technical properties were not affected when two common DNA extraction protocols were used. The average distance between SNPs in the 1k-RiCA was 1.5 cM, similar to the theoretical distance which would be expected between 1,000 uniformly distributed markers across the rice genome. The average minor allele frequencies on a panel of indica lines was 0.36 and polymorphic SNPs estimated on pairwise comparisons between indica by indica accessions and indica by japonica accessions were on average 430 and 450 respectively. The specific design parameters of the 1k-RiCA allow for a detailed view of genetic relationships and unambiguous molecular IDs within indica accessions and good cost vs. marker-density balance for genomic prediction applications in elite indica germplasm. Predictive abilities of Genomic Selection models for flowering time, grain yield, and plant height were on average 0.71, 0.36, and 0.65 respectively based on cross-validation analysis. Furthermore the inclusion of important trait markers associated with 11 different genes and QTL adds value to parental selection in crossing schemes and marker-assisted selection in forward breeding applications.

CONCLUSIONS

This study validated the marker quality and robustness of the 1k-RiCA genotypic platform for genotyping populations derived from indica rice subpopulation for genetic and breeding purposes including MAS and genomic selection. The 1k-RiCA has proven to be an alternative cost-effective genotyping system for breeding applications.

Development of introgression lines of AA genome Oryza species, O. glaberrima, O. rufipogon, and O. nivara, in the genetic background of O. sativa L. cv. Taichung 65

To evaluate and utilize potentially valuable quantitative trait loci or genes of wild relatives in the genetic background of domesticated crop species, chromosome segment substitution lines (CSSLs) are a valuable tool. CSSLs can be constructed through the exchange of chromosome segments of AA genome species of the genus Oryza with cultivated rice, Oryza sativa L. Here we report the development of three sets of CSSLs carrying segments of AA genome species closely related to Oryza sativa-O. glaberrima (IRGC 103777 from Mali), O. rufipogon (W1962 from China), and O. nivara (IRGC 105715 from Cambodia)-in the genetic background of ssp. japonica cultivar Taichung 65 through the use of 101 to 121 simple-sequence-repeat markers in whole-genome genotyping and marker-assisted selection. The materials are available via the National Bioresource Project (Rice) Oryzabase Web page.

Insight into the evolution and functional characteristics of the pan-genome assembly from sesame landraces and modern cultivars

Sesame (Sesamum indicum L.) is an important oil crop renowned for its high oil content and quality. Recently, genome assemblies for five sesame varieties including two landraces (S. indicum cv. Baizhima and Mishuozhima) and three modern cultivars (S. indicum var. Zhongzhi13, Yuzhi11 and Swetha), have become available providing a rich resource for comparative genomic analyses and gene discovery. Here, we employed a reference-assisted assembly approach to improve the draft assemblies of four of the sesame varieties. We then constructed a sesame pan-genome of 554.05 Mb. The pan-genome contained 26 472 orthologous gene clusters; 15 409 (58.21%) of them were core (present across all five sesame genomes), whereas the remaining 41.79% (11 063) clusters and the 15 890 variety-specific genes were dispensable. Comparisons between varieties suggest that modern cultivars from China and India display significant genomic variation. The gene families unique to the sesame modern cultivars contain genes mainly related to yield and quality, while those unique to the landraces contain genes involved in environmental adaptation. Comparative evolutionary analysis indicates that several genes involved in plant-pathogen interaction and lipid metabolism are under positive selection, which may be associated with sesame environmental adaption and selection for high seed oil content. This study of the sesame pan-genome provides insights into the evolution and genomic characteristics of this important oilseed and constitutes a resource for further sesame crop improvement.

Structural variants in 3000 rice genomes

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5′ UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

Rapid DNA-genotyping system targeting ten loci for resistance to blast disease in rice

The fungal pathogen Pyricularia oryzae causes blast, a severe disease of rice (Oryza sativa L.). Improving blast resistance is important in rice breeding programs. Inoculation tests have been used to select for resistance genotypes, with DNA marker-based selection becoming an efficient alternative. No comprehensive DNA marker system for race-specific resistance alleles in the Japanese rice breeding program has been developed because some loci contain multiple resistance alleles. Here, we used the Fluidigm SNP genotyping platform to determine a set of 96 single nucleotide polymorphism (SNP) markers for 10 loci with race-specific resistance. The markers were then used to evaluate the presence or absence of 24 resistance alleles in 369 cultivars; results were 93.5% consistent with reported inoculation test-based genotypes in japonica varieties. The evaluation system was successfully applied to high-yield varieties with indica genetic backgrounds. The system includes polymorphisms that distinguish the resistant alleles at the tightly linked Pita and Pita-2 loci, thereby confirming that all the tested cultivars with Pita-2 allele carry Pita allele. We also developed and validated insertion/deletion (InDel) markers for ten resistance loci. Combining SNP and InDel markers is an accurate and efficient strategy for selection for race-specific resistance to blast in breeding programs.

Prospects of pan-genomics in barley

The concept of a pan-genome refers to intraspecific diversity in genome content and structure, encompassing both genes and intergenic space. Pan-genomic studies employ a combination of de novo sequence assembly and reference-based alignment to discover and genotype structural variants. The large size and complex structure of Triticeae genomes were for a long time an obstacle for genomic research in barley and its relatives. Now that a reference genome is available, computational pipelines for high-quality sequence assembly are in place, and sequence costs continue to drop, investigations into the structural diversity of the barley genome seem within reach. Here, we review the recent progress on pan-genomics in the model grass Brachypodium distachyon, and the cereal crops rice and maize, and devise a multi-tiered strategy for a pan-genome project in barley. Our design involves: (1) the construction of high-quality de novo sequence assemblies for a small core set of representative genotypes, (2) short-read sequencing of a large diversity panel of genebank accessions to medium coverage and (3) the use of complementary methods such as chromosome-conformation capture sequencing and k-mer-based association genetics. The in silico representation of the barley pan-genome may inform about the mechanisms of structural genome evolution in the Triticeae and supplement quantitative genetics models of crop performance for better accuracy and predictive ability.

Tissue-Specific Regulation of Na+ and K+ Transporters Explains Genotypic Differences in Salinity Stress Tolerance in Rice

Rice (Oryza sativa) is a staple food that feeds more than half the world population. As rice is highly sensitive to soil salinity, current trends in soil salinization threaten global food security. To better understand the mechanistic basis of salinity tolerance in rice, three contrasting rice cultivars-Reiziq (tolerant), Doongara (moderately tolerant), and Koshihikari (sensitive)-were examined and the differences in operation of key ion transporters mediating ionic homeostasis in these genotypes were evaluated. Tolerant varieties had reduced Na⁺ translocation from roots to shoots. Electrophysiological and quantitative reverse transcription PCR experiments showed that tolerant genotypes possessed 2-fold higher net Na⁺ efflux capacity in the root elongation zone. Interestingly, this efflux was only partially mediated by the plasma membrane Na⁺/H⁺ antiporter (OsSOS1), suggesting involvement of some other exclusion mechanisms. No significant difference in Na⁺ exclusion from the mature root zones was found between cultivars, and the transcriptional changes in the salt overly sensitive signaling pathway genes in the elongation zone were not correlated with the genetic variability in salinity tolerance amongst genotypes. The most important hallmark of differential salinity tolerance was in the ability of the plant to retain K⁺ in both root zones. This trait was conferred by at least three complementary mechanisms: (1) its superior ability to activate H⁺-ATPase pump operation, both at transcriptional and functional levels; (2) reduced sensitivity of K⁺ efflux channels to reactive oxygen species; and (3) smaller upregulation in OsGORK and higher upregulation of OsAKT1 in tolerant cultivars in response to salt stress. These traits should be targeted in breeding programs aimed to improve salinity tolerance in commercial rice cultivars.

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.

Comparative whole genome re-sequencing analysis in upland New Rice for Africa: insights into the breeding history and respective genome compositions

BACKGROUND

Increasing rice demand is one of the consequences of the steadily improving socio-economic status of the African countries. New Rice for Africa (NERICA), which are interspecific hybrids between Asian and African rice varieties, are one of successful breeding products utilizing biodiversity across the two different rice crop species. Upland NERICA varieties (NU) exhibit agronomic traits of value for the harsh eco-geography, including shorter duration, higher yield and stress tolerance, compared to local African varieties. However, the molecular basis of the traits in NU varieties is largely unknown.

RESULTS

Whole genome re-sequencing was performed of four NU lines (3, 4, 5, and 7) and for the parental Oryza sativa WAB56-104 and Oryza glaberrima CG14. The k-mer analysis predicted large genomes for the four NU lines, most likely inherited from WAB56-104. Approximately 3.1, 0.10, and 0.40 million single nucleotide polymorphisms, multi nucleotide polymorphisms, and short insertions/deletions were mined between the parental lines, respectively. Integrated analysis with another four NU lines (1, 2, 8, and 9) showed that the ratios of the donor CG14 allelic sites in the NU lines ranged from 1.3 to 9.8%. High resolution graphical genotype indicated genome-level similarities and common genetic events during the breeding process: five xyloglucan fucosyltransferase from O. glaberrima were introgressed in common. Segregation of genic segments revealed potential causal genes for some agronomic traits including grain shattering, awnness, susceptibility to bacterial leaf bright, and salt tolerance. Analysis of unmapped sequences against the reference cultivar Nipponbare indicated existence of unique genes for pathogen and abiotic stress resistance in the NU varieties.

CONCLUSIONS

The results provide understanding of NU genomes for rice improvement for Africa reinforcing local capacity for food security and insights into molecular events in breeding of interspecific hybrid crops.

IR64: a high-quality and high-yielding mega variety

High-yielding varieties developed in the 1960s and 1970s at the International Rice Research Institute (IRRI) and elsewhere benefited farmers and the public, ultimately increasing yields and reducing the cost of rice to consumers. Most of these varieties, however, did not have the optimum cooking quality that was possessed by many of the traditional varieties they replaced. In 1985, the IRRI-developed indica variety IR64 was released in the Philippines. In addition to its high yield, early maturity and disease resistance, it had excellent cooking quality, matching that of the best varieties available. These merits resulted in its rapid spread and cultivation on over 10 million ha in the two decades after it was released. It has intermediate amylose content and gelatinization temperature, and good taste. It is resistant to blast and bacterial blight diseases, and to brown planthopper. Because of its success as a variety, it has been used extensively in scientific studies and has been well-characterized genetically. Many valuable genes have been introduced into IR64 through backcross breeding and it has been used in thousands of crosses. Its area of cultivation has declined in the past 10 years, but it has been replaced by a new generation of high-quality varieties that are mostly its progeny or relatives. Continued basic studies on IR64 and related varieties should help in unraveling the complex genetic control of yield and other desirable traits that are prized by rice farmers and consumers.

Exploring traditional aus-type rice for metabolites conferring drought tolerance

BACKGROUND

Traditional varieties and landraces belonging to the aus-type group of rice (Oryza sativa L.) are known to be highly tolerant to environmental stresses, such as drought and heat, and are therefore recognized as a valuable genetic resource for crop improvement. Using two aus-type (Dular, N22) and two drought intolerant irrigated varieties (IR64, IR74) an untargeted metabolomics analysis was conducted to identify drought-responsive metabolites associated with tolerance.

RESULTS

The superior drought tolerance of Dular and N22 compared with the irrigated varieties was confirmed by phenotyping plants grown to maturity after imposing severe drought stress in a dry-down treatment. Dular and N22 did not show a significant reduction in grain yield compared to well-watered control plants, whereas the intolerant varieties showed a significant reduction in both, total spikelet number and grain yield. The metabolomics analysis was conducted with shoot and root samples of plants at the tillering stage at the end of the dry-down treatment. The data revealed an overall higher accumulation of N-rich metabolites (amino acids and nucleotide-related metabolites allantoin and uridine) in shoots of the tolerant varieties. In roots, the aus-type varieties were characterised by a higher reduction of metabolites representative of glycolysis and the TCA cycle, such as malate, glyceric acid and glyceric acid-3-phosphate. On the other hand, the oligosaccharide raffinose showed a higher fold increase in both, shoots and roots of the sensitive genotypes. The data further showed that, for certain drought-responsive metabolites, differences between the contrasting rice varieties were already evident under well-watered control conditions.

CONCLUSIONS

The drought tolerance-related metabolites identified in the aus-type varieties provide a valuable set of protective compounds and an entry point for assessing genetic diversity in the underlying pathways for developing drought tolerant rice and other crops.

Genome Wide Association Mapping of Grain and Straw Biomass Traits in the Rice Bengal and Assam Aus Panel (BAAP) Grown Under Alternate Wetting and Drying and Permanently Flooded Irrigation

Growing demand for staple crops like rice will need to be achieved predominately through agricultural intensification and more efficient use of inputs. To meet this demand it is essential that the genetic diversity within rice is fully utilized. The aus subpopulation is considered an underappreciated resource within that diversity. A new rice panel, the Bengal and Assam Aus Panel (BAAP) of 266 aus accessions was generated with ∼2 million informative SNPs obtained using skim sequencing at ∼4× depth. The BAAP was grown in the field in Bangladesh in the 'boro' season under both continuously flooded and Alternate Wetting and Drying (AWD) irrigation during 2013 and 2014 in Mymensingh and during 2014 in Madhupur. Heading date, grain mass, straw biomass and harvest index were measured. The majority (94%) of BAAP accessions flowered within a relatively small window of 10 days. The AWD irrigation treatment generally caused an increase in grain mass, but no significant genotype by treatment interactions were detected for this trait. Shoot biomass was the only trait that showed evidence of genotype by treatment interaction. The average LD (Linkage Disequilibrium) decay across the genome was 243 Kbp. Genome wide association mapping revealed 115 quantitative trait loci (QTLs). There was little evidence of QTLs specific to the irrigation treatment, and only a few QTLs co-localized with known genes. However, some QTLs were detected across multiple sites and years. These QTLs should be targets for breeding, and include a region around 2.2 Mbp on chromosome 1, a large region in the middle of chromosome 7 and two regions on chromosome 11 (∼10 Mbp and ∼29 Mbp). The BAAP appears to be a valuable addition to the growing collection of GWA mapping populations of rice.

Cruz MT, Kumar A. Association Mapping of Yield and Yield-related Traits Under Reproductive Stage Drought Stress in Rice (Oryza sativa L.)

BACKGROUND

The identification and introgression of major-effect QTLs for grain yield under drought are some of the best and well-proven approaches for improving the drought tolerance of rice varieties. In the present study, we characterized Malaysian rice germplasm for yield and yield-related traits and identified significant trait marker associations by structured association mapping.

RESULTS

The drought screening was successful in screening germplasm with a yield reduction of up to 60% and heritability for grain yield under drought was up to 78%. There was a wider phenotypic and molecular diversity within the panel, indicating the suitability of the population for quantitative trait loci (QTL) mapping. Structure analyses clearly grouped the accessions into three subgroups with admixtures. Linkage disequilibrium (LD) analysis revealed that LD decreased with an increase in distance between marker pairs and the LD decay varied from 5-20 cM. The Mixed Linear model-based structured association mapping identified 80 marker trait associations (MTA) for grain yield (GY), plant height (PH) and days to flowering (DTF). Seven MTA were identified for GY under drought stress, four of these MTA were consistently identified in at least two of the three analyses. Most of these MTA identified were on chromosomes 2, 5, 10, 11 and 12, and their phenotypic variance (PV) varied from 5% to 19%. The in silico analysis of drought QTL regions revealed the association of several drought-responsive genes conferring drought tolerance. The major-effect QTLs are useful in marker-assisted QTL pyramiding to improve drought tolerance.

CONCLUSION

The results have clearly shown that structured association mapping is one of the feasible options to identify major-effect QTLs for drought tolerance-related traits in rice.

Assembly and comparison of two closely related Brassica napus genomes

As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor-bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor-bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.

Association Mapping of Yield and Yield-related Traits Under Reproductive Stage Drought Stress in Rice (Oryza sativa L.)

BACKGROUND

The identification and introgression of major-effect QTLs for grain yield under drought are some of the best and well-proven approaches for improving the drought tolerance of rice varieties. In the present study, we characterized Malaysian rice germplasm for yield and yield-related traits and identified significant trait marker associations by structured association mapping.

RESULTS

The drought screening was successful in screening germplasm with a yield reduction of up to 60% and heritability for grain yield under drought was up to 78%. There was a wider phenotypic and molecular diversity within the panel, indicating the suitability of the population for quantitative trait loci (QTL) mapping. Structure analyses clearly grouped the accessions into three subgroups with admixtures. Linkage disequilibrium (LD) analysis revealed that LD decreased with an increase in distance between marker pairs and the LD decay varied from 5-20 cM. The Mixed Linear model-based structured association mapping identified 80 marker trait associations (MTA) for grain yield (GY), plant height (PH) and days to flowering (DTF). Seven MTA were identified for GY under drought stress, four of these MTA were consistently identified in at least two of the three analyses. Most of these MTA identified were on chromosomes 2, 5, 10, 11 and 12, and their phenotypic variance (PV) varied from 5% to 19%. The in silico analysis of drought QTL regions revealed the association of several drought-responsive genes conferring drought tolerance. The major-effect QTLs are useful in marker-assisted QTL pyramiding to improve drought tolerance.

CONCLUSION

The results have clearly shown that structured association mapping is one of the feasible options to identify major-effect QTLs for drought tolerance-related traits in rice.

Assessing and Exploiting Functional Diversity in Germplasm Pools to Enhance Abiotic Stress Adaptation and Yield in Cereals and Food Legumes

There is a need to accelerate crop improvement by introducing alleles conferring host plant resistance, abiotic stress adaptation, and high yield potential. Elite cultivars, landraces and wild relatives harbor useful genetic variation that needs to be more easily utilized in plant breeding. We review genome-wide approaches for assessing and identifying alleles associated with desirable agronomic traits in diverse germplasm pools of cereals and legumes. Major quantitative trait loci and single nucleotide polymorphisms (SNPs) associated with desirable agronomic traits have been deployed to enhance crop productivity and resilience. These include alleles associated with variation conferring enhanced photoperiod and flowering traits. Genetic variants in the florigen pathway can provide both environmental flexibility and improved yields. SNPs associated with length of growing season and tolerance to abiotic stresses (precipitation, high temperature) are valuable resources for accelerating breeding for drought-prone environments. Both genomic selection and genome editing can also harness allelic diversity and increase productivity by improving multiple traits, including phenology, plant architecture, yield potential and adaptation to abiotic stresses. Discovering rare alleles and useful haplotypes also provides opportunities to enhance abiotic stress adaptation, while epigenetic variation has potential to enhance abiotic stress adaptation and productivity in crops. By reviewing current knowledge on specific traits and their genetic basis, we highlight recent developments in the understanding of crop functional diversity and identify potential candidate genes for future use. The storage and integration of genetic, genomic and phenotypic information will play an important role in ensuring broad and rapid application of novel genetic discoveries by the plant breeding community. Exploiting alleles for yield-related traits would allow improvement of selection efficiency and overall genetic gain of multigenic traits. An integrated approach involving multiple stakeholders specializing in management and utilization of genetic resources, crop breeding, molecular biology and genomics, agronomy, stress tolerance, and reproductive/seed biology will help to address the global challenge of ensuring food security in the face of growing resource demands and climate change induced stresses.

Fine mapping of RYMV3: a new resistance gene to Rice yellow mottle virus from Oryza glaberrima

A new resistance gene against Rice yellow mottle virus was identified and mapped in a 15-kb interval. The best candidate is a CC-NBS-LRR gene. Rice yellow mottle virus (RYMV) disease is a serious constraint to the cultivation of rice in Africa and selection for resistance is considered to be the most effective management strategy. The aim of this study was to characterize the resistance of Tog5307, a highly resistant accession belonging to the African cultivated rice species (Oryza glaberrima), that has none of the previously identified resistance genes to RYMV. The specificity of Tog5307 resistance was analyzed using 18 RYMV isolates. While three of them were able to infect Tog5307 very rapidly, resistance against the others was effective despite infection events attributed to resistance-breakdown or incomplete penetrance of the resistance. Segregation of resistance in an interspecific backcross population derived from a cross between Tog5307 and the susceptible Oryza sativa variety IR64 showed that resistance is dominant and is controlled by a single gene, named RYMV3. RYMV3 was mapped in an approximately 15-kb interval in which two candidate genes, coding for a putative transmembrane protein and a CC-NBS-LRR domain-containing protein, were annotated. Sequencing revealed non-synonymous polymorphisms between Tog5307 and the O. glaberrima susceptible accession CG14 in both candidate genes. An additional resistant O. glaberrima accession, Tog5672, was found to have the Tog5307 genotype for the CC-NBS-LRR gene but not for the putative transmembrane protein gene. Analysis of the cosegregation of Tog5672 resistance with the RYMV3 locus suggests that RYMV3 is also involved in Tog5672 resistance, thereby supporting the CC-NBS-LRR gene as the best candidate for RYMV3.