Population genomics unravels the Holocene history of bread wheat and its relatives

Wheat2035: Integrating pan-omics and advanced biotechnology for future wheat design

Wheat (Triticum aestivum) production is vital for global food security, providing energy and protein to millions of people worldwide. Recent advancements in wheat research have led to significant increases in production, fueled by technological and scientific innovation. Here, we summarize the major advancements in wheat research, particularly the integration of biotechnologies and a deeper understanding of wheat biology. The shift from multi-omics to pan-omics approaches in wheat research has greatly enhanced our understanding of the complex genome, genomic variations, and regulatory networks to decode complex traits. We also outline key scientific questions, potential research directions, and technological strategies for improving wheat over the next decade. Since global wheat production is expected to increase by 60% in 2050, continued innovation and collaboration are crucial. Integrating biotechnologies and a deeper understanding of wheat biology will be essential for addressing future challenges in wheat production, ensuring sustainable practices and improved productivity.

Pan-genome bridges wheat structural variations with habitat and breeding

Wheat is the second largest food crop with a very good breeding system and pedigree record in China. Investigating the genomic footprints of wheat cultivars will unveil potential avenues for future breeding efforts1,2. Here we report chromosome-level genome assemblies of 17 wheat cultivars that chronicle the breeding history of China. Comparative genomic analysis uncovered a wealth of structural rearrangements, identifying 249,976 structural variations with 49.03% (122,567) longer than 5 kb. Cultivars developed in 1980s displayed significant accumulations of structural variations, a pattern linked to the extensive incorporation of European and American varieties into breeding programmes of that era. We further proved that structural variations in the centromere-proximal regions are associated with a reduction of crossover events. We showed that common wheat evolved from spring to winter types via mutations and duplications of the VRN-A1 gene as an adaptation strategy to a changing environment. We confirmed shifts in wheat cultivars linked to dietary preferences, migration and cultural integration in Northwest China. We identified large presence or absence variations of pSc200 tandem repeats on the 1RS terminal, suggesting its own rapid evolution in the wheat genome. The high-quality genome assemblies of 17 representatives developed and their good complementarity to the 10+ pan-genomes offer a robust platform for future genomics-assisted breeding in wheat.

Molecular characterization of a novel photoperiod-insensitive allele Ppd-B1a.3 and its effect on heading date in Chinese wheat (Triticum aestivum) cultivar Qingchun 37

Breeders adjust wheat heading dates to improve regional adaptability and reduce or mitigate yield losses caused by meteorological disasters, pests and diseases. The Ppd-1 genes play a crucial role in determining wheat sensitivity to changes in day-length and serve as key regulators of heading dates once the vernalization requirement is satisfied. In this study, we identified a new allelic variant of the promoter region, Ppd-B1a.3, in the Chinese wheat cultivar Qingchun 37. Compared to the Ppd-B1b.1 (carried by Chihokukomugi), the main mutation sites in Ppd-B1a.3 include a substitution of C with G at the -505-bp, a T base insertion at the -625-bp, a mutation of TCG to GGT at the -632 to -634-bp, and a 163-bp insertion at the -691 bp. Analysis of F2 populations indicated that Ppd-B1a.3 promotes heading and flowering (approximately 6 days earlier in population 1 and 17 days in population 2) under short-day conditions in a greenhouse. However, the evaluation of Ppd-B1a.3's effect under field conditions may be influenced by the copy number of the Ppd-B1 locus inherited from the other parent in the F2 populations. Ppd-B1a.3 disrupts circadian rhythm expression and exhibits a stronger effect on heading and flowering than the three-copy Ppd-B1 allele carried by Jing 411. Origin analysis suggests that Ppd-B1a.3 may have derived from non-native germplasm. These results deepen our understanding of wheat photoperiod genes and provide useful genetic resources for fine-tuning wheat heading dates during breeding.

A comprehensive map of DNA-segment copy number variation in 491 genomes of common wheat uncovers genes associated with multiple agronomic traits

DNA-segment copy number variations (DSCNVs), such as deletions and duplications, are important sources of genomic structural variation. However, the types and sizes of DSCNVs, as well as their genome-wide distribution and potential functions, are poorly understood in wheat. Here, we identified 198 985 DSCNVs by investigating 491 genomes of common wheat and found that they account for 20% of the entire genome. Interestingly, approximately 38% of genes are linked to DSCNVs. The number of DSCNVs within each accession ranges from 47 366 to 96 342, and their total sizes vary from 421.3 to 1267.9 Mb. We found that 957 and 1304 DSCNVs have been favored by breeders in China and the United States, respectively. By conducting DSCNV-based genome-wide association studies for the principal components of plant developmental and yield-component traits, we identified 34 loci as directly or indirectly involved in controlling the formation of multiple traits. Notably, a newly discovered DSCNV covering TaFT-D1 is significantly associated with flowering time and other agronomic traits. Overall, our findings highlight the potential of DSCNVs to drive fundamental discoveries in plant science. The comprehensive DSCNV map and the DSCNV-associated genes will also facilitate future research efforts to improve wheat yield, quality, and adaptation.

Genetic diversity assessment of clonal plant Rosa persica in China

Rosa persica is considered a clonal plant because it is mainly propagated by clonal growth. Due to environmental degradation and habitat devastation, R. persica has been listed as a national second-class protected plant in China. However, the absence of research on wild populations of R. persica has impeded progress in formulating efficient conservation strategies. In this study, we investigated the clonal dispersal distance of R. persica to accurately determine the genetic diversity and population structure of the wild population in Xinjiang. We suggested that 20 m was the threshold distance with which to distinguish between different genets of plants. Based on this, we collated sequencing data from a total of 70 different genets of plants from 117 test samples. Eight populations of R. persica were primarily categorized into three subgroups: BL (Bole), TC (Tacheng) and CG (Changji). Of these, the CG subgroup exhibited the most genetic diversity. This research is the first to illustrates the clonal dispersal distance of R. persica, thus providing valuable reference guidelines for understanding the reproductive characteristics of clonal plants. In addition, the genetic diversity of R. persica provides a theoretical foundation for the formulation of conservation policies.

Genomic insights into the modifications of spike morphology traits during wheat breeding

Over the past century, environmental changes have significantly impacted wheat spike morphology, crucial for adaptation and grain yield. However, the changes in wheat spike modifications during this period remain largely unknown. This study examines 16 spike morphology traits in 830 accessions released from 1900 to 2020. It finds that spike weight, grain number per spike (GN), and thousand kernel weight have significantly increased, while spike length has no significant change. The increase in fertile spikelets is due to fewer degenerated spikelets, resulting in a higher GN. Genome-wide association studies identified 49,994 significant SNPs, grouped into 293 genomic regions. The accumulation of favorable alleles in these genomic regions indicates the genetic basis for modification in spike morphology traits. Genetic network analysis of these genomic regions reveals the genetic basis for phenotypic correlations among spike morphology traits. The haplotypes of the identified genomic regions display obvious geographical differentiation in global accessions and environmental adaptation over the past 120 years. In summary, we reveal the genetic basis of adaptive evolution and the interactions of spike morphology, offering valuable resources for the genetic improvement of spike morphology to enhance environmental adaptation.

On the evolution and genetic diversity of the bread wheat D genome

Bread wheat (Triticum aestivum) became a globally dominant crop after incorporating the D genome from the donor species Aegilops tauschii, but the evolutionary history that shaped the D genome during this process remains to be clarified. Here, we propose a renewed evolutionary model linking Ae. tauschii and the hexaploid wheat D genome by constructing an ancestral haplotype map covering 762 Ae. tauschii and hexaploid wheat accessions. We dissected the evolutionary trajectories of Ae. tauschii lineages and reported a few independent intermediate accessions, demonstrating that low-frequency inter-sublineage gene flow had enriched the diversity of Ae. tauschii. We discovered that the D genome of hexaploid wheat was inherited from a unified ancestral template, but with a mosaic composition that was highly mixed and derived mainly from three Ae. tauschii L2 sublineages located in the Caspian coastal region. This result suggests that early agricultural activities facilitated innovations in D-genome composition and finalized the success of hexaploidization. We found that the majority (51.4%) of genetic diversity was attributed to novel mutations absent in Ae. tauschii, and we identified large Ae. tauschii introgressions from various lineages, which expanded the diversity of the wheat D genome and introduced beneficial alleles. This work sheds light on the process of wheat hexaploidization and highlights the evolutionary significance of the multi-layered genetic diversity of the bread wheat D genome.

Convergence and divergence of diploid and tetraploid cotton genomes

Polyploidy is an important driving force in speciation and evolution; however, the genomic basis for parallel selection of a particular trait between polyploids and ancestral diploids remains unexplored. Here we construct graph-based pan-genomes for diploid (A2) and allotetraploid (AD1) cotton species, enabled by an assembly of 50 genomes of genetically diverse accessions. We delineate a mosaic genome map of tetraploid cultivars that illustrates genomic contributions from semi-wild forms into modern cultivars. Pan-genome comparisons identify syntenic and hyper-divergent regions of continued variation between diploid and tetraploid cottons, and suggest an ongoing process of sequence evolution potentially linked to the contrasting genome size change in two subgenomes. We highlight 43% of genetic regulatory relationships for gene expression in diploid encompassing sequence divergence after polyploidy, and specifically characterize six underexplored convergent genetic loci contributing to parallel selection of fiber quality. This study offers a framework for pan-genomic dissection of genetic regulatory components underlying parallel selection of desirable traits in organisms.

Structural variation discovery in wheat using PacBio high-fidelity sequencing

Structural variations (SVs) pervade plant genomes and contribute substantially to the phenotypic diversity. However, most SVs were ineffectively assayed due to their complex nature and the limitations of early genomic technologies. By applying the PacBio high-fidelity (HiFi) sequencing for wheat genomes, we performed a comprehensive evaluation of mainstream long-read aligners and SV callers in SV detection. The results indicated that the accuracy of deletion discovery is markedly influenced by callers, accounting for 87.73% of the variance, whereas both aligners (38.25%) and callers (49.32%) contributed substantially to the accuracy variance for insertions. Among the aligners, Winnowmap2 and NGMLR excelled in detecting deletions and insertions, respectively. For SV callers, SVIM achieved the best performance. We demonstrated that combining the aligners and callers mentioned above is optimal for SV detection. Furthermore, we evaluated the effect of sequencing depth on the accuracy of SV detection, revealing that low-coverage HiFi sequencing is sufficiently robust for high-quality SV discovery. This study thoroughly evaluated SV discovery approaches and established optimal workflows for investigating structural variations using low-coverage HiFi sequencing in the wheat genome, which will advance SV discovery and decipher the biological functions of SVs in wheat and many other plants.

Genome-wide identification of the E-class gene family in wheat: evolution, expression, and interaction

Introduction

Wheat (Triticum aestivum L.) is among themost important crop worldwide. Given a growing population and changing climate, enhancing wheat yield is of great importance. Yield is closely associated with flower and spike development, and E-class genes play important roles in the flower and kernel development of plants. Currently, the absence of systematic analysis on the E gene family hinders our comprehension of their roles in plant growth and development.

Methods

Identify E-class genes based on homologous sequence searches. Analyze the identified E-class genes through a series of gene family analyses. Determine the expression levels of wheat E-class genes by searching public databases. Validate the functions of these genes by transforming them into Arabidopsis. Finally, determine the interactions between the genes through yeast two-hybrid experiments.

Results

Fifteen E-class genes (TaEs) were identified in common wheat. Nine E-class genes were detected in five ancestral/closely related species, including one in Aegilops tauschii (AtE), one in T. Urartu (TuEs), two in T. turgidum (TtEs), two in T. dicoccoides (TdEs), and three in T. spelta (TsEs). The 24 E-class genes were classified into three subgroups using a phylogenetic approach. All genes were highly expressed in spikes, and most were only highly expressed at the floret meristem stage. The effects of TaSEP5-A on flowering and growth cycles were confirmed in homologous mutants and transgenic Arabidopsis thaliana. The E-class genes were able to regulate the growth cycle of Arabidopsis. Finally, we confirmed the interactions between TaSEP5-A and other wheat E-class genes based on yeast two-hybrid assays.

Discussion

Our findings provide information regarding the E-class genes in wheat and will potentially promote the application of these genes in wheat improvement.

Origin and evolution of the bread wheat D genome

Bread wheat (Triticum aestivum) is a globally dominant crop and major source of calories and proteins for the human diet. Compared with its wild ancestors, modern bread wheat shows lower genetic diversity, caused by polyploidisation, domestication and breeding bottlenecks1,2. Wild wheat relatives represent genetic reservoirs, and harbour diversity and beneficial alleles that have not been incorporated into bread wheat. Here we establish and analyse extensive genome resources for Tausch's goatgrass (Aegilops tauschii), the donor of the bread wheat D genome. Our analysis of 46 Ae. tauschii genomes enabled us to clone a disease resistance gene and perform haplotype analysis across a complex disease resistance locus, allowing us to discern alleles from paralogous gene copies. We also reveal the complex genetic composition and history of the bread wheat D genome, which involves contributions from genetically and geographically discrete Ae. tauschii subpopulations. Together, our results reveal the complex history of the bread wheat D genome and demonstrate the potential of wild relatives in crop improvement.

Flowering time genes branching out

Plants are sessile by nature, and as such they have evolved to sense changes in seasonality and their surrounding environment, and adapt to these changes. One prime example of this is the regulation of flowering time in angiosperms, which is precisely timed by the coordinated action of two proteins: FLOWERING LOCUS T (FT) and TERMINAL FLOWER 1 (TFL1). Both of these regulators are members of the PHOSPHATIDYLETHANOLAMINE BINDING PROTEIN (PEBP) family of proteins. These regulatory proteins do not interact with DNA themselves, but instead interact with transcriptional regulators, such as FLOWERING LOCUS D (FD). FT and TFL1 were initially identified as key regulators of flowering time, acting through binding with FD; however, PEBP family members are also involved in shaping plant architecture and development. In addition, PEBPs can interact with TCP transcriptional regulators, such as TEOSINTE BRANCHED 1 (TB1), a well-known regulator of plant architecture, and key domestication-related genes in many crops. Here, we review the role of PEBPs in flowering time, plant architecture, and development. As these are also key yield-related traits, we highlight examples from the model plant Arabidopsis as well as important food and feed crops such as, rice, barley, wheat, tomato, and potato.

Adaptation and the Geographic Spread of Crop Species

Crops are plant species that were domesticated starting about 11,000 years ago from several centers of origin, most prominently the Fertile Crescent, East Asia, and Mesoamerica. From their domestication centers, these crops spread across the globe and had to adapt to differing environments as a result of this dispersal. We discuss broad patterns of crop spread, including the early diffusion of crops associated with the rise and spread of agriculture, the later movement via ancient trading networks, and the exchange between the Old and New Worlds over the last ∼550 years after the European colonization of the Americas. We also examine the various genetic mechanisms associated with the evolutionary adaptation of crops to their new environments after dispersal, most prominently seasonal adaptation associated with movement across latitudes, as well as altitudinal, temperature, and other environmental factors.

Evolutionary genomics of two diploid goat grass species belonging to the section Sitopsis of Aegilops, Aegilops longissima, and Aegilops sharonensis

Aegilops longissima and Ae. sharonensis, being classified into the Sitopsis section of genus Aegilops, are distinct species both taxonomically and ecologically. Nevertheless, earlier observations indicate that the two species are not reproductively isolated to full extent and can inter-bred upon secondary contact. However, the genomic underpinnings of the morpho-ecological differentiation between the two foci species remained unexplored. Here, we resequenced 31 representative accessions of the two species and conducted in-depth comparative genomic analyses. We demonstrate recurrent and ongoing natural hybridizations between Ae. longissima and Ae. sharonensis, and depict features of genome composition of the resultant hybrids at both individual and population levels. We also delineate genomic regions and candidate genes potentially underpinning the differential morphological and edaphic adaptations of the two species. Intriguingly, a binary morphology was observed in the hybrids, suggesting existence of highly diverged genomic regions that remain uneroded by the admixtures. Together, our results provide new insights into the molding effects of interspecific hybridization on genome composition and mechanisms preventing merge of the two species.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation can empower the breeding of climate-resilient crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages, shows remarkable adaptability across diverse global environments, making it an excellent model for investigating species responses to climate change. In this study, we performed population genomic analyses using genome resequencing data from 702 accessions of 24 Medicago species to unravel alfalfa's climatic adaptation and genetic susceptibility to future climate change. We found that interspecific genetic exchange has contributed to the gene pool of alfalfa, particularly enriching defense and stress-response genes. Intersubspecific introgression between M. sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa's climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgressions from subsp. falcata. By integrating climate-associated variants and climate data, we identified populations that are vulnerable to future climate change, particularly in higher latitudes of the Northern Hemisphere. These findings serve as a clarion call for targeted conservation initiatives and breeding efforts. We also identified pre-adaptive populations that demonstrate heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. Collectively, this study enhances our understanding about the local adaptation mechanisms of alfalfa and facilitates the breeding of climate-resilient alfalfa cultivars, contributing to effective agricultural strategies for facing future climate change.

Variation of TaMyb10 and their function on grain color and pre-harvest sprouting resistance of wheat

Pre-harvest sprouting (PHS) is a significant threat to global food security due to its association with losses in both yield and quality. Among the genes involved in PHS resistance in wheat, PHS-3D (TaMyb10-D) plays a crucial role. Here, we characterized the sequence variations of TaMyb10 genes in 416 bread wheat and 302 Aegilops tauschii accessions. Within TaMyb10-A sequences, we identified a deletion ranging from 214 to 305 bp in the signal and amino acid coding region, present in 61.3% of the accessions. Similarly, 79.3% of the TaMyb10-B sequences within the third exon region exhibited a 19 bp deletion. Additionally, 40.8% of the accessions lacked the 2.4 Mb fragment (in/del mutations) on Chr3D, where TaMyb10-D/PHS-3D was located. Interestingly, the geographical distribution of accessions showed little correlation with the divergence of TaMyb10. TaMyb10-A-IIIDele, TaMyb10-B-IVDele, and TaMyb10-D-VDele genotypes were prevalent in wheat populations across continents. Despite their structural variations, the five distinct protein types exhibited comparable ability to bind the promoters of downstream genes in the flavonoid and ABA pathways, such as CHS, DFR, and NCED. Furthermore, the combination of TaMyb10 homologs was significantly associated with grain color and germination percentages. Accessions exclusively harboring TaMyb10-D displayed red seed color and reduced germination percentages, indicating the predominant role of TaMyb10-D compared to TaMyb10-A and TaMyb10-B. This comprehensive investigation enhances our understanding of the structural variations and functional divergence of TaMyb10, providing valuable insights and resources for improving PHS resistance in wheat.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Needs for a conceptual bridge between biological domestication and early food globalization

The past 15 y has seen much development in documentation of domestication of plants and animals as gradual traditions spanning millennia. There has also been considerable momentum in understanding the dispersals of major domesticated taxa across continents spanning thousands of miles. The two processes are often considered within different theoretical strains. What is missing from our repertoire of explanations is a conceptual bridge between the protracted process over millennia and the multiregional, globally dispersed nature of domestication. The evidence reviewed in this paper bears upon how we conceptualize domestication as an episode or a process. By bringing together the topics of crop domestication and crop movement, those complex, protracted, and continuous outcomes come more clearly into view.

Diversification of quantitative morphological traits in wheat

BACKGROUND AND AIMS

The development and morphology of crop plants have been profoundly altered by evolution under cultivation, initially through unconscious selection, without deliberate foresight, and later by directed breeding. Wild wheats remain an important potential source of variation for modern breeders; however, the sequence and timing of morphological changes during domestication are not fully resolved.

METHODS

We grew and measured 142 wheat accessions representing different stages in wheat evolution, including three independent domestication events, and compared their morphological traits to define the morphospace of each group.

KEY RESULTS

The results show that wild and domesticated wheats have overlapping morphospaces, but each also occupies a distinct area of morphospace from one another. Polyploid formation in wheat increased leaf biomass and seed weight but had its largest effects on tiller loss. Domestication continued to increase the sizes of wheat leaves and seeds and made wheat grow taller, with more erect architecture. Associated changes to the biomass of domesticated wheats generated more grains and achieved higher yields. Landrace improvement subsequently decreased the numbers of tillers and spikes, to focus resource allocation to the main stem, accompanied by a thicker main stem and larger flag leaves. During the Green Revolution, wheat height was reduced to increase the harvest index and therefore yield. Modern wheats also have more erect leaves and larger flower biomass proportions than landraces.

CONCLUSIONS

Quantitative trait history in wheat differs by trait. Some trait values show progressive changes in the same direction (e.g. leaf size, grain weight), whereas others change in a punctuated way at particular stages (e.g. canopy architecture), and other trait values switch directions during wheat evolution (e.g. plant height, flower biomass proportion). Agronomically valued domestication traits arose during different stages of wheat history, such that modern wheats are the product of >10 000 years of morphological evolution.

A comprehensive benchmark of graph-based genetic variant genotyping algorithms on plant genomes for creating an accurate ensemble pipeline

BACKGROUND

Although sequencing technologies have boosted the measurement of the genomic diversity of plant crops, it remains challenging to accurately genotype millions of genetic variants, especially structural variations, with only short reads. In recent years, many graph-based variation genotyping methods have been developed to address this issue and tested for human genomes. However, their performance in plant genomes remains largely elusive. Furthermore, pipelines integrating the advantages of current genotyping methods might be required, considering the different complexity of plant genomes.

RESULTS

Here we comprehensively evaluate eight such genotypers in different scenarios in terms of variant type and size, sequencing parameters, genomic context, and complexity, as well as graph size, using both simulated and real data sets from representative plant genomes. Our evaluation reveals that there are still great challenges to applying existing methods to plants, such as excessive repeats and variants or high resource consumption. Therefore, we propose a pipeline called Ensemble Variant Genotyper (EVG) that can achieve better genotyping performance in almost all experimental scenarios and comparably higher genotyping recall and precision even using 5× reads. Furthermore, we demonstrate that EVG is more robust with an increasing number of graphed genomes, especially for insertions and deletions.

CONCLUSIONS

Our study will provide new insights into the development and application of graph-based genotyping algorithms. We conclude that EVG provides an accurate, unbiased, and cost-effective way for genotyping both small and large variations and will be potentially used in population-scale genotyping for large, repetitive, and heterozygous plant genomes.

Genomic and pangenomic analyses provide insights into the population history and genomic diversification of bottle gourd

Bottle gourd (Lagenaria siceraria (Mol.) Strandl.) is an economically important vegetable crop and one of the earliest domesticated crops. However, the population history and genomic diversification of bottle gourd have not been extensively studied. We generated a comprehensive bottle gourd genome variation map from genome sequences of 197 world-wide representative accessions, which enables a genome-wide association study for identifying genomic loci associated with resistance to zucchini yellow mosaic virus, and constructed a bottle gourd pangenome that harbors 1534 protein-coding genes absent in the reference genome. Demographic analyses uncover that domesticated bottle gourd originated in Southern Africa c. 12 000 yr ago, and subsequently radiated to the New World via the Atlantic drift and to Eurasia through the efforts of early farmers in the initial Holocene. The identified highly differentiated genomic regions among different bottle gourd populations harbor many genes contributing to their local adaptations such as those related to disease resistance and stress tolerance. Presence/absence variation analysis of genes in the pangenome reveals numerous genes including those involved in abiotic/biotic stress responses that have been under selection during the world-wide expansion of bottle gourds. The bottle gourd variation map and pangenome provide valuable resources for future functional studies and genomics-assisted breeding.

Population genomic analysis unravels the evolutionary roadmap of pericarp color in rice

Pigmented rice stands out for its nutritional value and is gaining more and more attention. Wild rice, domesticated red rice, and weedy rice all have a red pericarp and a comprehensive genetic background in terms of the red-pericarp phenotype. We performed population genetic analyses using 5104 worldwide rice accessions, including 2794 accessions with red or black pericarps, 85 of which were newly sequenced in this study. The results suggested an evolutionary trajectory of red landraces originating from wild rice, and the split times of cultivated red and white rice populations were estimated to be within the past 3500 years. Cultivated red rice was found to feralize to weedy rice, and weedy rice could be further re-domesticated to cultivated red rice. A genome-wide association study based on the 2794 accessions with pigmented pericarps revealed several new candidate genes associated with the red-pericarp trait for further functional characterization. Our results provide genomic evidence for the origin of pigmented rice and a valuable genomic resource for genetic investigation and breeding of pigmented rice.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Whole-genome sequencing of diverse wheat accessions uncovers genetic changes during modern breeding in China and the United States

Breeding has dramatically changed the plant architecture of wheat (Triticum aestivum), resulting in the development of high-yielding varieties adapted to modern farming systems. However, how wheat breeding shaped the genomic architecture of this crop remains poorly understood. Here, we performed a comprehensive comparative analysis of a whole-genome resequencing panel of 355 common wheat accessions (representing diverse landraces and modern cultivars from China and the United States) at the phenotypic and genomic levels. The genetic diversity of modern wheat cultivars was clearly reduced compared to landraces. Consistent with these genetic changes, most phenotypes of cultivars from China and the United States were significantly altered. Of the 21 agronomic traits investigated, 8 showed convergent changes between the 2 countries. Moreover, of the 207 loci associated with these 21 traits, more than half overlapped with genomic regions that showed evidence of selection. The distribution of selected loci between the Chinese and American cultivars suggests that breeding for increased productivity in these 2 regions was accomplished by pyramiding both shared and region-specific variants. This work provides a framework to understand the genetic architecture of the adaptation of wheat to diverse agricultural production environments, as well as guidelines for optimizing breeding strategies to design better wheat varieties.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Wheat adaptation to environmental stresses under climate change: Molecular basis and genetic improvement

Wheat (Triticum aestivum) is a staple food for about 40% of the world's population. As the global population has grown and living standards improved, high yield and improved nutritional quality have become the main targets for wheat breeding. However, wheat production has been compromised by global warming through the more frequent occurrence of extreme temperature events, which have increased water scarcity, aggravated soil salinization, caused plants to be more vulnerable to diseases, and directly reduced plant fertility and suppressed yield. One promising option to address these challenges is the genetic improvement of wheat for enhanced resistance to environmental stress. Several decades of progress in genomics and genetic engineering has tremendously advanced our understanding of the molecular and genetic mechanisms underlying abiotic and biotic stress responses in wheat. These advances have heralded what might be considered a "golden age" of functional genomics for the genetic improvement of wheat. Here, we summarize the current knowledge on the molecular and genetic basis of wheat resistance to abiotic and biotic stresses, including the QTLs/genes involved, their functional and regulatory mechanisms, and strategies for genetic modification of wheat for improved stress resistance. In addition, we also provide perspectives on some key challenges that need to be addressed.

A high-resolution genotype-phenotype map identifies the TaSPL17 controlling grain number and size in wheat

BACKGROUND

Large-scale genotype-phenotype association studies of crop germplasm are important for identifying alleles associated with favorable traits. The limited number of single-nucleotide polymorphisms (SNPs) in most wheat genome-wide association studies (GWASs) restricts their power to detect marker-trait associations. Additionally, only a few genes regulating grain number per spikelet have been reported due to sensitivity of this trait to variable environments.

RESULTS

We perform a large-scale GWAS using approximately 40 million filtered SNPs for 27 spike morphology traits. We detect 132,086 significant marker-trait associations and the associated SNP markers are located within 590 associated peaks. We detect additional and stronger peaks by dividing spike morphology into sub-traits relative to GWAS results of spike morphology traits. We propose that the genetic dissection of spike morphology is a powerful strategy to detect signals for grain yield traits in wheat. The GWAS results reveal that TaSPL17 positively controls grain size and number by regulating spikelet and floret meristem development, which in turn leads to enhanced grain yield per plant. The haplotypes at TaSPL17 indicate geographical differentiation, domestication effects, and breeding selection.

CONCLUSION

Our study provides valuable resources for genetic improvement of spike morphology and a fast-forward genetic solution for candidate gene detection and cloning in wheat.

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

Einkorn genomics sheds light on history of the oldest domesticated wheat

Einkorn (Triticum monococcum) was the first domesticated wheat species, and was central to the birth of agriculture and the Neolithic Revolution in the Fertile Crescent around 10,000 years ago1,2. Here we generate and analyse 5.2-Gb genome assemblies for wild and domesticated einkorn, including completely assembled centromeres. Einkorn centromeres are highly dynamic, showing evidence of ancient and recent centromere shifts caused by structural rearrangements. Whole-genome sequencing analysis of a diversity panel uncovered the population structure and evolutionary history of einkorn, revealing complex patterns of hybridizations and introgressions after the dispersal of domesticated einkorn from the Fertile Crescent. We also show that around 1% of the modern bread wheat (Triticum aestivum) A subgenome originates from einkorn. These resources and findings highlight the history of einkorn evolution and provide a basis to accelerate the genomics-assisted improvement of einkorn and bread wheat.

Special Issue "Strategies to Develop High-Quality Gluten-Free Products Welcomed by Consumers"

Extensive and long-term efforts on wheat breeding [...].

Triticeae crop genome biology: an endless frontier

Triticeae, the wheatgrass tribe, includes several major cereal crops and their wild relatives. Major crops within the Triticeae are wheat, barley, rye, and oat, which are important for human consumption, animal feed, and rangeland protection. Species within this tribe are known for their large genomes and complex genetic histories. Powered by recent advances in sequencing technology, researchers worldwide have made progress in elucidating the genomes of Triticeae crops. In addition to assemblies of high-quality reference genomes, pan-genome studies have just started to capture the genomic diversities of these species, shedding light on our understanding of the genetic basis of domestication and environmental adaptation of Triticeae crops. In this review, we focus on recent signs of progress in genome sequencing, pan-genome analyses, and resequencing analysis of Triticeae crops. We also propose future research avenues in Triticeae crop genomes, including identifying genome structure variations, the association of genomic regions with desired traits, mining functions of the non-coding area, introgression of high-quality genes from wild Triticeae resources, genome editing, and integration of genomic resources.

Development and Characterization of Near-Isogenic Lines Derived from Synthetic Wheat Revealing the 2 kb Insertion in the PPD-D1 Gene Responsible for Heading Delay and Grain Number Improvement

Spikelet number and grain number per spike are two crucial and correlated traits for grain yield in wheat. Photoperiod-1 (Ppd-1) is a key regulator of inflorescence architecture and spikelet formation in wheat. In this study, near-isogenic lines derived from the cross of a synthetic hexaploid wheat and commercial cultivars generated by double top-cross and two-phase selection were evaluated for the number of days to heading and other agronomic traits. The results showed that heading time segregation was conferred by a single incomplete dominant gene PPD-D1, and the 2 kb insertion in the promoter region was responsible for the delay in heading. Meanwhile, slightly delayed heading plants and later heading plants obviously have advantages in grain number and spikelet number of the main spike compared with early heading plants. Utilization of PPD-D1 photoperiod sensitivity phenotype as a potential means to increase wheat yield potential.

Genetic basis of geographical differentiation and breeding selection for wheat plant architecture traits

BACKGROUND

Plant architecture associated with increased grain yield and adaptation to the local environments is selected during wheat (Triticum aestivum) breeding. The internode length of individual stems and tiller length of individual plants are important for the determination of plant architecture. However, few studies have explored the genetic basis of these traits.

RESULTS

Here, we conduct a genome-wide association study (GWAS) to dissect the genetic basis of geographical differentiation of these traits in 306 worldwide wheat accessions including both landraces and traditional varieties. We determine the changes of haplotypes for the associated genomic regions in frequency in 831 wheat accessions that are either introduced from other countries or developed in China from last two decades. We identify 83 loci that are associated with one trait, while the remaining 247 loci are pleiotropic. We also find 163 associated loci are under strong selective sweep. GWAS results demonstrate independent regulation of internode length of individual stems and consistent regulation of tiller length of individual plants. This makes it possible to obtain ideal haplotype combinations of the length of four internodes. We also find that the geographical distribution of the haplotypes explains the observed differences in internode length among the worldwide wheat accessions.

CONCLUSION

This study provides insights into the genetic basis of plant architecture. It will facilitate gene functional analysis and molecular design of plant architecture for breeding.