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.

Bread wheat (T. aestivum) variability: Phenotypic and genotypic data from 75 varieties

Most bread wheat is consumed after processing, which mainly depends on the quantity and quality of protein in the grain. Storage protein content and composition particularly influence the end use quality of milled grain products. Storage proteins are components of the gluten network that confer dough viscoelasticity, an essential property for processing. To explore grain storage protein diversity, 75 bread wheat accessions were grown with two replicates each at two locations. Grains were harvested at maturity and samples were phenotyped for each site and each replicate plant. Grain hardness, thousand-kernel weight and grain nitrogen content were measured. The protein composition of flour from each replicate was characterised by reverse phase-high performance liquid chromatography (RP-HPLC). The molecular distribution of flour polymers was determined by asymmetric flow field-flow fractionation (AF4) and dough technological properties were assessed using a Glutomatic system and a Chopin alveograph. In addition, the 75 accessions were genotyped by the BreedWheat 35k genotyping array (Axiom TaBW35K) containing 34,746 single nucleotide polymorphism markers (SNPs). The dataset produced by this work includes six files with raw data, two files with protocols and figures. Data show the genotypic and phenotypic variabilities of the material used and can be used to explore genetic and environmental effects on traits involved in grain protein quality. This dataset is associated to the research article "Differences in bread protein digestibility traced to wheat cultivar traits" [1].

Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number

Understanding meiotic crossover (CO) variation in crops like bread wheat (Triticum aestivum L.) is necessary as COs are essential to create new, original and powerful combinations of genes for traits of agronomical interest. We cytogenetically characterized a set of wheat aneuploid lines missing part or all of chromosome 3B to identify the most influential regions for chiasma formation located on this chromosome. We showed that deletion of the short arm did not change the total number of chiasmata genome-wide, whereas this latter was reduced by ~35% while deleting the long arm. Contrary to what was hypothesized in a previous study, deletion of the long arm does not disturb the initiation of the synaptonemal complex (SC) in early meiotic stages. However, progression of the SC is abnormal, and we never observed its completion when the long arm is deleted. By studying six different deletion lines (missing different parts of the long arm), we revealed that at least two genes located in both the proximal (C-3BL2-0.22) and distal (3BL7-0.63-1.00) deletion bins are involved in the control of chiasmata, each deletion reducing the number of chiasmata by ~15%. We combined sequence analyses of deletion bins with RNA-Seq data derived from meiotic tissues and identified a set of genes for which at least the homoeologous copy on chromosome 3B is expressed and which are involved in DNA processing. Among these genes, eight (CAP-E1/E2, DUO1, MLH1, MPK4, MUS81, RTEL1, SYN4, ZIP4) are known to be involved in the recombination pathway.

Long-read and chromosome-scale assembly of the hexaploid wheat genome achieves high resolution for research and breeding

BACKGROUND

The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost.

RESULTS

Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan.

CONCLUSIONS

We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

SPO11.2 is essential for programmed double-strand break formation during meiosis in bread wheat (Triticum aestivum L.)

Meiotic recombination is initiated by formation of DNA double-strand breaks (DSBs). This involves a protein complex that includes in plants the two similar proteins, SPO11-1 and SPO11-2. We analysed the sequences of SPO11-2 in hexaploid bread wheat (Triticum aestivum), as well as in its diploid and tetraploid progenitors. We investigated its role during meiosis using single, double and triple mutants. The three homoeologous SPO11-2 copies of hexaploid wheat exhibit high nucleotide and amino acid similarities with those of the diploids, tetraploids and Arabidopsis. Interestingly, however, two nucleotides deleted in exon-2 of the A copy lead to a premature stop codon and suggest that it encodes a non-functional protein. Remarkably, the mutation was absent from the diploid A-relative Triticum urartu, but present in the tetraploid Triticum dicoccoides and in different wheat cultivars indicating that the mutation occurred after the first polyploidy event and has since been conserved. We further show that triple mutants with all three copies (A, B, D) inactivated are sterile. Cytological analyses of these mutants show synapsis defects, accompanied by severe reductions in bivalent formation and numbers of DMC1 foci, thus confirming the essential role of TaSPO11-2 in meiotic recombination in wheat. In accordance with its 2-nucleotide deletion in exon-2, double mutants for which only the A copy remained are also sterile. Notwithstanding, some DMC1 foci remain visible in this mutant, suggesting a residual activity of the A copy, albeit not sufficient to restore fertility.

Structural Variations Affecting Genes and Transposable Elements of Chromosome 3B in Wheats

Structural variations (SVs) such as copy number and presence–absence variations are polymorphisms that are known to impact genome composition at the species level and are associated with phenotypic variations. In the absence of a reference genome sequence, their study has long been hampered in wheat. The recent production of new wheat genomic resources has led to a paradigm shift, making possible to investigate the extent of SVs among cultivated and wild accessions. We assessed SVs affecting genes and transposable elements (TEs) in a Triticeae diversity panel of 45 accessions from seven tetraploid and hexaploid species using high-coverage shotgun sequencing of sorted chromosome 3B DNA and dedicated bioinformatics approaches. We showed that 23% of the genes are variable within this panel, and we also identified 330 genes absent from the reference accession Chinese Spring. In addition, 60% of the TE-derived reference markers were absent in at least one accession, revealing a high level of intraspecific and interspecific variability affecting the TE space. Chromosome extremities are the regions where we observed most of the variability, confirming previous hypotheses made when comparing wheat with the other grasses. This study provides deeper insights into the genomic variability affecting the complex Triticeae genomes at the intraspecific and interspecific levels and suggests a phylogeny with independent hybridization events leading to different hexaploid species.

High throughput genotyping of structural variations in a complex plant genome using an original Affymetrix® axiom® array

Background

Insertions/deletions (InDels) and more specifically presence/absence variations (PAVs) are pervasive in several species and have strong functional and phenotypic effect by removing or drastically modifying genes. Genotyping of such variants on large panels remains poorly addressed, while necessary for approaches such as association mapping or genomic selection.

Results

We have developed, as a proof of concept, a new high-throughput and affordable approach to genotype InDels. We first identified 141,000 InDels by aligning reads from the B73 line against the genome of three temperate maize inbred lines (F2, PH207, and C103) and reciprocally. Next, we designed an Affymetrix® Axiom® array to target these InDels, with a combination of probes selected at breakpoint sites (13%) or within the InDel sequence, either at polymorphic (25%) or non-polymorphic sites (63%) sites. The final array design is composed of 662,772 probes and targets 105,927 InDels, including PAVs ranging from 35 bp to 129kbp. After Affymetrix® quality control, we successfully genotyped 86,648 polymorphic InDels (82% of all InDels interrogated by the array) on 445 maize DNA samples with 422,369 probes. Genotyping InDels using this approach produced a highly reliable dataset, with low genotyping error (~ 3%), high call rate (~ 98%), and high reproducibility (> 95%). This reliability can be further increased by combining genotyping of several probes calling the same InDels (< 0.1% error rate and > 99.9% of call rate for 5 probes). This “proof of concept” tool was used to estimate the kinship matrix between 362 maize lines with 57,824 polymorphic InDels. This InDels kinship matrix was highly correlated with kinship estimated using SNPs from Illumina 50 K SNP arrays.

Conclusions

We efficiently genotyped thousands of small to large InDels on a sizeable number of individuals using a new Affymetrix® Axiom® array. This powerful approach opens the way to studying the contribution of InDels to trait variation and heterosis in maize. The approach is easily extendable to other species and should contribute to decipher the biological impact of InDels at a larger scale.

Worldwide phylogeography and history of wheat genetic diversity

Since its domestication in the Fertile Crescent ~8000 to 10,000 years ago, wheat has undergone a complex history of spread, adaptation, and selection. To get better insights into the wheat phylogeography and genetic diversity, we describe allele distribution through time using a set of 4506 landraces and cultivars originating from 105 different countries genotyped with a high-density single-nucleotide polymorphism array. Although the genetic structure of landraces is collinear to ancient human migration roads, we observe a reshuffling through time, related to breeding programs, with the appearance of new alleles enriched with structural variations that may be the signature of introgressions from wild relatives after 1960.

Linking the International Wheat Genome Sequencing Consortium bread wheat reference genome sequence to wheat genetic and phenomic data

The Wheat@URGI portal has been developed to provide the international community of researchers and breeders with access to the bread wheat reference genome sequence produced by the International Wheat Genome Sequencing Consortium. Genome browsers, BLAST, and InterMine tools have been established for in-depth exploration of the genome sequence together with additional linked datasets including physical maps, sequence variations, gene expression, and genetic and phenomic data from other international collaborative projects already stored in the GnpIS information system. The portal provides enhanced search and browser features that will facilitate the deployment of the latest genomics resources in wheat improvement.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Exploiting the Repetitive Fraction of the Wheat Genome for High-Throughput Single-Nucleotide Polymorphism Discovery and Genotyping

Transposable elements (TEs) account for more than 80% of the wheat genome. Although they represent a major obstacle for genomic studies, TEs are also a source of polymorphism and consequently of molecular markers such as insertion site-based polymorphism (ISBP) markers. Insertion site-based polymorphisms have been found to be a great source of genome-specific single-nucleotide polymorphism (SNPs) in the hexaploid wheat ( L.) genome. Here, we report on the development of a high-throughput SNP discovery approach based on sequence capture of ISBP markers. By applying this approach to the reference sequence of chromosome 3B from hexaploid wheat, we designed 39,077 SNPs that are evenly distributed along the chromosome. We demonstrate that these SNPs can be efficiently scored with the KASPar (Kompetitive allele-specific polymerase chain reaction) genotyping technology. Finally, through genetic diversity and genome-wide association studies, we also demonstrate that ISBP-derived SNPs can be used in marker-assisted breeding programs.