Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

An Accurate Representation of the Number of bZIP Transcription Factors in the Triticum aestivum (Wheat) Genome and the Regulation of Functional Genes during Salt Stress

Climate change is dramatically increasing the overall area of saline soils around the world, which is increasing by approximately two million hectares each year. Soil salinity decreases crop yields and, thereby, makes farming less profitable, potentially causing increased poverty and hunger in many areas. A solution to this problem is increasing the salt tolerance of crop plants. Transcription factors (TFs) within crop plants represent a key to understanding salt tolerance, as these proteins play important roles in the regulation of functional genes linked to salt stress. The basic leucine zipper (bZIP) TF has a well-documented role in the regulation of salt tolerance. To better understand how bZIP TFs are linked to salt tolerance, we performed a genome-wide analysis in wheat using the Chinese spring wheat genome, which has been assembled by the International Wheat Genome Sequencing Consortium. We identified 89 additional bZIP gene sequences, which brings the total of bZIP gene sequences in wheat to 237. The majority of these 237 sequences included a single bZIP protein domain; however, different combinations of five other domains also exist. The bZIP proteins are divided into ten subfamily groups. Using an in silico analysis, we identified five bZIP genes (ABF2, ABF4, ABI5, EMBP1, and VIP1) that were involved in regulating salt stress. By scrutinizing the binding properties to the 2000 bp upstream region, we identified putative functional genes under the regulation of these TFs. Expression analyses of plant tissue that had been treated with or without 100 mM NaCl revealed variable patterns between the TFs and functional genes. For example, an increased expression of ABF4 was correlated with an increased expression of the corresponding functional genes in both root and shoot tissues, whereas VIP1 downregulation in root tissues strongly decreased the expression of two functional genes. Identifying strategies to sustain the expression of the functional genes described in this study could enhance wheat's salt tolerance.

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Wheat powdery mildew resistance gene Pm13 encodes a mixed lineage kinase domain-like protein

Wheat powdery mildew is one of the most destructive diseases threatening global wheat production. The wild relatives of wheat constitute rich sources of diversity for powdery mildew resistance. Here, we report the map-based cloning of the powdery mildew resistance gene Pm13 from the wild wheat species Aegilops longissima. Pm13 encodes a mixed lineage kinase domain-like (MLKL) protein that contains an N-terminal-domain of MLKL (MLKL_NTD) domain in its N-terminus and a C-terminal serine/threonine kinase (STK) domain. The resistance function of Pm13 is validated by mutagenesis, gene silencing, transgenic assay, and allelic association analyses. The development of introgression lines with significantly reduced chromosome segments of Ae. longissima encompassing Pm13 enables widespread deployment of this gene into wheat cultivars. The cloning of Pm13 may provide valuable insights into the molecular mechanisms underlying Pm13-mediated powdery mildew resistance and highlight the important roles of kinase fusion proteins (KFPs) in wheat immunity.

A Glu209Lys substitution in DRG1/TaACT7, which disturbs F-actin organization, reduces plant height and grain length in bread wheat

Plant height and grain size are two important agronomic traits that are closely related to crop yield. Numerous dwarf and grain-shape mutants have been studied to identify genes that can be used to increase crop yield and improve breeding programs. In this study, we characterized a dominant mutant, dwarf and round grain 1 (drg1-D), in bread wheat (Triticum aestivum L.). drg1-D plants exhibit multiple phenotypic changes, including dwarfism, round grains, and insensitivity to brassinosteroids (BR). Cell structure observation in drg1-D mutant plants showed that the reduced organ size is due to irregular cell shape. Using map-based cloning and verification in transgenic plants, we found that a Glu209Lys substitution in the DRG1 protein is responsible for the irregular cell size and arrangement in the drg1-D mutant. DRG1/TaACT7 encodes an actin family protein that is essential for polymerization stability and microfilament (MF) formation. In addition, the BR response and vesicular transport were altered by the abnormal actin cytoskeleton in drg1-D mutant plants. Our study demonstrates that DRG1/TaACT7 plays an important role in wheat cell shape determination by modulating actin organization and intracellular material transport, which could in the longer term provide tools to better understand the polymerization of actin and its assembly into filaments and arrays.

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.

Fine mapping of two recessive powdery mildew resistance genes from Aegilops tauschii accession CIae8

KEY MESSAGE

Two recessive powdery mildew resistance loci pmAeCIae8_2DS and pmAeCIae8_7DS from Aegilops tauschii were mapped and two synthesized hexaploid wheat lines were developed by distant hybridization. Wheat powdery mildew (Pm), one of the worldwide destructive fungal diseases, causes significant yield loss up to 30%. The identification of new Pm resistance genes will enrich the genetic diversity of wheat breeding for Pm resistance. Aegilops tauschii is the ancestor donor of sub-genome D of hexaploid wheat. It provides beneficial genes that can be easily transferred into wheat by producing synthetic hexaploid wheat followed by genetic recombination. We assessed the Pm resistance level of 35 Ae. tauschii accessions from different origins. Accession CIae8 exhibited high Pm resistance. Inheritance analysis and gene mapping were performed using F2 and F2:3 populations derived from the cross between CIae8 and a Pm susceptible accession PI574467. The Pm resistance of CIae8 was controlled by two independent recessive genes. Bulked segregate analysis using a 55 K SNP array revealed the SNPs were mainly enriched into genome regions, i.e. 2DS (13.5-20 Mb) and 7DS (4.0-15.5 Mb). The Pm resistance loci were named as pmAeCIae8_2DS and pmAeCIae8_7DS, respectively. By recombinant screening, we narrowed the pmAeCIae8_2DS into a 370-kb interval flanked by markers CINAU-AE7800 (14.89 Mb) and CINAU-AE20 (15.26 Mb), and narrowed the pmAeCIae8_7DS into a 260-kb interval flanked by markers CINAU-AE58 (4.72 Mb) and CINAU-AE25 (4.98 Mb). The molecular markers closely linked with the resistance loci were developed, and two synthesized hexaploid wheat (SHW) lines were produced. These laid the foundation for cloning of the two resistance loci and for transferring the resistance into common wheat.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15 Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T. aestivum cultivars (6x = AABBDD) and a single genome for Triticum durum (4x = AABB), Triticum dicoccoides (4x = AABB), Triticum urartu (2x = AA), and Aegilops tauschii (2x = DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

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.

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.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

Transcriptome Analysis of Roots from Wheat (Triticum aestivum L.) Varieties in Response to Drought Stress

Under climate change, drought is one of the most limiting factors that influences wheat (Triticum aestivum L.) production. Exploring stress-related genes is vital for wheat breeding. To identify genes related to the drought tolerance response, two common wheat cultivars, Zhengmai 366 (ZM366) and Chuanmai 42 (CM42), were selected based on their obvious difference in root length under 15% PEG-6000 treatment. The root length of the ZM366 cultivar was significantly longer than that of CM42. Stress-related genes were identified by RNA-seq in samples treated with 15% PEG-6000 for 7 days. In total, 11,083 differentially expressed genes (DEGs) and numerous single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) were identified. GO enrichment analysis revealed that the upregulated genes were mainly related to the response to water, acidic chemicals, oxygen-containing compounds, inorganic substances, and abiotic stimuli. Among the DEGs, the expression levels of 16 genes in ZM366 were higher than those in CM42 after the 15% PEG-6000 treatment based on RT-qPCR. Furthermore, EMS-induced mutants in Kronos (T. turgidum L.) of 4 representative DEGs possessed longer roots than the WT after the 15% PEG-6000 treatment. Altogether, the drought stress genes identified in this study represent useful gene resources for wheat breeding.

Breeding and Genomic Approaches towards Development of Fusarium Wilt Resistance in Chickpea

Chickpea is an important leguminous crop with potential to provide dietary proteins to both humans and animals. It also ameliorates soil nitrogen through biological nitrogen fixation. The crop is affected by an array of biotic and abiotic factors. Among different biotic stresses, a major fungal disease called Fusarium wilt, caused by Fusarium oxysporum f. sp. ciceris (FOC), is responsible for low productivity in chickpea. To date, eight pathogenic races of FOC (race 0, 1A, and 1B/C, 2-6) have been reported worldwide. The development of resistant cultivars using different conventional breeding methods is very time consuming and depends upon the environment. Modern technologies can improve conventional methods to solve these major constraints. Understanding the molecular response of chickpea to Fusarium wilt can help to provide effective management strategies. The identification of molecular markers closely linked to genes/QTLs has provided great potential for chickpea improvement programs. Moreover, omics approaches, including transcriptomics, metabolomics, and proteomics give scientists a vast viewpoint of functional genomics. In this review, we will discuss the integration of all available strategies and provide comprehensive knowledge about chickpea plant defense against Fusarium wilt.

Exon shuffling potentiates a diverse repertoire of brown algal NB-ARC-TPR candidate immune receptor proteins via alternative splicing

Like other organisms, brown algae are subject to diseases caused by bacteria, fungi, and viruses. Brown algal immunity mechanisms are not well characterized; however, there is evidence suggesting that pathogen receptors exist in brown algae. One key protein family likely associated with brown algal innate immunity possesses an NB-ARC domain analogous to innate immune proteins in plants and animals. In this study, we conducted an extensive survey of NB-ARC genes in brown algae and obtained insights into the domain organization and evolutionary history of the encoded proteins. Our data show that brown algae possess an ancient NB-ARC-tetratricopeptide repeat (NB-TPR) domain architecture. We identified an N-terminal effector domain, the four-helix bundle, which was not previously found associated with NB-ARC domains. The phylogenetic tree including NB-ARC domains from all kingdoms of life suggests the three clades of brown algal NB-TPRs are likely monophyletic, whereas their TPRs seem to have distinct origins. One group of TPRs exhibit intense exon shuffling, with various alternative splicing and diversifying selection acting on them, suggesting exon shuffling is an important mechanism for evolving ligand-binding specificities. The reconciliation of gene duplication and loss events of the NB-ARC genes reveals that more independent gene gains than losses have occurred during brown algal evolution, and that tandem duplication has played a major role in the expansion of NB-ARC genes. Our results substantially enhance our understanding of the evolutionary history and exon shuffling mechanisms of the candidate innate immune repertoire of brown algae.

The wheat leaf delayed virescence of mutant dv4 is associated with the abnormal photosynthetic and antioxidant systems

Chlorophyll (Chl) is a key pigment for wheat (Triticum aestivum L.) photosynthesis, consequently impacts grain yield. A wheat mutant named as delayed virescence 4 (dv4) was obtained from cultivar Guomai 301 (wild type, WT) treated with ethyl methane sulfonate (EMS). The seedling leaves of dv4 were shallow yellow, apparently were chlorophyll deficient. They started to turn green at the jointing stage and returned to almost ordinary green at the heading stage. Leaf transcriptome comparison of Guomai 301 and dv4 at the jointing stage showed that most differentially expressed genes (DEGs) of transcription and translation were highly expressed in dv4, one key gene nicotianamine aminotransferase A (NAAT-A) involved in the synthesis and metabolism pathways of tyrosine, methionine and phenylalanine was significantly lowly expressed. The expression levels of the most photosynthesis related genes, such as photosystem I (PS I), ATPase and light-harvesting chlorophyll protein complex-related homeotypic genes, and protochlorophyllide reductase A (PORA) were lower; but macromolecule degradation and hypersensitivity response (HR) related gene heat shock protein 82 (HSP82) was highly expressed. Compared to WT, the contents of macromolecules such as proteins and sugars were reduced; the contents of Chl a, Chl b, total Chl, and carotenoids in leaves of dv4 were significantly less at the jointing stage, while the ratio of Chl a / Chl b was the same as that of WT. The net photosynthetic rate, stomatal conductance and transpiration rate of dv4 were significantly lower. The H₂O₂ content were higher, while the contents of total phenol and malondialdehyde (MDA), antioxidant enzyme activities were lower in leaves of dv4. In conclusion, the reduced contents of macromolecules and photosynthetic pigments, the abnormal photosynthetic and antioxidant systems were closely related to the phenotype of dv4.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

Advances in the Mining of Disease Resistance Genes from Aegilops tauschii and the Utilization in Wheat

Aegilops tauschii is one of the malignant weeds that affect wheat production and is also the wild species ancestor of the D genome of hexaploid wheat (Triticum aestivum, AABBDD). It contains many disease resistance genes that have been lost in the long-term evolution of wheat and is an important genetic resource for the mining and utilization of wheat disease resistance genes. In recent years, the genome sequence of Aegilops tauschii has been preliminarily completed, which has laid a good foundation for the further exploration of wheat disease resistance genes in Aegilops tauschii. There are many studies on disease resistance genes in Aegilops tauschii; in order to provide better help for the disease resistance breeding of wheat, this paper analyzes and reviews the relationship between Aegilops tauschii and wheat, the research progress of Aegilops tauschii, the discovery of disease resistance genes from Aegilops tauschii, and the application of disease resistance genes from Aegilops tauschii to modern wheat breeding, providing a reference for the further exploration and utilization of Aegilops tauschii in wheat disease resistance breeding.

Specific-Locus Amplified Fragment Sequencing (SLAF-Seq)

Specific length amplified fragment sequencing (SLAF-seq) technology is a simplified genome sequencing technology based on next-generation sequencing. SLAF-seq technology has several distinguishing characteristics: 1. Deep sequencing to ensure accuracy of genotyping; 2. Effectively reduce sequencing costs; 3. Pre-designed simplified representation scheme to optimize marker efficiency; 4. Doubled barcode system for large populations. The advantages and technical process of SLAF-seq are described briefly with summarized results for the application of SLAF-seq in development of molecular markers, construction of high-density genetic map and gene mapping in ornamental plants. Finally, the difficulties and prospects of this method are discussed in application.

Genetic Improvement and Application Practices of Synthetic Hexaploid Wheat

Synthetic hexaploid wheat (SHW) is a useful genetic resource that can be used to improve the performance of common wheat by transferring favorable genes from a wide range of tetraploid or diploid donors. From the perspectives of physiology, cultivation, and molecular genetics, the use of SHW has the potential to increase wheat yield. Moreover, genomic variation and recombination were enhanced in newly formed SHW, which could generate more genovariation or new gene combinations compared to ancestral genomes. Accordingly, we presented a breeding strategy for the application of SHW-the 'large population with limited backcrossing method'-and we pyramided stripe rust resistance and big-spike-related QTLs/genes from SHW into new high-yield cultivars, which represents an important genetic basis of big-spike wheat in southwestern China. For further breeding applications of SHW-derived cultivars, we used the 'recombinant inbred line-based breeding method' that combines both phenotypic and genotypic evaluations to pyramid multi-spike and pre-harvest sprouting resistance QTLs/genes from other germplasms to SHW-derived cultivars; consequently, we created record-breaking high-yield wheat in southwestern China. To meet upcoming environmental challenges and continuous global demand for wheat production, SHW with broad genetic resources from wild donor species will play a major role in wheat breeding.

Mapping and characterization of a novel adult-plant leaf rust resistance gene LrYang16G216 via bulked segregant analysis and conventional linkage method

A novel adult-plant leaf rust resistance gene LrYang16G216 on wheat chromosome 6BL was identified and mapped to a 0.59 cM genetic interval by BSA and conventional linkage method. Leaf rust (Puccinia triticina) is one of the most devastating fungal diseases of wheat (Triticum aestivum L.). Discovery and identification of new resistance genes is essential to develop disease-resistant cultivars. An advanced breeding line Yang16G216 was previously identified to confer adult-plant resistance (APR) to leaf rust. In this research, a recombinant inbred line (RIL) population was constructed from the cross between Yang16G216 and a highly susceptible line Yang16M6393, and genotyped with exome capture sequencing and 55 K SNP array. Through bulked segregant analysis (BSA) and genetic linkage mapping, a stable APR gene, designated as LrYang16G216, was detected and mapped to the distal region of chromosome arm 6BL with a genetic interval of 2.8 cM. For further verification, another RIL population derived from the cross between Yang16G216 and a susceptible wheat variety Yangmai 29 was analyzed using the enriched markers in the target interval, and LrYang16G216 was further narrowed to a 0.59 cM genetic interval flanked by the KASP markers Ax109403980 and Ax95083494, corresponding to the physical position 712.34-713.94 Mb in the Chinese Spring reference genome, in which twenty-six disease resistance-related genes were annotated. Based on leaf rust resistance spectrum, mapping data and physical location, LrYang16G216 was identified to be a novel and effective APR gene. The LrYang16G216 with linked markers will be useful for marker-assisted selection in wheat resistance breeding.

Involvement of miRNAs in Metabolic Herbicide Resistance to Bispyribac-Sodium in Echinochloa crus-galli (L.) P. Beauv

Several mechanisms involved in weed herbicide resistance are unknown, particularly those acting at the epigenetic level, such as the capacity of small-non-coding RNAs (sncRNAs) to target messenger RNAs of genes involved in herbicide detoxification. The transcription of these sncRNAs is stimulated by epigenetic factors, thereby affecting gene expression. This study was carried out in order to evaluate, for the first time in Echinochloa crus-galli (L.) P. Beauv. (barnyardgrass), the capacity of miRNAs to regulate the expression of genes associated with bispyribac-sodium detoxification. The expression profiles of eight miRNAs with a high degree of complementarity (≥80%) with mRNAs of genes involved in herbicide detoxification (CYP450, GST and eIF4B) were determined by qRT-PCR before and after herbicide spraying. Five of the miRNAs studied (gra-miR7487c, gma-miR396f, gra-miR8759, osa-miR395f, ath-miR847) showed an increased expression after herbicide application in both susceptible and resistant biotypes. All the miRNAs, except gra-miR8759, were more highly expressed in the herbicide-resistant biotypes. In specimens with increased expression of miRNAs, we observed reduced expression of the target genes. The remaining three miRNAs (ata-miR166c-5p, ath-miR396b-5p and osa-miR5538) showed no over-expression after herbicide treatment, and no difference in expression was recorded between susceptible and resistant biotypes. Our results represent a first overview of the capacity of miRNAs to regulate the expression of genes involved in bispyribac-sodium detoxification in the genus Echinochloa. Further research is required to identify novel miRNAs and target genes to develop more focused and sustainable strategies of weed control.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

The genome and gene editing system of sea barleygrass provide a novel platform for cereal domestication and stress tolerance studies

The tribe Triticeae provides important staple cereal crops and contains elite wild species with wide genetic diversity and high tolerance to abiotic stresses. Sea barleygrass (Hordeum marinum Huds.), a wild Triticeae species, thrives in saline marshlands and is well known for its high tolerance to salinity and waterlogging. Here, a 3.82-Gb high-quality reference genome of sea barleygrass is assembled de novo, with 3.69 Gb (96.8%) of its sequences anchored onto seven chromosomes. In total, 41 045 high-confidence (HC) genes are annotated by homology, de novo prediction, and transcriptome analysis. Phylogenetics, non-synonymous/synonymous mutation ratios (Ka/Ks), and transcriptomic and functional analyses provide genetic evidence for the divergence in morphology and salt tolerance among sea barleygrass, barley, and wheat. The large variation in post-domestication genes (e.g. IPA1 and MOC1) may cause interspecies differences in plant morphology. The extremely high salt tolerance of sea barleygrass is mainly attributed to low Na⁺ uptake and root-to-shoot translocation, which are mainly controlled by SOS1, HKT, and NHX transporters. Agrobacterium-mediated transformation and CRISPR/Cas9-mediated gene editing systems were developed for sea barleygrass to promote its utilization for exploration and functional studies of hub genes and for the genetic improvement of cereal crops.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

Epistatic interaction effect between chromosome 1BL (Yr29) and a novel locus on 2AL facilitating resistance to stripe rust in Chinese wheat Changwu 357-9

Four stable QTL for adult plant resistance were identified in wheat line Changwu 357-9, including a new QTL on 2AL showing significant interaction with Yr29 to reduce stripe rust severity. Stripe rust (yellow rust) is a serious disease of bread wheat (Triticum aestivum L.) worldwide. Genetic resistance is considered the most economical, effective and environmentally friendly method to control the disease and to minimize the use of fungicides. The current study focused on characterizing the components of stripe rust resistance and understanding the interactions in Changwu 357-9 (CW357-9)/Avocet S RIL population. A genetic linkage map constructed using a new GenoBaits Wheat 16K Panel and the 660K SNP array had 5104 polymorphic SNP markers spanning 3533.11 cM. Four stable QTL, consistently identified across five environments, were detected on chromosome arms 1BL, 2AL, 3DS, and 6BS in Changwu357-9. The most effective QTL QYrCW357-1BL was Yr29. The 6BS QTL was identified as Yr78, which has been combined with the 1BL QTL in many wheat cultivars and breeding lines. The novel QTL on 2AL with moderate effect showed a stable and significant epistatic interaction with Yr29. The QTL on 3DL should be same as QYrsn.nwafu-3DL and enriches the overall stripe rust resistance gene pool for breeding. Polymorphisms of flanking AQP markers AX-110020417 (for QYrCW357-1BL), AX-110974948 (for QYrCW357-2AL), AX-109466386 (for QYrCW357-3DL), and AX-109995005 (for QYrCW357-6BS) were evaluated in a diversity panel including 225 wheat cultivars and breeding lines. These results suggested that these high-throughput markers could be used to introduce QYrCW357-1BL, QYrCW357-2AL, QYrCW357-3DL, and QYrCW357-6BS into commercial wheat cultivars. Combinations of these genes with other APR QTL should lead to higher levels of stripe rust resistance along with the beneficial effects of multi-disease resistance gene Yr29 on improving resistance to other diseases.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Whole-genome resequencing of the wheat A subgenome progenitor Triticum urartu provides insights into its demographic history and geographic adaptation

Triticum urartu is the progenitor of the A subgenome in tetraploid and hexaploid wheat. Uncovering the landscape of genetic variations in T. urartu will help us understand the evolutionary and polyploid characteristics of wheat. Here, we investigated the population genomics of T. urartu by genome-wide sequencing of 59 representative accessions collected around the world. A total of 42.2 million high-quality single-nucleotide polymorphisms and 3 million insertions and deletions were obtained by mapping reads to the reference genome. The ancient T. urartu population experienced a significant reduction in effective population size (Ne) from ∼3 000 000 to ∼140 000 and subsequently split into eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations during the Younger Dryas period. A map of allelic drift paths displayed splits and mixtures between different geographic groups, and a strong genetic drift towards hexaploid wheat was also observed, indicating that the direct donor of the A subgenome originated from northwestern Syria. Genetic changes were revealed between the eastern Mediterranean coastal and Mesopotamian-Transcaucasian populations in genes orthologous to those regulating plant development and stress responses. A genome-wide association study identified two single-nucleotide polymorphisms in the exonic regions of the SEMI-DWARF 37 ortholog that corresponded to the different T. urartu ecotype groups. Our study provides novel insights into the origin and genetic legacy of the A subgenome in polyploid wheat and contributes a gene repertoire for genomics-enabled improvements in wheat breeding.

Metabolome dynamics during wheat domestication

One of the most important crops worldwide is wheat. Wheat domestication took place about 10,000 years ago. Not only that its wild progenitors have been discovered and phenotypically characterized, but their genomes were also sequenced and compared to modern wheat. While comparative genomics is essential to track genes that contribute to improvement in crop yield, comparative analyses of functional biological end-products, such as metabolites, are still lacking. With the advent of rigorous mass-spectrometry technologies, it is now possible to address that problem on a big-data scale. In attempt to reveal classes of metabolites, which are associated with wheat domestication, we analyzed the metabolomes of wheat kernel samples from various wheat lines. These wheat lines represented subspecies of tetraploid wheat along primary and secondary domestications, including wild emmer, domesticated emmer, landraces durum, and modern durum. We detected that the groups of plant metabolites such as plant-defense metabolites, antioxidants and plant hormones underwent significant changes during wheat domestication. Our data suggest that these metabolites may have contributed to the improvement in the agricultural fitness of wheat. Closer evaluation of specific metabolic pathways may result in the future in genetically-engineered high-yield crops.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

A new strategy for using historical imbalanced yield data to conduct genome-wide association studies and develop genomic prediction models for wheat breeding

Using imbalanced historical yield data to predict performance and select new lines is an arduous breeding task. Genome-wide association studies (GWAS) and high throughput genotyping based on sequencing techniques can increase prediction accuracy. An association mapping panel of 227 Texas elite (TXE) wheat breeding lines was used for GWAS and a training population to develop prediction models for grain yield selection. An imbalanced set of yield data collected from 102 environments (year-by-location) over 10 years, through testing yield in 40-66 lines each year at 6-14 locations with 38-41 lines repeated in the test in any two consecutive years, was used. Based on correlations among data from different environments within two adjacent years and heritability estimated in each environment, yield data from 87 environments were selected and assigned to two correlation-based groups. The yield best linear unbiased estimation (BLUE) from each group, along with reaction to greenbug and Hessian fly in each line, was used for GWAS to reveal genomic regions associated with yield and insect resistance. A total of 74 genomic regions were associated with grain yield and two of them were commonly detected in both correlation-based groups. Greenbug resistance in TXE lines was mainly controlled by Gb3 on chromosome 7DL in addition to two novel regions on 3DL and 6DS, and Hessian fly resistance was conferred by the region on 1AS. Genomic prediction models developed in two correlation-based groups were validated using a set of 105 new advanced breeding lines and the model from correlation-based group G2 was more reliable for prediction. This research not only identified genomic regions associated with yield and insect resistance but also established the method of using historical imbalanced breeding data to develop a genomic prediction model for crop improvement.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01287-8.

Wheat breeding history reveals synergistic selection of pleiotropic genomic sites for plant architecture and grain yield

Diversity surveys of crop germplasm are important for gaining insights into the genomic basis for plant architecture and grain yield improvement, which is still poorly understood in wheat. In this study, we exome sequenced 287 wheat accessions that were collected in the past 100 years. Population genetics analysis identified that 6.7% of the wheat genome falls within the selective sweeps between landraces and cultivars, which harbors the genes known for yield improvement. These regions were asymmetrically distributed on the A and B subgenomes with regulatory genes being favorably selected. Genome-wide association study (GWAS) identified genomic loci associated with traits for yield potential, and two underlying genes, TaARF12 encoding an auxin response factor and TaDEP1 encoding the G-protein γ-subunit, were located and characterized to pleiotropically regulate both plant height and grain weight. Elite single-nucleotide haplotypes with increased allele frequency in cultivars relative to the landraces were identified and found to have accumulated over the course of breeding. Interestingly, we found that TaARF12 and TaDEP1 function in epistasis with the classical plant height Rht-1 locus, leading to propose a "Green Revolution"-based working model for historical wheat breeding. Collectively, our study identifies selection signatures that fine-tune the gibberellin pathway during modern wheat breeding and provides a wealth of genomic diversity resources for the wheat research community.

Characterization of TaSPP-5A gene associated with sucrose content in wheat (Triticum aestivum L.)

BACKGROUND

Sucrose, the major product of photosynthesis and the primary sugar transported as a soluble carbohydrate via the phloem, is a critical determinant for harvest yield in wheat crops. Sucrose-phosphatase (SPP) catalyzes the final step in the sucrose biosynthesis pathway, implying its essential role in the plant.

RESULT

In this study, wheat SPP homologs genes were isolated from chromosomes 5A, 5B, and 5D, designated as TaSPP-5A, TaSPP-5B, and TaSPP-5D, respectively. Sequence alignment showed one 1-bp Insertion-deletion (InDel) and three single nucleotide polymorphisms (SNPs) at TaSPP-5A coding region, forming two haplotypes, TaSPP-5Aa and TaSPP-5Ab, respectively. A derived cleaved amplified polymorphism sequence (dCAPS) marker, TaSPP-5A-dCAPS, was developed to discriminate allelic variation based on the polymorphism at position 1242 (C-T). A total of 158 varieties were used to perform a TaSPP-5A marker-trait association analysis, where two haplotypes were significantly associated with sucrose content in two environments and with thousand-grain weight (TGW) and grain length (GL) in three environments. Quantitative real-time PCR further revealed that TaSPP-5Aa showed relatively higher expression than TaSPP-5Ab in wheat seedling leaves, generally associating with increased sucrose content and TGW. The expression of TaSPP-5A and sucrose content in TaSPP-5Aa haplotypes were also higher than those in TaSPP-5Ab haplotypes under both 20% PEG-6000 and 100 μM ABA treatment. Sequence alignment showed that the two TaSPP-5A haplotypes comprised 11 SNPs from -395 to -1962 bp at TaSPP-5A promoter locus, participating in the formation of several conserved sequences, may account for the high expression of TaSPP-5A in TaSPP-5Aa haplotypes. In addition, the distribution analysis of TaSPP-5A haplotypes revealed that TaSPP-5Aa was preferred in the natural wheat population, being strongly positively selected in breeding programs.

CONCLUSION

According to the SNPs detected in the TaSPP-5A sequence, two haplotypes, TaSPP-5Aa and TaSPP-5Ab, were identified among wheat accessions, which potential value for sucrose content selection was validated by association analysis. Our results indicate that the favorable allelic variation TaSPP-5Aa should be valuable in enhancing grain yield by improving the sucrose content. Furthermore, a functional marker, TaSPP-5A-dCAPS, can be used for marker-assisted selection to improve grain weight in wheat and provides insights into the biological function of TaSPP-5A gene.

Significance of the Diversification of Wheat Species for the Assembly and Functioning of the Root-Associated Microbiome

Wheat, one of the major crops in the world, has had a complex history that includes genomic hybridizations between Triticum and Aegilops species and several domestication events, which resulted in various wild and domesticated species (especially Triticum aestivum and Triticum durum), many of them still existing today. The large body of information available on wheat-microbe interactions, however, was mostly obtained without considering the importance of wheat evolutionary history and its consequences for wheat microbial ecology. This review addresses our current understanding of the microbiome of wheat root and rhizosphere in light of the information available on pre- and post-domestication wheat history, including differences between wild and domesticated wheats, ancient and modern types of cultivars as well as individual cultivars within a given wheat species. This analysis highlighted two major trends. First, most data deal with the taxonomic diversity rather than the microbial functioning of root-associated wheat microbiota, with so far a bias toward bacteria and mycorrhizal fungi that will progressively attenuate thanks to the inclusion of markers encompassing other micro-eukaryotes and archaea. Second, the comparison of wheat genotypes has mostly focused on the comparison of T. aestivum cultivars, sometimes with little consideration for their particular genetic and physiological traits. It is expected that the development of current sequencing technologies will enable to revisit the diversity of the wheat microbiome. This will provide a renewed opportunity to better understand the significance of wheat evolutionary history, and also to obtain the baseline information needed to develop microbiome-based breeding strategies for sustainable wheat farming.

Assessing the Potential of Using the Langdon 5D(5B) Substitution Line for the Introgression of Aegilops tauschii Into Durum Wheat

Aegilops tauschii, the D-genome donor of hexaploid wheat, provides a source of genetic variation that could be used for tetraploid (durum) wheat improvement. In addition to the genes for wheat quality on the D-genome, which differentiate between bread and durum wheats in terms of end-use properties, genes coding for resistances to biotic and abiotic stresses are also present on the D-genome which would be useful in durum wheat. The introgression of Ae. tauschii into durum wheat, however, requires cytogenetic manipulation to induce homoeologous chromosome pairing to promote recombination. For this purpose, the introgression of Ae. tauschii into durum wheat was performed through a bridge cross of the wild species to the Langdon 5D(5B) disomic substitution line that lacks the Ph1 locus present on chromosome 5B, followed by a cross of the F₁ to the durum wheat cultivar Om Rabi 5. Subsequent generations were self-fertilized, and these were screened for D-genome introgressions using (i) D-genome-specific Kompetitive Allele-Specific PCR (KASP) markers and (ii) KASP markers polymorphic between the 5D chromosomes of wheat, present in the Langdon 5D(5B) substitution line, and of Ae. tauschii. Homozygous introgression lines were confirmed using genomic and fluorescence in situ hybridization. The results showed that the use of the Langdon 5D(5B) disomic substitution line did not promote D-genome introgression across all linkage groups with only a limited success in the introgression of Ae. tauschii 5D segments into durum wheat.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.

Differentially Amplified Repetitive Sequences Among Aegilops tauschii Subspecies and Genotypes

Genomic repetitive sequences commonly show species-specific sequence type, abundance, and distribution patterns, however, their intraspecific characteristics have been poorly described. We quantified the genomic repetitive sequences and performed single nucleotide polymorphism (SNP) analysis between 29 Ae. tauschii genotypes and subspecies using publicly available raw genomic Illumina sequence reads and used fluorescence in situ hybridization (FISH) to experimentally analyze some repeats. The majority of the identified repetitive sequences had similar contents and proportions between anathera, meyeri, and strangulata subspecies. However, two Ty3/gypsy retrotransposons (CL62 and CL87) showed significantly higher abundances, and CL1, CL119, CL213, CL217 tandem repeats, and CL142 retrotransposon (Ty1/copia type) showed significantly lower abundances in subspecies strangulata compared with the subspecies anathera and meyeri. One tandem repeat and 45S ribosomal DNA (45S rDNA) abundances showed a high variation between genotypes but their abundances were not subspecies specific. Phylogenetic analysis using the repeat abundances of the aforementioned clusters placed the strangulata subsp. in a distinct clade but could not discriminate anathera and meyeri. A near complete differentiation of anathera and strangulata subspecies was observed using SNP analysis; however, var. meyeri showed higher genetic diversity. FISH using major tandem repeats couldn't detect differences between subspecies, although (GAA)₁₀ signal patterns generated two different karyotype groups. Taken together, the different classes of repetitive DNA sequences have differentially accumulated between strangulata and the other two subspecies of Ae. tauschii that is generally in agreement with spike morphology, implying that factors affecting repeatome evolution are variable even among highly closely related lineages.

Long-range assembly of sequences helps to unravel the genome structure and small variation of the wheat-Haynaldia villosa translocated chromosome 6VS.6AL

Genomics studies in wild species of wheat have been limited due to the lack of references; however, new technologies and bioinformatics tools have much potential to promote genomic research. The wheat-Haynaldia villosa translocation line T6VS·6AL has been widely used as a backbone parent of wheat breeding in China. Therefore, revealing the genome structure of translocation chromosome 6VS·6AL will clarify how this chromosome formed and will help to determine how it affects agronomic traits. In this study, chromosome flow sorting, NGS sequencing and Chicago long-range linkage assembly were innovatively used to produce the assembled sequences of 6VS·6AL, and gene prediction and genome structure characterization at the molecular level were effectively performed. The analysis discovered that the short arm of 6VS·6AL was actually composed of a large distal segment of 6VS, a small proximal segment of 6AS and the centromere of 6A, while the collinear region in 6VS corresponding to 230-260 Mb of 6AS-Ta was deleted when the recombination between 6VS and 6AS occurred. In addition to the molecular mechanism of the increased grain weight and enhanced spike length produced by the translocation chromosome, it may be correlated with missing GW2-V and an evolved NRT-V cluster. Moreover, a fine physical bin map of 6VS was constructed by the high-throughput developed 6VS-specific InDel markers and a series of newly identified small fragment translocation lines involving 6VS. This study will provide essential information for mining of new alien genes carried by the 6VS·6AL translocation chromosome.

Characterization of Nucleotide Binding -Site-Encoding Genes in Sweetpotato, Ipomoea batatas(L.) Lam., and Their Response to Biotic and Abiotic Stresses

Sweetpotato, Ipomoea batatas (L.) Lam., is an important and widely grown crop, yet its production is affected severely by biotic and abiotic stresses. The nucleotide binding site (NBS)-encoding genes have been shown to improve stress tolerance in several plant species. However, the characterization of NBS-encoding genes in sweetpotato is not well-documented to date. In this study, a comprehensive analysis of NBS-encoding genes has been conducted on this species by using bioinformatics and molecular biology methods. A total of 315 NBS-encoding genes were identified, and 260 of them contained all essential conserved domains while 55 genes were truncated. Based on domain architectures, the 260 NBS-encoding genes were grouped into 6 distinct categories. Phylogenetic analysis grouped these genes into 3 classes: TIR, CC (I), and CC (II). Chromosome location analysis revealed that the distribution of NBS-encoding genes in chromosomes was uneven, with a number ranging from 1 to 34. Multiple stress-related regulatory elements were detected in the promoters, and the NBS-encoding genes' expression profiles under biotic and abiotic stresses were obtained. According to the bioinformatics analysis, 9 genes were selected for RT-qPCR analysis. The results revealed that IbNBS75, IbNBS219, and IbNBS256 respond to stem nematode infection; Ib-NBS240, IbNBS90, and IbNBS80 respond to cold stress, while IbNBS208, IbNBS71, and IbNBS159 respond to 30% PEG treatment. We hope these results will provide new insights into the evolution of NBS-encoding genes in the sweetpotato genome and contribute to the molecular breeding of sweetpotato in the future.

Present and future prospects for wheat improvement through genome editing and advanced technologies

Wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) is one of the most important staple food crops in the world. Despite the fact that wheat production has significantly increased over the past decades, future wheat production will face unprecedented challenges from global climate change, increasing world population, and water shortages in arid and semi-arid lands. Furthermore, excessive applications of diverse fertilizers and pesticides are exacerbating environmental pollution and ecological deterioration. To ensure global food and ecosystem security, it is essential to enhance the resilience of wheat production while minimizing environmental pollution through the use of cutting-edge technologies. However, the hexaploid genome and gene redundancy complicate advances in genetic research and precision gene modifications for wheat improvement, thus impeding the breeding of elite wheat cultivars. In this review, we first introduce state-of-the-art genome-editing technologies in crop plants, especially wheat, for both functional genomics and genetic improvement. We then outline applications of other technologies, such as GWAS, high-throughput genotyping and phenotyping, speed breeding, and synthetic biology, in wheat. Finally, we discuss existing challenges in wheat genome editing and future prospects for precision gene modifications using advanced genome-editing technologies. We conclude that the combination of genome editing and other molecular breeding strategies will greatly facilitate genetic improvement of wheat for sustainable global production.

High-quality chromosome-level genomes of Cucumis metuliferus and Cucumis melo provide insight into Cucumis genome evolution

Cucumis metuliferus (African horned cucumber), a wild relative of Cucumis sativus (cucumber) and Cucumis melo (melon), displays high-level resistance to several important plant pathogens (e.g., root-knot nematodes and several viruses). Here, we report a chromosome-level genome assembly for C. metuliferus, with a 316 Mb genome sequence comprising 29 039 genes. Phylogenetic analysis of related species in family Cucurbitaceae indicated that the divergence time between C. metuliferus and melon was 17.8 million years ago. Comparisons between the C. metuliferus and melon genomes revealed large structural variations (inversions and translocations >1 Mb) in eight chromosomes of these two species. Gene family comparison showed that C. metuliferus has the largest number of resistance-related nucleotide-binding site leucine-rich repeat (NBS-LRR) genes in Cucurbitaceae. The loss of NBS-LRR loci caused by large insertions or deletions (indels) and pseudogenization caused by small indels explained the loss of NBS-LRR genes in Cucurbitaceae. Population structure analysis suggested that C. metuliferus originated in Zimbabwe, then spread to other southern African regions where it likely underwent similar domestic selection as melon. This C. metuliferus reference sequence will accelerate the understanding of the molecular evolution of resistance-related genes and enhance cucurbit crop improvement efforts.

High-resolution mapping of the Mov-1 locus in wheat by combining radiation hybrid (RH) and recombination-based mapping approaches

This work reports a quick method that integrates RH mapping and genetic mapping to map the dominant Mov-1 locus to a 1.1-Mb physical interval with a small number of candidate genes. Bread wheat is an important crop for global human population. Identification of genes and alleles controlling agronomic traits is essential toward sustainably increasing crop production. The unique multi-ovary (MOV) trait in wheat holds potential for improving yields and is characterized by the formation of 2-3 grains per spikelet. The genetic basis of the multi-ovary trait is known to be monogenic and dominant in nature. Its precise mapping and functional characterization is critical to utilizing this trait in a feasible manner. Previous mapping efforts of the locus controlling multiple ovary/pistil formation in the hexaploid wheat have failed to produce a consensus for a particular chromosome. We describe a mapping strategy integrating radiation hybrid mapping and high-resolution genetic mapping to locate the chromosomal position of the Mov-1 locus in hexaploid wheat. We used RH mapping approach using a panel of 188 lines to map the Mov-1 locus in the terminal part of long arm of wheat chromosome 2D with a map resolution of 1.67 Mb/cR₁₅₀₀. Then using a genetic population of MOV × Synthetic wheat of F₂ lines, we delineated the Mov-1 locus to a 1.1-Mb physical region with a small number of candidate genes. This demonstrates the value of this integrated strategy to mapping dominant genes in wheat.

The global significance of biodiversity science in China: an overview

Biodiversity science in China has seen rapid growth over recent decades, ranging from baseline biodiversity studies to understanding the processes behind evolution across dynamic regions such as the Qinghai-Tibetan Plateau. We review research, including species catalogues; biodiversity monitoring; the origins, distributions, maintenance and threats to biodiversity; biodiversity-related ecosystem function and services; and species and ecosystems' responses to global change. Next, we identify priority topics and offer suggestions and priorities for future biodiversity research in China. These priorities include (i) the ecology and biogeography of the Qinghai-Tibetan Plateau and surrounding mountains, and that of subtropical and tropical forests across China; (ii) marine and inland aquatic biodiversity; and (iii) effective conservation and management to identify and maintain synergies between biodiversity and socio-economic development to fulfil China's vision for becoming an ecological civilization. In addition, we propose three future strategies: (i) translate advanced biodiversity science into practice for biodiversity conservation; (ii) strengthen capacity building and application of advanced technologies, including high-throughput sequencing, genomics and remote sensing; and (iii) strengthen and expand international collaborations. Based on the recent rapid progress of biodiversity research, China is well positioned to become a global leader in biodiversity research in the near future.

Distribution of MITE family Monkey King in rapeseed (Brassica napus L) and its influence on gene expression

Miniature inverted-repeat transposable elements (MITEs) are a group of class II transposable elements. The MITE Monkey King (MK) was first discovered upstream of BnFLC.A10. In this study, genome resequencing of four selected B. napus accessions, revealed more than 4000 distributed copies of MKs constituting ~2.4 Mb of the B. napus genomic sequence and caused 677 polymorphisms among the four accessions. MK -polymorphism-related markers across 128 natural and 58 synthetic accessions revealed more polymorphic MKs in natural than synthetic accessions. Ten MK -induced indels significantly affected the expression levels of the nearest gene based on RNAseq analysis, six of these effects were subsequently confirmed using qRT-PCR. Decreased expression pattern of MK -derived miRNA-bna-miR6031 was also observed under various stress treatments. Further research focused on the MITE families should promote not only our understanding of gene regulatory networks but also inform crop improvement efforts.

Introgressing the Aegilops tauschii genome into wheat as a basis for cereal improvement

Increasing crop production is necessary to feed the world's expanding population, and crop breeders often utilize genetic variations to improve crop yield and quality. However, the narrow diversity of the wheat D genome seriously restricts its selective breeding. A practical solution is to exploit the genomic variations of Aegilops tauschii via introgression. Here, we established a rapid introgression platform for transferring the overall genetic variations of A. tauschii to elite wheats, thereby enriching the wheat germplasm pool. To accelerate the process, we assembled four new reference genomes, resequenced 278 accessions of A. tauschii and constructed the variation landscape of this wheat progenitor species. Genome comparisons highlighted diverse functional genes or novel haplotypes with potential applications in wheat improvement. We constructed the core germplasm of A. tauschii, including 85 accessions covering more than 99% of the species' overall genetic variations. This was crossed with elite wheat cultivars to generate an A. tauschii-wheat synthetic octoploid wheat (A-WSOW) pool. Laboratory and field analysis with two examples of the introgression lines confirmed its great potential for wheat breeding. Our high-quality reference genomes, genomic variation landscape of A. tauschii and the A-WSOW pool provide valuable resources to facilitate gene discovery and breeding in wheat.

Novel Insights Into Genetic Responses for Waterlogging Stress in Two Local Wheat Cultivars in Yangtze River Basin

Wheat (Triticum aestivum L.), the most widely cultivated crop, is affected by waterlogging that limited the wheat production. Given the incompleteness of its genome annotation, PacBio SMRT plus Illumina short-read sequencing strategy provided an efficient approach to investigate the genetic regulation of waterlogging stress in wheat. A total of 947,505 full-length non-chimetric (FLNC) sequences were obtained with two wheat cultivars (XM55 and YM158) with PacBio sequencing. Of these, 5,309 long-non-coding RNAs, 1,574 fusion genes and 739 transcription factors were identified with the FLNC sequences. These full-length transcripts were derived from 49,368 genes, including 47.28% of the genes annotated in IWGSC RefSeq v1.0 and 40.86% genes encoded two or more isoforms, while 27.31% genes in the genome annotation of IWGSC RefSeq v1.0 were multiple-exon genes encoding two or more isoforms. Meanwhile, the individuals with waterlogging treatments (WL) and control group (CK) were selected for Illumina sequencing. Totally, 6,829 differentially expressed genes (DEGs) were detected from four pairwise comparisons. Notably, 942 DEGs were overlapped in the two comparisons (i.e., XM55-WL vs. YM158-WL and YM158-WL vs. YM158-CK). Undoubtedly, the genes involved in photosynthesis were downregulated after waterlogging treatment in two cultivars. Notably, the genes related to steroid biosynthesis, steroid hormone biosynthesis, and downstream plant hormone signal transduction were significantly upregulated after the waterlogging treatment, and the YM158 variety revealed different genetic regulation patterns compared with XM55. The findings provided valuable insights into unveiling regulation mechanisms of waterlogging stress in wheat at anthesis and contributed to molecular selective breeding of new wheat cultivars in future.

Comprehensive molecular dissection of TIFY Transcription factors reveal their dynamic responses to biotic and abiotic stress in wheat (Triticum aestivum L.)

The plant specific TIFY (previously known as ZIM) transcription factor (TF) family plays crucial roles in cross talk between Jasmonic Acid and other phytohormones like gibberellins, salicylic acid, abscisic acid, auxin, and ethylene signaling pathways. Wheat yield is severely affected by rust diseases and many abiotic stresses, where different phytohormone signaling pathways are involved. TIFYs have been studied in many plants yet reports describing their molecular structure and function in wheat are lacking. In the present study, we have identified 23 novel TIFY genes in wheat genome using in silico approaches. The identified proteins were characterized based on their conserved domains and phylogenetically classified into nine subfamilies. Chromosomal localization of the identified TIFY genes showed arbitrary distribution. Forty cis-acting elements including phytohormone, stress and light receptive elements were detected in the upstream regions of TIFY genes. Seventeen wheat microRNAs targeted the identified wheat TIFY genes. Gene ontological studies revealed their major contribution in defense response and phytohormone signaling. Secondary structure of TIFY proteins displayed the characteristic alpha-alpha-beta fold. Synteny analyses indicated all wheat TIFY genes had orthologous sequences in sorghum, rice, maize, barley and Brachypodium indicating presence of similar TIFY domains in monocot plants. Six TIFY genes had been cloned from wheat genomic and cDNA. Sequence characterization revealed similar characteristics as the in silico identified novel TIFY genes. Tertiary structures predicted the active sites in these proteins to play critical roles in DNA binding. Expression profiling of TIFY genes showed their contribution during incompatible and compatible leaf rust infestation. TIFY genes were also highly expressed during the initial hours of phytohormone induced stress. This study furnishes fundamental information on characterization and putative functions of TIFY genes in wheat.