Structural evolution of wheat chromosomes 4A, 5A, and 7B and its impact on recombination

Ancestral genome reconstruction and the evolution of chromosomal rearrangements in Triticeae

Chromosomal rearrangements (CRs) often cause phenotypic variations. Although several major rearrangements have been identified in Triticeae, a comprehensive study of the order, timing, and breakpoints of CRs has not been conducted. Here, we reconstruct high-quality ancestral genomes for the most recent common ancestor (MRCA) of the Triticeae, and the MRCA of the wheat lineage (Triticum and Aegilops). The protogenes of MRCA of the Triticeae and the wheat lineage are 22,894 and 29,060, respectively, which were arranged in their ancestral order. By partitioning modern Triticeae chromosomes into sets of syntenic regions and linking each to the corresponding protochromosomes, we revisit the rye chromosome structural evolution and propose alternative evolutionary routes. The previously identified 4L/5L reciprocal translocation in rye and Triticum urartu are found to have occurred independently and are unlikely the result of chromosomal introgression following distant hybridization. We also clarify that the 4AL/7BS translocation in tetraploid wheat was a bidirectional rather than unidirectional translocation event. Lastly, we identify several breakpoints in protochromosomes that independently reoccur following Triticeae evolution, representing potential CR hotspots. This study demonstrates that these reconstructed ancestral genomes can serve as special comparative references and facilitate a better understanding of the evolution of structural rearrangements in Triticeae.

A highly potential Zn biofortification tool: MTP1 in Triticum aestivum

TaMTPs belong to metal tolerance proteins (MTPs) family in common wheat and have significant potential to address the "hidden hunger" caused by inadequate dietary intake of a key micronutrient (Zn). In this study, a total of 33 MTP members in Triticum aestivum were identified, among which six TaMTP1-likes were closely related to Arabidopsis thaliana MTP1 and were designated as TaMTP1-A/B/D and TaMTP1.1-A/B/D. When heterologously expressed in yeast mutants, TaMTP1-likes complemented their hypersensitivity to Zn and Co, and three of the most metal-resistant members, TaMTP1-A, TaMTP1-D and TaMTP1.1-B, were selected for further subcellular localization and functional experiment in Arabidopsis and rice. The results showed that all three proteins were localized in the vacuole membrane, that TaMTP1-D was more resistant to Zn and less resistant to Co than other TaMTP1-like members, and that TaMTP1-D was expressed at a higher level in the endosperm than other members. All results reveal that the use of TaMTP1-D for biofortification can substantially increase the content of Zn in the edible part of wheat and avoid the overaccumulation of Co, suggesting that TaMTP1-D is a potential Zn biofortifier / bioreinforcement.

Dissecting the genetic basis of resistance to Soil-borne cereal mosaic virus (SBCMV) in durum wheat by bi-parental mapping and GWAS

Soil-borne cereal mosaic virus (SBCMV), the causative agent of wheat mosaic, is a Furovirus challenging wheat production all over Europe. Differently from bread wheat, durum wheat shows greater susceptibility and stronger yield penalties, so identification and genetic characterization of resistance sources are major targets for durum genetics and breeding. The Sbm1 locus providing high level of resistance to SBCMV was mapped in bread wheat to the 5DL chromosome arm (Bass in Genome 49:1140-1148, 2006). This excluded the direct use of Sbm1 for durum wheat improvement. Only one major QTL has been mapped in durum wheat, namely QSbm.ubo-2B, on the 2BS chromosome region coincident with Sbm2, already known in bread wheat as reported (Bayles in HGCA Project Report, 2007). Therefore, QSbm.ubo-2B = Sbm2 is considered a pillar for growing durum in SBCMV-affected areas. Herein, we report the fine mapping of Sbm2 based on bi-parental mapping and GWAS, using the Infinium 90 K SNP array and high-throughput KASP®. Fine mapping pointed out a critical haploblock of 3.2 Mb defined by concatenated SNPs successfully converted to high-throughput KASP® markers coded as KUBO. The combination of KUBO-27, wPt-2106-ASO/HRM, KUBO-29, and KUBO-1 allows unequivocal tracing of the Sbm2-resistant haplotype. The interval harbors 52 high- and 41 low-confidence genes, encoding 17 cytochrome p450, three receptor kinases, two defensins, and three NBS-LRR genes. These results pave the way for Sbm2 positional cloning. Importantly, the development of Sbm2 haplotype tagging KASP® provides a valuable case study for improving efficacy of the European variety testing system and, ultimately, the decision-making process related to varietal characterization and choice.

Genome wide characterization and expression analysis of CrRLK1L gene family in wheat unravels their roles in development and stress-specific responses

Catharanthus roseus receptor-like kinase 1-like (CrRLK1L) genes encode a subfamily of receptor-like kinases (RLK) that regulate diverse processes during plant growth, development, and stress responses. The first CrRLK1L was identified from the Catharanthus roseus, commonly known as Madagascar periwinkle. Subsequently, CrRLK1L gene families have been characterized in many plants. The genome of T. aestivum encodes 15 CrRLK1L genes with 43 paralogous copies, with three homeologs each, except for -2-D and -7-A, which are absent. Chromosomal localization analysis revealed a markedly uneven distribution of CrRLK1L genes across seven different chromosomes, with chromosome 4 housing the highest number of genes, while chromosome 6 lacked any CrRLK1L genes. Tissue-specific gene expression analysis revealed distinct expression patterns among the gene family members, with certain members exhibiting increased expression in reproductive tissues. Gene expression analysis in response to various abiotic and biotic stress conditions unveiled differential regulation of gene family members. Cold stress induces CrRLK1Ls -4-B and -15-A while downregulating -3-A and -7B. Drought stress upregulates -9D, contrasting with the downregulation of -7D. CrRLK1L-15-B and -15-D were highly induced in response to 1 hr of heat, and combined drought and heat stress, whereas -10-B is downregulated. Similarly, in response to NaCl stress, only CrRLK1L1 homeologs were induced. Fusarium graminearum and Claviceps purpurea inoculation induces homeologs of CrRLK1L-6 and -7. The analysis of cis-acting elements in the promoter regions identified elements crucial for plant growth and developmental processes. This comprehensive genome-wide analysis and expression study provides valuable insights into the essential functions of CrRLK1L members in wheat.

A chromosome level genome assembly of Pseudoroegneria Libanotica reveals a key Kcs gene involves in the cuticular wax elongation for drought resistance

BACKGROUND

The genus Pseudoroegneria (Nevski) Löve (Triticeae, Poaceae), whose genome symbol was designed as "St", accounts for more than 60% of perennial Triticeae species. The diploid species Psudoroegneria libanotica (2n = 14) contains the most ancient St genome, exhibited strong drought resistance, and was morphologically covered by cuticular wax on the aerial part. Therefore, the St-genome sequencing data could provide fundamental information for studies of genome evolution and reveal its mechanisms of cuticular wax and drought resistance.

RESULTS

In this study, we reported the chromosome-level genome assembly for the St genome of Pse. libanotica, with a total size of 2.99 Gb. 46,369 protein-coding genes annotated and 71.62% was repeat sequences. Comparative analyses revealed that the genus Pseudoroegneria diverged during the middle and late Miocene. During this period, unique genes, gene family expansion, and contraction in Pse. libanotica were enriched in biotic and abiotic stresses, such as fatty acid biosynthesis which may greatly contribute to its drought adaption. Furthermore, we investigated genes associated with the cuticular wax formation and water deficit and found a new Kcs gene evm.TU.CTG175.54. It plays a critical role in the very long chain fatty acid (VLCFA) elongation from C18 to C26 in Pse. libanotica. The function needs more evidence to be verified.

CONCLUSIONS

We sequenced and assembled the St genome in Triticeae and discovered a new KCS gene that plays a role in wax extension to cope with drought. Our study lays a foundation for the genome diversification of Triticeae species and deciphers cuticular wax formation genes involved in drought resistance.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

Chromosome translocation affects multiple phenotypes, causes genome-wide dysregulation of gene expression, and remodels metabolome in hexaploid wheat

Chromosomal rearrangements (CRs) may occur in newly formed polyploids due to compromised meiotic fidelity. Moreover, CRs can be more readily tolerated in polyploids allowing their longer-term retention and hence potential spreading/fixation within a lineage. The direct functional consequences of CRs in plant polyploids remain unexplored. Here, we identified a heterozygous individual from a synthetic allohexaploid wheat in which the terminal parts of the long-arms of chromosomes 2D (approximately 193 Mb) and 4A (approximately 167 Mb) were reciprocally translocated. Five homogeneous translocation lines including both unbalanced and balanced types were developed by selfing fertilization of the founder mutant (RT [2DL; 4AL]-ter/1, reciprocal translocation). We investigated impacts of these translocations on phenotype, genome-wide gene expression and metabolome. We find that, compared with sibling wild-type, CRs in the form of both unbalanced and balanced translocations induced substantial changes of gene expression primarily via trans-regulation in the nascent allopolyploid wheat. The CRs also manifested clear phenotypic and metabolic consequences. In particular, the genetically balanced, stable reciprocal translocations lines showed immediate enhanced reproductive fitness relative to wild type. Our results underscore the profound impact of CRs on gene expression in nascent allopolyploids with wide-ranging phenotypic and metabolic consequences, suggesting CRs are an important source of genetic variation that can be exploited for crop breeding.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Historical Selection, Adaptation Signatures, and Ambiguity of Introgressions in Wheat

Wheat was one of the crops domesticated in the Fertile Crescent region approximately 10,000 years ago. Despite undergoing recent polyploidization, hull-to-free-thresh transition events, and domestication bottlenecks, wheat is now grown in over 130 countries and accounts for a quarter of the world's cereal production. The main reason for its widespread success is its broad genetic diversity that allows it to thrive in different environments. To trace historical selection and hybridization signatures, genome scans were performed on two datasets: approximately 113K SNPs from 921 predominantly bread wheat accessions and approximately 110K SNPs from about 400 wheat accessions representing all ploidy levels. To identify environmental factors associated with the loci, a genome-environment association (GEA) was also performed. The genome scans on both datasets identified a highly differentiated region on chromosome 4A where accessions in the first dataset were dichotomized into a group (n = 691), comprising nearly all cultivars, wild emmer, and most landraces, and a second group (n = 230), dominated by landraces and spelt accessions. The grouping of cultivars is likely linked to their potential ancestor, bread wheat cv. Norin-10. The 4A region harbored important genes involved in adaptations to environmental conditions. The GEA detected loci associated with latitude and temperature. The genetic signatures detected in this study provide insight into the historical selection and hybridization events in the wheat genome that shaped its current genetic structure and facilitated its success in a wide spectrum of environmental conditions. The genome scans and GEA approaches applied in this study can help in screening the germplasm housed in gene banks for breeding, and for conservation purposes.

Single-gene FISH maps and major chromosomal rearrangements in Elymus sibiricus and E. nutans

BACKGROUND

Chromosomal variations have been revealed in both E. sibiricus and E. nutans, but chromosomal structural variations, such as intra-genome translocations and inversions, are still not recognized due to the cytological limitations of previous studies. Furthermore, the syntenic relationship between both species and wheat chromosomes remains unknown.

RESULTS

Fifty-nine single-gene fluorescence in situ hybridization (FISH) probes, including 22 single-gene probes previously mapped on wheat chromosomes and other newly developed probes from the cDNA of Elymus species, were used to characterize the chromosome homoeologous relationship and collinearity of both E. sibiricus and E. nutans with those of wheat. Eight species-specific chromosomal rearrangements (CRs) were exclusively identified in E. sibiricus, including five pericentric inversions in 1H, 2H, 3H, 6H and 2St; one possible pericentric inversion in 5St; one paracentric inversion in 4St; and one reciprocal 4H/6H translocation. Five species-specific CRs were identified in E. nutans, including one possible pericentric inversion in 2Y, three possible pericentric multiple-inversions in 1H, 2H and 4Y, and one reciprocal 4Y/5Y translocation. Polymorphic CRs were detected in three of the six materials in E. sibiricus, which were mainly represented by inter-genomic translocations. More polymorphic CRs were identified in E. nutans, including duplication and insertion, deletion, pericentric inversion, paracentric inversion, and intra- or inter-genomic translocation in different chromosomes.

CONCLUSIONS

The study first identified the cross-species homoeology and the syntenic relationship between E. sibiricus, E. nutans and wheat chromosomes. There are distinct different species-specific CRs between E. sibiricus and E. nutans, which may be due to their different polyploidy processes. The frequencies of intra-species polymorphic CRs in E. nutans were higher than that in E. sibiricus. To conclude, the results provide new insights into genome structure and evolution and will facilitate the utilization of germplasm diversity in both E. sibiricus and E. nutans.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in wheat (Triticum aestivum L.)

Calmodulin (CaM) and calmodulin-like (CML) genes are widely involved in plant growth and development and mediating plant stress tolerance. However, the whole genome scale studies about CaM and CML gene families have not been done in wheat, and the possible functions of most wheat CaM/CML gene members are still unknown. In this study, a total of 18 TaCaM and 230 TaCML gene members were identified in wheat genome. Among these genes, 28 TaCaM/CML gene members have 74 duplicated copies, while 21 genes have 48 transcript variants, resulting in 321 putative TaCaM/CML transcripts totally. Phylogenetic tree analysis showed that they can be classified into 7 subfamilies. Similar gene structures and protein domains can be found in members of the same gene cluster. The TaCaM/CML genes were spread among all 21 chromosomes with unbalanced distributions, while most of the gene clusters contained 3 homoeologous genes located in the same homoeologous chromosome group. Synteny analysis showed that most of TaCaM/CMLs gene members can be found with 1-4 paralogous genes in T. turgidum and Ae. Tauschii. High numbers of cis-acting elements related to plant hormones and stress responses can be observed in the promoters of TaCaM/CMLs. The spatiotemporal expression patterns showed that most of the TaCaM/TaCML genes can be detected in at least one tissue. The expression levels of TaCML17, 21, 30, 50, 59 and 75 in the root or shoot can be up-regulated by abiotic stresses, suggesting that TaCML17, 21, 30, 50, 59 and 75 may be related with responses to abiotic stresses in wheat. The spatiotemporal expression patterns of TaCaM/CML genes indicated they may be involved widely in wheat growth and development. Our results provide important clues for exploring functions of TaCaMs/CMLs in growth and development as well as responses to abiotic stresses in wheat in the future.

Pseudorogneria libanotica Intraspecific Genetic Polymorphism Revealed by Fluorescence In Situ Hybridization with Newly Identified Tandem Repeats and Wheat Single-Copy Gene Probes

The genus Pseudoroegneria (Nevski) Löve (Triticeae, Poaceae) with its genome abbreviated 'St' accounts for more than 60% of perennial Triticeae species. The diploid species Psudoroegneria libanotica (2n = 14) contains the most ancient St genome. Therefore, investigating its chromosomes could provide some fundamental information required for subsequent studies of St genome evolution. Here, 24 wheat cDNA probes covering seven chromosome groups were mapped in P. libanotica to distinguish homoelogous chromosomes, and newly identified tandem repeats were performed to differentiate seven chromosome pairs. Using these probes, we investigated intraspecific population chromosomal polymorphism of P. libanotica. We found that (i) a duplicated fragment of the 5St long arm was inserted into the short arm of 2St; (ii) asymmetrical fluorescence in situ hybridization (FISH) hybridization signals among 2St, 5St, and 7St homologous chromosome pairs; and (iii) intraspecific population of polymorphism in P. libanotica. These observations established the integrated molecular karyotype of P. libanotica. Moreover, we suggested heterozygosity due to outcrossing habit and adaptation to the local climate of P. libanotica. Specifically, the generated STlib_96 and STlib_98 repeats showed no cross-hybridization signals with wheat chromosomes, suggesting that they are valuable for identifying alien chromosomes or introgressed fragments of wild relatives in wheat.

Comparative genomic and transcriptomic analyses uncover the molecular basis of high nitrogen-use efficiency in the wheat cultivar Kenong 9204

Studying the regulatory mechanisms that drive nitrogen-use efficiency (NUE) in crops is important for sustainable agriculture and environmental protection. In this study, we generated a high-quality genome assembly for the high-NUE wheat cultivar Kenong 9204 and systematically analyzed genes related to nitrogen uptake and metabolism. By comparative analyses, we found that the high-affinity nitrate transporter gene family had expanded in Triticeae. Further studies showed that subsequent functional differentiation endowed the expanded family members with saline inducibility, providing a genetic basis for improving the adaptability of wheat to nitrogen deficiency in various habitats. To explore the genetic and molecular mechanisms of high NUE, we compared genomic and transcriptomic data from the high-NUE cultivar Kenong 9204 (KN9204) and the low-NUE cultivar Jing 411 and quantified their nitrogen accumulation under high- and low-nitrogen conditions. Compared with Jing 411, KN9204 absorbed significantly more nitrogen at the reproductive stage after shooting and accumulated it in the shoots and seeds. Transcriptome data analysis revealed that nitrogen deficiency clearly suppressed the expression of genes related to cell division in the young spike of Jing 411, whereas this suppression of gene expression was much lower in KN9204. In addition, KN9204 maintained relatively high expression of NPF genes for a longer time than Jing 411 during seed maturity. Physiological and transcriptome data revealed that KN9204 was more tolerant of nitrogen deficiency than Jing 411, especially at the reproductive stage. The high NUE of KN9204 is an integrated effect controlled at different levels. Taken together, our data provide new insights into the molecular mechanisms of NUE and important gene resources for improving wheat cultivars with a higher NUE trait.

Chromosome-scale assembly and annotation of the perennial ryegrass genome

BACKGROUND

The availability of chromosome-scale genome assemblies is fundamentally important to advance genetics and breeding in crops, as well as for evolutionary and comparative genomics. The improvement of long-read sequencing technologies and the advent of optical mapping and chromosome conformation capture technologies in the last few years, significantly promoted the development of chromosome-scale genome assemblies of model plants and crop species. In grasses, chromosome-scale genome assemblies recently became available for cultivated and wild species of the Triticeae subfamily. Development of state-of-the-art genomic resources in species of the Poeae subfamily, which includes important crops like fescues and ryegrasses, is lagging behind the progress in the cereal species.

RESULTS

Here, we report a new chromosome-scale genome sequence assembly for perennial ryegrass, obtained by combining PacBio long-read sequencing, Illumina short-read polishing, BioNano optical mapping and Hi-C scaffolding. More than 90% of the total genome size of perennial ryegrass (approximately 2.55 Gb) is covered by seven pseudo-chromosomes that show high levels of collinearity to the orthologous chromosomes of Triticeae species. The transposon fraction of perennial ryegrass was found to be relatively low, approximately 35% of the total genome content, which is less than half of the genome repeat content of cultivated cereal species. We predicted 54,629 high-confidence gene models, 10,287 long non-coding RNAs and a total of 8,393 short non-coding RNAs in the perennial ryegrass genome.

CONCLUSIONS

The new reference genome sequence and annotation presented here are valuable resources for comparative genomic studies in grasses, as well as for breeding applications and will expedite the development of productive varieties in perennial ryegrass and related species.

Introgression of the Triticum timopheevii Genome Into Wheat Detected by Chromosome-Specific Kompetitive Allele Specific PCR Markers

Triticum timopheevii (2n = 28, A t A t GG) is a tetraploid wild relative species with great potential to increase the genetic diversity of hexaploid wheat Triticum aestivum (2n = 42, AABBDD) for various important agronomic traits. A breeding scheme that propagated advanced backcrossed populations of wheat-T. timopheevii introgression lines through further backcrossing and self-fertilisation resulted in the generation of 99 introgression lines (ILs) that carried 309 homozygous segments from the A t and G subgenomes of T. timopheevii. These introgressions contained 89 and 74 unique segments from the A t and G subgenomes, respectively. These overlapping segments covered 98.9% of the T. timopheevii genome that has now been introgressed into bread wheat cv. Paragon including the entirety of all T. timopheevii chromosomes via varying sized segments except for chromosomes 3A t , 4G, and 6G. Homozygous ILs contained between one and eight of these introgressions with an average of three per introgression line. These homozygous introgressions were detected through the development of a set of 480 chromosome-specific Kompetitive allele specific PCR (KASP) markers that are well-distributed across the wheat genome. Of these, 149 were developed in this study based on single nucleotide polymorphisms (SNPs) discovered through whole genome sequencing of T. timopheevii. A majority of these KASP markers were also found to be T. timopheevii subgenome specific with 182 detecting A t subgenome and 275 detecting G subgenome segments. These markers showed that 98% of the A t segments had recombined with the A genome of wheat and 74% of the G genome segments had recombined with the B genome of wheat with the rest recombining with the D genome of wheat. These results were validated through multi-colour in situ hybridisation analysis. Together these homozygous wheat-T. timopheevii ILs and chromosome-specific KASP markers provide an invaluable resource to wheat breeders for trait discovery to combat biotic and abiotic stress factors affecting wheat production due to climate change.

Dissection of Structural Reorganization of Wheat 5B Chromosome Associated With Interspecies Recombination Suppression

Chromosomal rearrangements that lead to recombination suppression can have a significant impact on speciation, and they are also important for breeding. The regions of recombination suppression in wheat chromosome 5B were identified based on comparisons of the 5B map of a cross between the Chinese Spring (CS) variety of hexaploid wheat and CS-5Bdic (genotype CS with 5B substituted with its homologue from tetraploid Triticum dicoccoides) with several 5B maps of tetraploid and hexaploid wheat. In total, two regions were selected in which recombination suppression occurred in cross CS × CS-5Bdic when compared with other maps: one on the short arm, 5BS_RS, limited by markers BS00009810/BS00022336, and the second on the long arm, 5BL_RS, between markers Ra_c10633_2155 and BS00087043. The regions marked as 5BS_RS and 5BL_RS, with lengths of 5 Mb and 3.6 Mb, respectively, were mined from the 5B pseudomolecule of CS and compared to the homoeologous regions (7.6 and 3.8 Mb, respectively) of the 5B pseudomolecule of Zavitan (T. dicoccoides). It was shown that, in the case of 5BS_RS, the local heterochromatin islands determined by the satellite DNA (119.2) and transposable element arrays, as well as the dissimilarity caused by large insertions/deletions (chromosome rearrangements) between 5BSs aestivum/dicoccoides, are likely the key determinants of recombination suppression in the region. Two major and two minor segments with significant loss of similarity were recognized within the 5BL_RS region. It was shown that the loss of similarity, which can lead to suppression of recombination in the 5BL_RS region, is caused by chromosomal rearrangements, driven by the activity of mobile genetic elements (both DNA transposons and long terminal repeat retrotransposons) and their divergence during evolution. It was noted that the regions marked as 5BS_RS and 5BL_RS are associated with chromosomal rearrangements identified earlier by С-banding analysis of intraspecific polymorphism of tetraploid emmer wheat. The revealed divergence in 5BS_RS and 5BL_RS may be a consequence of interspecific hybridization, plant genetic adaptation, or both.

Characterization of the heterotrimeric G protein gene families in Triticum aestivum and related species

This study characterizes the heterotrimeric G protein gene families in Triticum aestivum, their tissue-specific expression patterns during development and in response to biotic and abiotic stress conditions. There are three Gα genes, three Gβ and 12 Gγ genes, totaling 18 genes encoding heterotrimeric G proteins in the hexaploid wheat genome. Each haploid genome of the hexaploid T. aestivum has a single gene encoding the α subunit of the heterotrimeric G protein complex, GA1, a single Gβ and four Gγ genes. Each gene has three homeologous copies in the A, B and D genomes. The physical interaction between the Gβ (Gpb) and two Gγ subunits, Gpg1 and Gpg2, was shown through bimolecular fluorescence complementation (BiFC). The gene expression in response to biotic and abiotic stresses showed both up-regulation and down-regulation of members of the gene families. Gγ2-B and Gγ2-D are significantly upregulated during heat stress, GA1-D is upregulated by cold stress and Gγ1-A and Gγ1-D were upregulated by Fusarium graminearum inoculation in a F. graminearum resistant cultivar. This suggests that these members may play roles in biotic and abiotic signaling pathways and the roles of these genes within these pathways need further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03156-9.

Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines

Many wild-relative species are being used in prebreeding programs to increase the genetic diversity of wheat (Triticum aestivum L.). Genotyping tools such as single nucleotide polymorphism (SNP)-based arrays and molecular markers have been widely used to characterize wheat-wild relative introgression lines. However, due to the polyploid nature of the recipient wheat genome, it is difficult to develop SNP-based Kompetitive allele-specific polymerase chain reaction (KASP) markers that are codominant to track the introgressions from the wild species. Previous attempts to develop KASP markers have involved both exome- and polymerase chain reaction (PCR)-amplicon-based sequencing of the wild species. But chromosome-specific KASP assays have been hindered by homoeologous SNPs within the wheat genome. This study involved whole genome sequencing of the diploid wheat wild relative Amblyopyrum muticum (Boiss.) Eig and development of a de novo SNP discovery pipeline that generated ∼38,000 SNPs in unique wheat genome sequences. New assays were designed to increase the density of Am. muticum polymorphic KASP markers. With a goal of one marker per 60 Mbp, 335 new KASP assays were validated as diagnostic for Am. muticum in a wheat background. Together with assays validated in previous studies, 498 well distributed chromosome-specific markers were used to recharacterize previously genotyped wheat-Am. muticum doubled haploid (DH) introgression lines. The chromosome-specific nature of the KASP markers allowed clarification of which wheat chromosomes were involved with recombination events or substituted with Am. muticum chromosomes and the higher density of markers allowed detection of new small introgressions in these DH lines.

Molecular Cytogenetic Characterization of Novel 1RS.1BL Translocation and Complex Chromosome Translocation Lines with Stripe Rust Resistance

Rye is the most important source for the genetic improvement of wheat. In this study, two stable wheat-rye primary 1RS.1BL translocation lines, RT855-13 and RT855-14, were selected and identified by acid polyacrylamide gel electrophoresis (A-PAGE), co-dominant PCR, and multi-color fluorescence in situ hybridization (MC-FISH) from the progeny of the crossing of the wheat cultivar Mianyang11 and a Chinese rye Weining. When more than two independent, simple reciprocal translocations are involved in a carrier, they are defined as complex chromosome translocations (CCT). The MC-FISH results also indicated that CCT occurred in RT855-13; namely that, besides 1RS.1BL translocation chromosomes, there are other two pairs of balanced reciprocal translocations. It was demonstrated that the interchange between a distal segment of 4B and long arm of 3D occurred in the RT855-13. The novel translocation chromosomes in wheat were recorded as 3DS.4BS^DS and 3DL-4BS^PS.4BL. Reports about CCT as a genetic resource in plant breeding programs are scarce. Both lines expressed high resistance to Puccinia striiformis f. sp. tritici, which are prevalent in China and are virulent on Yr9, and the CCT line RT855-13 retained better resistance as adult plants compared with RT855-14 in the field. Both lines, especially the CCT line RT855-13, exhibited better agronomic traits than their wheat parent, Mianyang11, indicating that both translocation lines could potentially be used for wheat improvement. The results also indicated that the position effects of CCT can lead to beneficial variations in agronomic and resistant traits, making them a valuable genetic resource to wheat breeding programs.

Wheat heat tolerance is impaired by heightened deletions in the distal end of 4AL chromosomal arm

Heat stress (HS) causes substantial damages to worldwide crop production. As a cool season crop, wheat (Triticum aestivum) is sensitive to HS-induced damages. To support the genetic improvement of wheat HS tolerance (HST), we conducted fine mapping of TaHST1, a locus required for maintaining wheat vegetative and reproductive growth under elevated temperatures. TaHST1 was mapped to the distal terminus of 4AL chromosome arm using genetic populations derived from two BC₆ F₆ breeding lines showing tolerance (E6015-4T) or sensitivity (E6015-3S) to HS. The 4AL region carrying TaHST1 locus was approximately 0.949 Mbp and contained the last 19 high confidence genes of 4AL according to wheat reference genome sequence. Resequencing of E6015-3S and E6015-4T and haplotype analysis of 3087 worldwide wheat accessions revealed heightened deletion polymorphisms in the distal 0.949 Mbp region of 4AL, which was confirmed by the finding of frequent gene losses in this region in eight genome-sequenced hexaploid wheat cultivars. The great majority (86.36%) of the 3087 lines displayed different degrees of nucleotide sequence deletions, with only 13.64% of them resembling E6015-4T in this region. These deletions can impair the presence and/or function of TaHST1 and surrounding genes, thus rendering global wheat germplasm vulnerable to HS or other environmental adversities. Therefore, conscientious and urgent efforts are needed in global wheat breeding programmes to optimize the structure and function of 4AL distal terminus by ensuring the presence of TaHST1 and surrounding genes. The new information reported here will help to accelerate the ongoing global efforts in improving wheat HST.

Generation of Doubled Haploid Wheat-Triticum urartu Introgression Lines and Their Characterisation Using Chromosome-Specific KASP Markers

Wheat is one of the most important food and protein sources in the world and although, in recent years wheat breeders have achieved yield gains, they are not sufficient to meet the demands of an ever-growing population. Development of high yielding wheat varieties, resilient to abiotic and biotic stress resulting from climate change, has been limited by wheat's narrow genetic base. In contrast to wheat, the wild relatives of wheat provide a vast reservoir of genetic variation for most, if not all, agronomic traits. Previous studies by the authors have shown the transfer of genetic variation from T. urartu into bread wheat. However, before the introgression lines can be exploited for trait analysis, they are required to have stable transmission of the introgressions to the next generation. In this work, we describe the generation of 86 doubled haploid (DH) wheat-T. urartu introgression lines that carry homozygous introgressions which are stably inherited. The DH lines were characterised using the Axiom^® Wheat Relative Genotyping Array and 151 KASP markers to identify 65 unique T. urartu introgressions in a bread wheat background. DH production has helped accelerate the breeding process and facilitated the early release of homozygous wheat-T. urartu introgression lines. Together with the KASP markers, this valuable resource could greatly advance identification of beneficial alleles that can be used in wheat improvement.

A Collinearity-Incorporating Homology Inference Strategy for Connecting Emerging Assemblies in the Triticeae Tribe as a Pilot Practice in the Plant Pangenomic Era

Plant genome sequencing has dramatically increased, and some species even have multiple high-quality reference versions. Demands for clade-specific homology inference and analysis have increased in the pangenomic era. Here we present a novel method, GeneTribe (https://chenym1.github.io/genetribe/), for homology inference among genetically similar genomes that incorporates gene collinearity and shows better performance than traditional sequence-similarity-based methods in terms of accuracy and scalability. The Triticeae tribe is a typical allopolyploid-rich clade with complex species relationships that includes many important crops, such as wheat, barley, and rye. We built Triticeae-GeneTribe (http://wheat.cau.edu.cn/TGT/), a homology database, by integrating 12 Triticeae genomes and 3 outgroup model genomes and implemented versatile analysis and visualization functions. With macrocollinearity analysis, we were able to construct a refined model illustrating the structural rearrangements of the 4A-5A-7B chromosomes in wheat as two major translocation events. With collinearity analysis at both the macro- and microscale, we illustrated the complex evolutionary history of homologs of the wheat vernalization gene Vrn2, which evolved as a combined result of genome translocation, duplication, and polyploidization and gene loss events. Our work provides a useful practice for connecting emerging genome assemblies, with awareness of the extensive polyploidy in plants, and will help researchers efficiently exploit genome sequence resources.

Dissection and physical mapping of wheat chromosome 7B by inducing meiotic recombination with its homoeologues in Aegilops speltoides and Thinopyrum elongatum

We constructed a homoeologous recombination-based bin map of wheat chromosome 7B, providing a unique physical framework for further study of chromosome 7B and its homoeologues in wheat and its relatives. Homoeologous recombination leads to the dissection and diversification of the wheat genome. Advances in genome sequencing and genotyping have dramatically improved the efficacy and throughput of homoeologous recombination-based genome studies and alien introgression in wheat and its relatives. In this study, we aimed to physically dissect and map wheat chromosome 7B by inducing meiotic recombination of chromosome 7B with its homoeologues 7E in Thinopyrum elongatum and 7S in Aegilops speltoides. The special genotypes, which were double monosomic for chromosomes 7B' + 7E' or 7B' + 7S' and homozygous for the ph1b mutant, were produced to enhance 7B - 7E and 7B - 7S recombination. Chromosome-specific DNA markers were developed and used to pre-screen the large recombination populations for 7B - 7E and 7B - 7S recombinants. The DNA marker-mediated preselections were verified by fluorescent genomic in situ hybridization (GISH). In total, 29 7B - 7E and 61 7B - 7S recombinants and multiple chromosome aberrations were recovered and delineated by GISH and the wheat 90 K SNP assay. Integrated GISH and SNP analysis of the recombinants physically mapped the recombination breakpoints and partitioned wheat chromosome 7B into 44 bins with 523 SNPs assigned within. A composite bin map was constructed for chromosome 7B, showing the bin size and physical distribution of SNPs. This provides a unique physical framework for further study of chromosome 7B and its homoeologues. In addition, the 7B - 7E and 7B - 7S recombinants extend the genetic variability of wheat chromosome 7B and represent useful germplasm for wheat breeding. Thereby, this genomics-enabled chromosome engineering approach facilitates wheat genome study and enriches the gene pool of wheat improvement.

Chromosome genomics uncovers plant genome organization and function

The identification of causal genomic loci and their interactions underlying various traits in plants has been greatly aided by progress in understanding the organization of the nuclear genome. This provides clues to the responses of plants to environmental stimuli at the molecular level. Apart from other uses, these insights are needed to fully explore the potential of new breeding techniques that rely on genome editing. However, genome analysis and sequencing is not straightforward in the many agricultural crops and their wild relatives that possess large and complex genomes. Chromosome genomics streamlines this task by dissecting the genome to single chromosomes whose DNA is then used instead of nuclear DNA. This results in a massive and lossless reduction in DNA sample complexity, reduces the time and cost of the experiment, and simplifies data interpretation. Flow cytometric sorting of condensed mitotic chromosomes makes it possible to purify single chromosomes in large quantities, and as the DNA remains intact this process can be coupled successfully with many techniques in molecular biology and genomics. Since the first experiments with flow cytometric sorting in the late 1980s, numerous applications have been developed, and chromosome genomics has been having a significant impact in many areas of research, including the sequencing of complex genomes of important crops and gene cloning. This review discusses these applications, describes their contribution to advancements in plant genome analysis and gene cloning, and outlines future directions.

A new mode of NPR1 action via an NB-ARC-NPR1 fusion protein negatively regulates the defence response in wheat to stem rust pathogen

NPR1 has been found to be a key transcriptional regulator in some plant defence responses. There are nine NPR1 homologues (TaNPR1) in wheat, but little research has been done to understand the function of those NPR1-like genes in the wheat defence response against stem rust (Puccinia graminis f. sp. tritici) pathogens. We used bioinformatics and reverse genetics approaches to study the expression and function of each TaNPR1. We found six members of TaNPR1 located on homoeologous group 3 chromosomes (designated as TaG3NPR1) and three on homoeologous group 7 chromosomes (designated as TaG7NPR1). The group 3 NPR1 proteins regulate transcription of SA-responsive PR genes. Downregulation of all the TaNPR1 homologues via virus-induced gene co-silencing resulted in enhanced resistance to stem rust. More specifically downregulating TaG7NPR1 homeologues or Ta7ANPR1 expression resulted in stem rust resistance phenotype. By contrast, knocking down TaG3NPR1 alone did not show visible phenotypic changes in response to the rust pathogen. Knocking out Ta7ANPR1 enhanced resistance to stem rust. The Ta7ANPR1 locus is alternatively spliced under pathogen inoculated conditions. We discovered a new mode of NPR1 action in wheat at the Ta7ANPR1 locus through an NB-ARC-NPR1 fusion protein negatively regulating the defence to stem rust infection.

Accurate calling of homeoallelic genotypes of iSelect markers using inbred structured populations

MOTIVATION

Genetic map construction is a foundational step in the analysis of structured experimental populations. For markers that hybridize to several genetically similar locations, or where several alleles are present (such as in multiparental populations), current methods often discard the marker or incorrectly call the genotypes. These errors result in information loss, or incorrect genotypes that can corrupt map construction.

RESULTS

We present a new approach for simultaneously performing genetic map construction and marker calling. Our new approach allows the calling of a larger number of markers, a larger number of unique alleles per marker and the correct use of markers which hybridize to multiple genetically similar locations. We demonstrate our new approach using simulations, a biparental wheat population and an eight-parent population of spring bread wheat. Applying our method to the eight-parent population increased the number of mapped markers by 71%. We show that the new genetic map allows the investigation of synteny in ways that were not previously possible in that dataset.

AVAILABILITY AND IMPLEMENTATION

The method described in this article has been incorporated into R package mpMap2. It is available from CRAN and also from https://github.com/rohan-shah/mpMap2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential chromatin accessibility landscape reveals structural and functional features of the allopolyploid wheat chromosomes

Background

Our understanding of how the complexity of the wheat genome influences the distribution of chromatin states along the homoeologous chromosomes is limited. Using a differential nuclease sensitivity assay, we investigate the chromatin states of the coding and repetitive regions of the allopolyploid wheat genome.

Results

Although open chromatin is found to be significantly enriched around genes, the majority of MNase-sensitive regions are located within transposable elements (TEs). Chromatin of the smaller D genome is more accessible than that of the larger A and B genomes. Chromatin states of different TEs vary among families and are influenced by the TEs’ chromosomal position and proximity to genes. While the chromatin accessibility of genes is influenced by proximity to TEs, and not by their position on the chromosomes, we observe a negative chromatin accessibility gradient along the telomere-centromere axis in the intergenic regions, positively correlated with the distance between genes. Both gene expression levels and homoeologous gene expression bias are correlated with chromatin accessibility in promoter regions. The differential nuclease sensitivity assay accurately predicts previously detected centromere locations. SNPs located within more accessible chromatin explain a higher proportion of genetic variance for a number of agronomic traits than SNPs located within more closed chromatin.

Conclusions

Chromatin states in the wheat genome are shaped by the interplay of repetitive and gene-encoding regions that are predictive of the functional and structural organization of chromosomes, providing a powerful framework for detecting genomic features involved in gene regulation and prioritizing genomic variation to explain phenotypes.

AutoCloner: automatic homologue-specific primer design for full-gene cloning in polyploids

Background

Polyploid organisms such as wheat complicate even the simplest of procedures in molecular biology. Whilst knowledge of genomic sequences in crops is increasing rapidly, the scientific community is still a long way from producing a full pan-genome for every species. Polymerase chain reaction and Sanger sequencing therefore remain widely used as methods for characterizing gene sequences in many varieties of crops. High sequence similarity between genomes in polyploids means that if primers are not homeologue-specific via the incorporation of a SNP at the 3’ tail, sequences other than the target sequence will also be amplified. Current consensus for gene cloning in wheat is to manually perform many steps in a long bioinformatics pipeline.

Results

Here we present AutoCloner (www.autocloner.com), a fully automated pipeline for crop gene cloning that includes a free-to-use web interface for users. AutoCloner takes a sequence of interest from the user and performs a basic local alignment search tool (BLAST) search against the genome assembly for their particular polyploid crop. Homologous sequences are then compiled with the input sequence into a multiple sequence alignment which is mined for single-nucleotide polymorphisms (SNPs). Various combinations of potential primers that cover the entire gene of interest are then created and evaluated by Primer3; the set of primers with the highest score, as well as all possible primers at every SNP location, are then returned to the user for polymerase chain reaction (PCR). We have successfully used AutoCloner to clone various genes of interest in the Apogee wheat variety, which has no current genome sequence. In addition, we have successfully run the pipeline on ~ 80,000 high-confidence gene models from a wheat genome assembly.

Conclusion

AutoCloner is the first tool to fully-automate primer design for gene cloning in polyploids, where previously the consensus within the wheat community was to perform this process manually. The web interface for AutoCloner provides a simple and effective polyploid primer-design method for gene cloning, with no need for researchers to download software or input any other details other than their sequence of interest.

Shaping Durum Wheat for the Future: Gene Expression Analyses and Metabolites Profiling Support the Contribution of BCAT Genes to Drought Stress Response

Global climate change, its implications for agriculture, and the complex scenario presented by the scientific community are of worldwide concern. Drought is a major abiotic stress that can restrict plants growth and yields, thus the identification of genotypes with higher adaptability to drought stress represents one of the primary goals in breeding programs. During abiotic stress, metabolic adaptation is crucial for stress tolerance, and accumulation of specific amino acids and/or as secondary metabolites deriving from amino acid metabolism may correlate with the increased tolerance to adverse environmental conditions. This work, focused on the metabolism of branched chain-amino acids (BCAAs) in durum wheat and the role of branched-chain amino acid aminotransferases (BCATs) in stress response. The role of BCATs in plant response to drought was previously proposed for Arabidopsis, where the levels of BCAAs were altered at the transcriptional level under drought conditions, triggering the onset of defense response metabolism. However, in wheat the role of BCAAs as a trigger of the onset of the drought defense response has not been elucidated. A comparative genomic approach elucidated the composition of the BCAT gene family in durum wheat. Here we demonstrate a tissue and developmental stage specificity of BCATs regulation in the drought response. Moreover, a metabolites profiling was performed on two contrasting durum wheat cultivars Colosseo and Cappelli resulting in the detection of a specific pattern of metabolites accumulated among genotypes and, in particular, in an enhanced BCAAs accumulation in the tolerant cv Cappelli further supporting a role of BCAAs in the drought defense response. The results support the use of gene expression and target metabolomic in modern breeding to shape new cultivars more resilient to a changing climate.

Exploiting the genome of Thinopyrum elongatum to expand the gene pool of hexaploid wheat

One hundred and thirty four introgressions from Thinopyrum elongatum have been transferred into a wheat background and were characterised using 263 SNP markers. Species within the genus Thinopyrum have been shown to carry genetic variation for a very wide range of traits including biotic and abiotic stresses and quality. Research has shown that one of the species within this genus, Th. elongatum, has a close relationship with the genomes of wheat making it a highly suitable candidate to expand the gene pool of wheat. Homoeologous recombination, in the absence of the Ph1 gene, has been exploited to transfer an estimated 134 introgressions from Th. elongatum into a hexaploid wheat background. The introgressions were detected and characterised using 263 single nucleotide polymorphism markers from a 35 K Axiom^® Wheat-Relative Genotyping Array, spread across seven linkage groups and validated using genomic in situ hybridisation. The genetic map had a total length of 187.8 cM and the average chromosome length was 26.8 cM. Comparative analyses of the genetic map of Th. elongatum and the physical map of hexaploid wheat confirmed previous work that indicated good synteny at the macro-level, although Th. elongatum does not contain the 4A/5A/7B translocation found in wheat.

Genome-Scale Analysis of Homologous Genes among Subgenomes of Bread Wheat (Triticum aestivum L.)

Determining the distribution and correspondence of genome-scale homologous genes in wheat are effective ways to uncover chromosome rearrangement that has occurred during crop evolution and domestication, which can contribute to improvements in crop breeding. High-resolution and comprehensive analysis of the wheat genome by the International Wheat Genome Sequencing Consortium (IWGSC) revealed a total of 88,733 high-confidence homologous genes of four major types (1:1:1, 1:1:0, 0:1:1 and 1:0:1) among the A, B and D subgenomes of wheat. This data was used to compare homologous gene densities among chromosomes, clarify their distribution and correspondence relationship, and compare their functional enrichment. The average density of 1:1:1 homologous genes was about 10 times more than the density of the other three types of homologous genes, although the homologous gene densities of the various chromosomes were similar within each homologous type. Three regions of exceptional density were detected in 1:1:1 homologous genes, the isolate peak on the tail of chromosome 4A, and the desert regions at the start of chromosome 7A and 7D. The correspondence between homologous genes of the wheat subgenomes demonstrated translocation between the tail segments of chromosome 4A and 5A, and the inversion of the segment of original 5A and 7B into the tail of 4A. The homologous genes on the inserting segments of 5A and 7B to 4A were highly enriched in nitrogen, primary metabolite and small molecular metabolism processes, compared with genes on other regions of the original 4A chromosome. This study provides a refined genome-scale reference of homologous genes for wheat molecular research and breeding, which will help to broaden the application of the wheat genome and can be used as a template for research on other polyploid plants.

Identification and characterization of large-scale genomic rearrangements during wheat evolution

Following allopolyploidization, nascent polyploid wheat species react with massive genomic rearrangements, including deletion of transposable element-containing sequences. While such massive rearrangements are considered to be a prominent process in wheat genome evolution and speciation, their structure, extent, and underlying mechanisms remain poorly understood. In this study, we retrieved ~3500 insertions of a specific variant of Fatima, one of the most dynamic gypsy long-terminal repeat retrotransposons in wheat from the recently available high-quality genome drafts of Triticum aestivum (bread wheat) and Triticum turgidum ssp. dicoccoides or wild emmer, the allotetraploid mother of all modern wheats. The dynamic nature of Fatima facilitated the identification of large (i.e., up to ~ 1 million bases) Fatima-containing insertions/deletions (indels) upon comparison of bread wheat and wild emmer genomes. We characterized 11 such indels using computer-assisted analysis followed by PCR validation, and found that they might have occurred via unequal intra-strand recombination or double-strand break (DSB) events. Additionally, we observed one case of introgression of novel DNA fragments from an unknown source into the wheat genome. Our data thus indicate that massive large-scale DNA rearrangements might play a prominent role in wheat speciation.

Genome-wide introgression from a bread wheat × Lophopyrum elongatum amphiploid into wheat

We introgressed wheatgrass germplasm from the octoploid amphiploid Triticum aestivum× Lophopyrum elongatum into wheat by manipulating the wheat Ph1 gene and discovered and characterized 130 introgression lines harboring single or, in various combinations, complete and recombined L. elongatum chromosomes. Diploid wheatgrass Lophopyrum elongatum (genomes EE) possesses valuable traits for wheat genetics and breeding. We evaluated several strategies for introgression of this germplasm into wheat. To detect it, we developed and validated multiplexed sets of Sequenom MassARRAY single nucleotide polymorphism (SNP) markers, which differentiated disomic and monosomic L. elongatum chromosomes from wheat chromosomes. We identified 130 introgression lines (ILs), which harbored 108 complete and 89 recombined L. elongatum chromosomes. Of the latter, 59 chromosomes were recombined by one or more crossovers and 30 were involved in centromeric (Robertsonian) translocations or were telocentric. To identify wheat chromosomes substituted for or recombined with L. elongatum chromosomes, we genotyped the ILs with the wheat 90-K Infinium SNP array. We found that most of the wheat 90-K probes correctly detected their targets in the L. elongatum genome and showed that some wheat SNPs are ancient and had originated prior to the divergence of the wheat and L. elongatum lineages. Of the 130 ILs, 52% were homozygous for Ph1 deletion and thus are staged to be recombined further. We failed to detect in the L. elongatum genome the 4/5 reciprocal translocation that has been reported in Thinopyrum bessarabicum and several other Triticeae genomes.

A roadmap for gene functional characterisation in crops with large genomes: Lessons from polyploid wheat

Understanding the function of genes within staple crops will accelerate crop improvement by allowing targeted breeding approaches. Despite their importance, a lack of genomic information and resources has hindered the functional characterisation of genes in major crops. The recent release of high-quality reference sequences for these crops underpins a suite of genetic and genomic resources that support basic research and breeding. For wheat, these include gene model annotations, expression atlases and gene networks that provide information about putative function. Sequenced mutant populations, improved transformation protocols and structured natural populations provide rapid methods to study gene function directly. We highlight a case study exemplifying how to integrate these resources. This review provides a helpful guide for plant scientists, especially those expanding into crop research, to capitalise on the discoveries made in Arabidopsis and other plants. This will accelerate the improvement of crops of vital importance for food and nutrition security.

Genomic analysis of Spanish wheat landraces reveals their variability and potential for breeding

Background

One of the main goals of the plant breeding in the twenty-first century is the development of crop cultivars that can maintain current yields in unfavorable environments. Landraces that have been grown under varying local conditions include genetic diversity that will be essential to achieve this objective. The Center of Plant Genetic Resources of the Spanish Institute for Agriculture Research maintains a broad collection of wheat landraces. These accessions, which are locally adapted to diverse eco-climatic conditions, represent highly valuable materials for breeding. However, their efficient use requires an exhaustive genetic characterization. The overall aim of this study was to assess the diversity and population structure of a selected set of 380 Spanish landraces and 52 reference varieties of bread and durum wheat by high-throughput genotyping.

Results

The DArTseq GBS approach generated 10 K SNPs and 40 K high-quality DArT markers, which were located against the currently available bread and durum wheat reference genomes. The markers with known locations were distributed across all chromosomes with relatively well-balanced genome-wide coverage. The genetic analysis showed that the Spanish wheat landraces were clustered in different groups, thus representing genetic pools providing a range of allelic variation. The subspecies had a major impact on the population structure of the durum wheat landraces, with three distinct clusters that corresponded to subsp. durum, turgidum and dicoccon being identified. The population structure of bread wheat landraces was mainly biased by geographic origin.

Conclusions

The results showed broader genetic diversity in the landraces compared to a reference set that included commercial varieties, and higher divergence between the landraces and the reference set in durum wheat than in bread wheat. The analyses revealed genomic regions whose patterns of variation were markedly different in the landraces and reference varieties, indicating loci that have been under selection during crop improvement, which could help to target breeding efforts. The results obtained from this work will provide a basis for future genome-wide association studies.

Genome-wide association mapping revealed syntenic loci QFhb-4AL and QFhb-5DL for Fusarium head blight resistance in common wheat (Triticum aestivum L.)

BACKGROUND

Fusarium head blight (FHB), primarily caused by Fusarium graminearum, is a major threat to wheat production and food security worldwide. Breeding stably and durably resistant cultivars is the most effective approach for managing and controlling the disease. The success of FHB resistance breeding relies on identification of an effective resistant germplasm. We conducted a genome-wide association study (GWAS) using the high-density wheat 90 K single nucleotide polymorphism (SNP) assays to better understand the genetic basis of FHB resistance in natural population and identify associated molecular markers.

RESULTS

The resistance to FHB fungal spread along the rachis (Type II resistance) was evaluated on 171 wheat cultivars in the 2016-2017 (abbr. as 2017) and 2017-2018 (abbr. as 2018) growing seasons. Using Illumina Infinum iSelect 90 K SNP genotyping data, a genome-wide association study (GWAS) identified 26 loci (88 marker-trait associations), which explained 6.65-14.18% of the phenotypic variances. The associated loci distributed across all chromosomes except 2D, 6A, 6D and 7D, with those on chromosomes 1B, 4A, 5D and 7A being detected in both years. New loci for Type II resistance were found on syntenic genomic regions of chromsome 4AL (QFhb-4AL, 621.85-622.24 Mb) and chromosome 5DL (QFhb-5DL, 546.09-547.27 Mb) which showed high collinearity in gene content and order. SNP markers wsnp_JD_c4438_5568170 and wsnp_CAP11_c209_198467 of 5D, reported previously linked to a soil-borne wheat mosaic virus (SBWMV) resistance gene, were also associated with FHB resistance in this study.

CONCLUSION

The syntenic FHB resistant loci and associated SNP markers identified in this study are valuable for FHB resistance breeding via marker-assisted selection.

Cloning, characterization of TaGS3 and identification of allelic variation associated with kernel traits in wheat (Triticum aestivum L.)

Background

Grain weight is an important yield component. Selection of advanced lines with heavy grains show high grain sink potentials and strong sink activity, which is an increasingly important objective in wheat breeding programs. Rice OsGS3 has been identified as a major quantitative trait locus for both grain weight and grain size. However, allelic variation of GS3 has not been characterized previously in hexaploid wheat.

Results

We cloned 2445, 2393, and 2409 bp sequences of the homologs TaGS3-4A, TaGS3-7A, and TaGS3-7D in wheat ‘Changzhi 6406’, a cultivar that shows high grain weight. The TaGS3 genes each contained five exons and four introns, and encoded a deduced protein of 170, 169, and 169 amino acids, respectively. Phylogenetic analysis of plant GS3 protein sequences revealed GS3 to be a monocotyledon-specific gene and the GS3 proteins were resolved into three classes. The length of the atypical Gγ domain and the cysteine-rich region was conserved within each class and not conserved between classes. A single-nucleotide polymorphism in the fifth exon (at position 1907) of TaGS3-7A leads to an amino acid change (ALA/THR) and showed different frequencies in two pools of Chinese wheat accessions representing extremes in grain weight. Association analysis indicated that the TaGS3-7A-A allele was associated with higher grain weight in the natural population. The TaGS3-7A-A allele was favoured in global modern wheat cultivars but the allelic frequency varied among different wheat-production regions of China, which indicated that this allele is of potential utility to improve wheat grain weight in certain wheat-production areas of China.

Conclusions

The novel molecular information on wheat GS3 homologs and the KASP functional marker designed in this study may be useful in marker-assisted breeding for genetic improvement of wheat.

Development of a New Am -Genome-Specific Single Nucleotide Polymorphism Marker Set for the Molecular Characterization of Wheat-Triticum monococcum Introgression Lines

We identified 1247 polymorphic single nucleotide polymorphisms between Triticum monococcum and wheat. We identified 191 markers validated across all seven chromosomes of T. monococcum. Detected a T. monococcum introgression in leaf-rust-resistant lines. Cultivated einkorn wheat (Triticum monococcum L. subsp. monococcum, 2n = 2x = 14, A^m A^m ) and its wild relative T. monococcum subsp. aegilopoides are important sources of economically useful genes that can be exploited for wheat (Triticum aestivum L.) breeding. Einkorn has excellent resistance to fungal diseases and gene transfer is relatively simple via standard breeding methods. To fulfill the growing demand by modern prebreeding programs for a cost-effective high-throughput procedure for accurately detecting introgressed chromosomes or chromosome segments from T. monococcum into wheat, we used the Axiom Wheat-Relative Genotyping Array and developed a set of A^m genome-specific exome-based single nucleotide polymorphism (SNP) markers suitable for rapid identification of T. monococcum chromatin in a wheat background. We identified 1247 polymorphic SNPs between T. monococcum and wheat. We identified 191 markers across all seven chromosomes of T. monococcum that are also present on an existing Triticum urartu Thum. ex Gandil. genetic map and potentially ordered them on the basis of the high macrocollinearity and conservation of marker order between T. monococcum and T. urartu. The marker set has been tested on leaf-rust-resistant BC₃ F₄ progenies of wheat-T. monococcum hybrids. Two markers (AX-94492165, AX-95073542) placed on the distal end of the chromosome arm 7AL detected a T. monococcum introgression into wheat. The SNP marker set thus proved highly effective in the identification of T. monococcum chromatin in a wheat background, offering a reliable method for screening and selecting wheat-T. monococcum introgression lines, a procedure that could significantly speed up prebreeding programs.

Genome-Wide Identification and Analysis of the AP2 Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

The AP2 transcription factors play important roles in regulating plant growth and development. However, limited data are available on the contributions of AP2 transcription factors in wheat (Triticum aestivum L.). In the present study, a total of 62 AP2 genes were identified in wheat from a genome-wide search against the latest wheat genome data. Phylogenetic and sequence alignment analyses divided the wheat AP2 genes into 3 clusters, euAP2, euANT, and basalANT. Chromosomal distribution, gene structure and duplication, and motif composition were subsequently investigated. The 62 TaAP2 genes were unevenly distributed on 21 chromosomes. Twenty-four homologous gene sets among A, B, and D sub-genomes were detected, which contributed to the expansion of the wheat AP2 gene family. The expression levels of TaAP2 genes were examined using the WheatExp database; most detected genes exhibited tissue-specific expression patterns. The transcript levels of 9 randomly selected TaAP2 genes were validated through qPCR analyses. Overexpression of TaAP2-10-5D, the most likely homolog of Arabidopsis ANT gene, increased organ sizes in Arabidopsis. Our results extend our knowledge of the AP2 gene family in wheat, and contribute to further functional characterization of AP2s during wheat development with the ultimate goal of improving crop production.

Unraveling the genetic architecture of grain size in einkorn wheat through linkage and homology mapping and transcriptomic profiling

Understanding the genetic architecture of grain size is a prerequisite to manipulating grain development and improving the potential crop yield. In this study, we conducted a whole genome-wide quantitative trait locus (QTL) mapping of grain-size-related traits by constructing a high-density genetic map using 109 recombinant inbred lines of einkorn wheat. We explored the candidate genes underlying QTLs through homologous analysis and RNA sequencing. The high-density genetic map spanned 1873 cM and contained 9937 single nucleotide polymorphism markers assigned to 1551 bins on seven chromosomes. Strong collinearity and high genome coverage of this map were revealed by comparison with physical maps of wheat and barley. Six grain size-related traits were surveyed in five environments. In total, 42 QTLs were identified; these were assigned to 17 genomic regions on six chromosomes and accounted for 52.3-66.7% of the phenotypic variation. Thirty homologous genes involved in grain development were located in 12 regions. RNA sequencing identified 4959 genes differentially expressed between the two parental lines. Twenty differentially expressed genes involved in grain size development and starch biosynthesis were mapped to nine regions that contained 26 QTLs, indicating that the starch biosynthesis pathway plays a vital role in grain development in einkorn wheat. This study provides new insights into the genetic architecture of grain size in einkorn wheat; identification of the underlying genes enables understanding of grain development and wheat genetic improvement. Furthermore, the map facilitates quantitative trait mapping, map-based cloning, genome assembly, and comparative genomics in wheat taxa.

Unraveling the genetic architecture of grain size in einkorn wheat through linkage and homology mapping and transcriptomic profiling

Understanding the genetic architecture of grain size is a prerequisite to manipulating grain development and improving the potential crop yield. In this study, we conducted a whole genome-wide quantitative trait locus (QTL) mapping of grain-size-related traits by constructing a high-density genetic map using 109 recombinant inbred lines of einkorn wheat. We explored the candidate genes underlying QTLs through homologous analysis and RNA sequencing. The high-density genetic map spanned 1873 cM and contained 9937 single nucleotide polymorphism markers assigned to 1551 bins on seven chromosomes. Strong collinearity and high genome coverage of this map were revealed by comparison with physical maps of wheat and barley. Six grain size-related traits were surveyed in five environments. In total, 42 QTLs were identified; these were assigned to 17 genomic regions on six chromosomes and accounted for 52.3-66.7% of the phenotypic variation. Thirty homologous genes involved in grain development were located in 12 regions. RNA sequencing identified 4959 genes differentially expressed between the two parental lines. Twenty differentially expressed genes involved in grain size development and starch biosynthesis were mapped to nine regions that contained 26 QTLs, indicating that the starch biosynthesis pathway plays a vital role in grain development in einkorn wheat. This study provides new insights into the genetic architecture of grain size in einkorn wheat; identification of the underlying genes enables understanding of grain development and wheat genetic improvement. Furthermore, the map facilitates quantitative trait mapping, map-based cloning, genome assembly, and comparative genomics in wheat taxa.

Comparative linkage mapping of diploid, tetraploid, and hexaploid Avena species suggests extensive chromosome rearrangement in ancestral diploids

The genus Avena (oats) contains diploid, tetraploid and hexaploid species that evolved through hybridization and polyploidization. Four genome types (named A through D) are generally recognized. We used GBS markers to construct linkage maps of A genome diploid (Avena strigosa x A. wiestii, 2n = 14), and AB genome tetraploid (A. barbata 2n = 28) oats. These maps greatly improve coverage from older marker systems. Seven linkage groups in the tetraploid showed much stronger homology and synteny with the A genome diploids than did the other seven, implying an allopolyploid hybrid origin of A. barbata from distinct A and B genome diploid ancestors. Inferred homeologies within A. barbata revealed that the A and B genomes are differentiated by several translocations between chromosomes within each subgenome. However, no translocation exchanges were observed between A and B genomes. Comparison to a consensus map of ACD hexaploid A. sativa (2n = 42) revealed that the A and D genomes of A. sativa show parallel rearrangements when compared to the A genomes of the diploids and tetraploids. While intergenomic translocations are well known in polyploid Avena, our results are most parsimoniously explained if translocations also occurred in the A, B and D genome diploid ancestors of polyploid Avena.

Genome mapping of quantitative trait loci (QTL) controlling domestication traits of intermediate wheatgrass (Thinopyrum intermedium)

Allohexaploid (2n = 6x = 42) intermediate wheatgrass (Thinopyrum intermedium), abbreviated IWG, is an outcrossing perennial grass belonging to the tertiary gene pool of wheat. Perenniality would be valuable option for grain production, but attempts to introgress this complex trait from wheat-Thinopyrum hybrids have not been commercially successful. Efforts to breed IWG itself as a dual-purpose forage and grain crop have demonstrated useful progress and applications, but grain yields are significantly less than wheat. Therefore, genetic and physical maps have been developed to accelerate domestication of IWG. Herein, these maps were used to identify quantitative trait loci (QTLs) and candidate genes associated with IWG grain production traits in a family of 266 full-sib progenies derived from two heterozygous parents, M26 and M35. Transgressive segregation was observed for 17 traits related to seed size, shattering, threshing, inflorescence capacity, fertility, stem size, and flowering time. A total of 111 QTLs were detected in 36 different regions using 3826 genotype-by-sequence markers in 21 linkage groups. The most prominent QTL had a LOD score of 15 with synergistic effects of 29% and 22% over the family means for seed retention and percentage of naked seeds, respectively. Many QTLs aligned with one or more IWG gene models corresponding to 42 possible domestication orthogenes including the wheat Q and RHT genes. A cluster of seed-size and fertility QTLs showed possible alignment to a putative Z self-incompatibility gene, which could have detrimental grain-yield effects when genetic variability is low. These findings elucidate pathways and possible hurdles in the domestication of IWG.

Genome-Wide Identification and Homoeologous Expression Analysis of PP2C Genes in Wheat (Triticum aestivum L.)

Plant protein phosphatase 2Cs (PP2Cs) play crucial roles in phytohormone signaling, developmental processes, and both biotic and abiotic stress responses. However, little research has been conducted on the PP2C gene family in hexaploid wheat (Triticum aestivum L.), which is an important cereal crop. In this study, a genome-wide investigation of TaPP2C gene family was performed. A total of 257 homoeologs of 95 TaPP2C genes were identified, of which 80% of genes had all the three homoeologs across A, B, and D subgenomes. Domain analysis indicated that all the TaPP2C homoeologs harbored the type 2C phosphatase domains. Based on the phylogenetic analysis, TaPP2Cs were divided into 13 groups (A-M) and 4 single branches, which corresponded to the results of gene structure and protein motif analyses. Results of chromosomal location and synteny relationship analysis of TaPP2C homoeologs revealed that known chromosome translocation events and pericentromeric inversions were responsible for the formation of TaPP2C gene family. Expression patterns of TaPP2C homoeologs in various tissues and under diverse stress conditions were analyzed using publicly available RNA-seq data. The results suggested that TaPP2C genes regulate wheat developmental processes and stress responses. Homoeologous expression patterns of TaPP2C triad homoeologs from A, B, and D subgenomes, revealed expression bias within triads under the normal condition, and variability in expression under different stress treatments. Quantitative real-time PCR (qRT-PCR) analysis of eight TaPP2C genes in group A revealed that they were all up-regulated after abscisic acid treatment. Some genes in group A also responded to other phytohormones such as methyl jasmonate and gibberellin. Yeast two-hybrid assays showed that group A TaPP2Cs also interacted with TaSnRK2.1 and TaSnRK2.2 from subclass II, besides with subclass III TaSnRK2s. TaPP2C135 in group A was transformed into Arabidopsis and germination assay revealed that ectopic expression of TaPP2C135 in Arabidopsis enhanced its tolerance to ABA. Overall, these results enhance our understanding of the function of TaPP2Cs in wheat, and provide novel insights into the roles of group A TaPP2Cs. This information will be useful for in-depth functional analysis of TaPP2Cs in future studies and for wheat breeding.

Genome-wide polymorphisms from RNA sequencing assembly of leaf transcripts facilitate phylogenetic analysis and molecular marker development in wild einkorn wheat

A survey of genome-wide polymorphisms between closely related species is required to understand the molecular basis of the evolutionary differentiation of their genomes. Two wild diploid wheat species, namely Triticum monococcum ssp. aegilopoides and T. urartu, are closely related and harbour the A^m and A genomes, respectively. The A-genome donor of tetraploid and common wheat is T. urartu, and T. monococcum ssp. monococcum is the cultivated form derived from the wild einkorn wheat subspecies aegilopoides. Although subspecies aegilopoides has been a useful genetic resource in wheat breeding, genome-wide molecular markers for this subspecies have not been sufficiently developed. Here, we describe the detection of genome-wide polymorphisms such as single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) from RNA sequencing (RNA-seq) data of leaf transcripts in 15 accessions of the two diploid wheat species. The SNPs and indels, detected using the A genome of common wheat as the reference genome, covered the entire chromosomes of these species. The polymorphism information facilitated a comparison of the genetic diversity of einkorn wheat with that of two related diploid Aegilops species, namely, Ae. tauschii and Ae. umbellulata. Cleaved amplified polymorphic sequence (CAPS) markers converted from the SNP data were efficiently developed to confirm the addition of aegilopoides subspecies chromosomes to tetraploid wheat in nascent allohexaploid lines with AABBA^mA^m genomes. In addition, the CAPS markers permitted linkage map construction in mapping populations of aegilopoides subspecies accessions. Therefore, these RNA-seq data provide information for further breeding of closely related species with no reference genome sequence data.

Development and characterisation of interspecific hybrid lines with genome-wide introgressions from Triticum timopheevii in a hexaploid wheat background

BACKGROUND

Triticum timopheevii (2n = 4x = 28; AtAtGG), is an important source for new genetic variation for wheat improvement with genes for potential disease resistance and salt tolerance. By generating a range of interspecific hybrid lines, T. timopheevii can contribute to wheat's narrow gene-pool and be practically utilised in wheat breeding programmes. Previous studies that have generated such introgression lines between wheat and its wild relatives have been unable to use high-throughput methods to detect the presence of wild relative segments in such lines.

RESULTS

A whole genome introgression approach, exploiting homoeologous recombination in the absence of the Ph1 locus, has resulted in the transfer of different chromosome segments from both the At and G genomes of T. timopheevii into wheat. These introgressions have been detected and characterised using single nucleotide polymorphism (SNP) markers present on a high-throughput Axiom® Genotyping Array. The analysis of these interspecific hybrid lines has resulted in the detection of 276 putative unique introgressions from T. timopheevii, thereby allowing the generation of a genetic map of T. timopheevii containing 1582 SNP markers, spread across 14 linkage groups representing each of the seven chromosomes of the At and G genomes of T. timopheevii. The genotyping of the hybrid lines was validated through fluorescence in situ hybridisation (FISH). Comparative analysis of the genetic map of T. timopheevii and the physical map of the hexaploid wheat genome showed that synteny between the two species is highly conserved at the macro-level and confirmed the presence of inter- and intra-genomic translocations within the At and G genomes of T. timopheevii that have been previously only detected through cytological techniques.

CONCLUSIONS

In this work, we report a set of SNP markers present on a high-throughput genotyping array, able to detect the presence of T. timopheevii in a hexaploid wheat background making it a potentially valuable tool for marker assisted selection (MAS) in wheat pre-breeding programs. These valuable resources of high-density molecular markers and wheat-T. timopheevii hybrid lines will greatly enhance the work being undertaken for wheat improvement through wild relative introgressions.

Development and validation of an exome-based SNP marker set for identification of the St, Jr and Jvs genomes of Thinopyrym intermedium in a wheat background

KEY MESSAGE

Cytogenetic analysis and array-based SNP genotyping of wheat- Th. intermedium introgression lines allowed identification of 634 chromosome-specific SNP markers across all twenty-one chromosomes of Th. intermedium (StJ r J vs , 2 n  = 6 x  = 42). Thinopyrum intermedium (2n = 6x = 42, StJrJvs) is one of the most promising reservoirs of useful genes including tolerance to abiotic stresses, perenniality and disease resistance not available in the cultivated bread wheat. The transfer of genetic diversity from wild species to wheat offers valuable responses to the effects of climate change. The new array-based single-nucleotide polymorphism (SNP) marker technology provides cheap and easy-to-use molecular markers for marker-assisted selection (MAS) in wheat breeding programmes. Here, we focus on the generation of a new chromosome-specific SNP marker set that can be used to characterize and identify the Th. intermedium chromosomes or chromosome segments transferred into wheat. A progressive investigation of marker development was conducted using 187 various newly developed wheat-Th. intermedium introgression lines and the Axiom® Wheat-Relative Genotyping array. We employed molecular cytogenetic techniques to clarify the genome constitution of the Th. intermedium parental lines and validated 634 chromosome-specific SNPs. Our data confirmed the allohexaploid nature of Th. intermedium and demonstrated that the St genome-specific GISH signal and markers are present at the centromeric regions of chromosomes 1Jvs, 2Jvs, 3Jvs and 7Jvs. The SNP markers presented here will be introduced into current wheat improvement programmes, offering a significant speed-up in wheat breeding and making it possible to deal with the transfer of the full genetic potential of Th. intermedium into wheat.