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.

The transcriptional landscape of polyploid wheat

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

A physical map of the 1-gigabase bread wheat chromosome 3B

As the staple food for 35% of the world's population, wheat is one of the most important crop species. To date, sequence-based tools to accelerate wheat improvement are lacking. As part of the international effort to sequence the 17-billion-base-pair hexaploid bread wheat genome (2n = 6x = 42 chromosomes), we constructed a bacterial artificial chromosome (BAC)-based integrated physical map of the largest chromosome, 3B, that alone is 995 megabases. A chromosome-specific BAC library was used to assemble 82% of the chromosome into 1036 contigs that were anchored with 1443 molecular markers, providing a major resource for genetic and genomic studies. This physical map establishes a template for the remaining wheat chromosomes and demonstrates the feasibility of constructing physical maps in large, complex, polyploid genomes with a chromosome-based approach.

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.

Molecular-genetic maps for group 1 chromosomes of Triticeae species and their relation to chromosomes in rice and oat

Group 1 chromosomes of the Triticeae tribe have been studied extensively because many important genes have been assigned to them. In this paper, chromosome 1 linkage maps of Triticum aestivum, T. tauschii, and T. monococcum are compared with existing barley and rye maps to develop a consensus map for Triticeae species and thus facilitate the mapping of agronomic genes in this tribe. The consensus map that was developed consists of 14 agronomically important genes, 17 DNA markers that were derived from known-function clones, and 76 DNA markers derived from anonymous clones. There are 12 inconsistencies in the order of markers among seven wheat, four barley, and two rye maps. A comparison of the Triticeae group 1 chromosome consensus map with linkage maps of homoeologous chromosomes in rice indicates that the linkage maps for the long arm and the proximal portion of the short arm of group 1 chromosomes are conserved among these species. Similarly, gene order is conserved between Triticeae chromosome 1 and its homoeologous chromosome in oat. The location of the centromere in rice and oat chromosomes is estimated from its position in homoeologous group 1 chromosomes of Triticeae.

Megabase level sequencing reveals contrasted organization and evolution patterns of the wheat gene and transposable element spaces

To improve our understanding of the organization and evolution of the wheat (Triticum aestivum) genome, we sequenced and annotated 13-Mb contigs (18.2 Mb) originating from different regions of its largest chromosome, 3B (1 Gb), and produced a 2x chromosome survey by shotgun Illumina/Solexa sequencing. All regions carried genes irrespective of their chromosomal location. However, gene distribution was not random, with 75% of them clustered into small islands containing three genes on average. A twofold increase of gene density was observed toward the telomeres likely due to high tandem and interchromosomal duplication events. A total of 3222 transposable elements were identified, including 800 new families. Most of them are complete but showed a highly nested structure spread over distances as large as 200 kb. A succession of amplification waves involving different transposable element families led to contrasted sequence compositions between the proximal and distal regions. Finally, with an estimate of 50,000 genes per diploid genome, our data suggest that wheat may have a higher gene number than other cereals. Indeed, comparisons with rice (Oryza sativa) and Brachypodium revealed that a high number of additional noncollinear genes are interspersed within a highly conserved ancestral grass gene backbone, supporting the idea of an accelerated evolution in the Triticeae lineages.

New Secale cereale (rye) DNA derivatives for the detection of rye chromosome segments in wheat

Subcloning of a clone of the 120-bp family of rye, pSc119, has produced two extremely useful probes. pSc119.1 assays rye-specific dispersed repetitive sequence families. It is present on all seven rye chromosomes and hybridizes to the entire length of each chromosome, with the exception of some telomeres and the nucleolar organiser region. pSc119.2, in contrast, hybridizes predominantly to the telomeric regions of rye chromosomes, with some interstitial sites. Unlike pSc119.1, it assays similar repetitive sequence families in both wheat and rye chromosomes.

A workshop report on wheat genome sequencing: International Genome Research on Wheat Consortium

Sponsored by the National Science Foundation and the U.S. Department of Agriculture, a wheat genome sequencing workshop was held November 10-11, 2003, in Washington, DC. It brought together 63 scientists of diverse research interests and institutions, including 45 from the United States and 18 from a dozen foreign countries (see list of participants at http://www.ksu.edu/igrow). The objectives of the workshop were to discuss the status of wheat genomics, obtain feedback from ongoing genome sequencing projects, and develop strategies for sequencing the wheat genome. The purpose of this report is to convey the information discussed at the workshop and provide the basis for an ongoing dialogue, bringing forth comments and suggestions from the genetics community.

Resequencing of 145 Landmark Cultivars Reveals Asymmetric Sub-genome Selection and Strong Founder Genotype Effects on Wheat Breeding in China

Controlled pedigrees and the multi-decade timescale of national crop plant breeding programs offer a unique experimental context for examining how selection affects plant genomes. More than 3000 wheat cultivars have been registered, released, and documented since 1949 in China. In this study, a set of 145 elite cultivars selected from historical points of wheat breeding in China were re-sequenced. A total of 43.75 Tb of sequence data were generated with an average read depth of 17.94× for each cultivar, and more than 60.92 million SNPs and 2.54 million InDels were captured, based on the Chinese Spring RefSeq genome v1.0. Seventy years of breeder-driven selection led to dramatic changes in grain yield and related phenotypes, with distinct genomic regions and phenotypes targeted by different breeders across the decades. There are very clear instances illustrating how introduced Italian and other foreign germplasm was integrated into Chinese wheat programs and reshaped the genomic landscape of local modern cultivars. Importantly, the resequencing data also highlighted significant asymmetric breeding selection among the three sub-genomes: this was evident in both the collinear blocks for homeologous chromosomes and among sets of three homeologous genes. Accumulation of more newly assembled genes in newer cultivars implied the potential value of these genes in breeding. Conserved and extended sharing of linkage disequilibrium (LD) blocks was highlighted among pedigree-related cultivars, in which fewer haplotype differences were detected. Fixation or replacement of haplotypes from founder genotypes after generations of breeding was related to their breeding value. Based on the haplotype frequency changes in LD blocks of pedigree-related cultivars, we propose a strategy for evaluating the breeding value of any given line on the basis of the accumulation (pyramiding) of beneficial haplotypes. Collectively, our study demonstrates the influence of "founder genotypes" on the output of breeding efforts over many decades and also suggests that founder genotype perspectives are in fact more dynamic when applied in the context of modern genomics-informed breeding.

Frequent intra- and inter-species introgression shapes the landscape of genetic variation in bread wheat

BACKGROUND

Bread wheat is one of the most important and broadly studied crops. However, due to the complexity of its genome and incomplete genome collection of wild populations, the bread wheat genome landscape and domestication history remain elusive.

RESULTS

By investigating the whole-genome resequencing data of 93 accessions from worldwide populations of bread wheat and its diploid and tetraploid progenitors, together with 90 published exome-capture data, we find that the B subgenome has more variations than A and D subgenomes, including SNPs and deletions. Population genetics analyses support a monophyletic origin of domesticated wheat from wild emmer in northern Levant, with substantial introgressed genomic fragments from southern Levant. Southern Levant contributes more than 676 Mb in AB subgenomes and enriched in the pericentromeric regions. The AB subgenome introgression happens at the early stage of wheat speciation and partially contributes to their greater genetic diversity. Furthermore, we detect massive alien introgressions that originated from distant species through natural and artificial hybridizations, resulting in the reintroduction of ~ 709 Mb and ~ 1577 Mb sequences into bread wheat landraces and varieties, respectively. A large fraction of these intra- and inter-introgression fragments are associated with quantitative trait loci of important traits, and selection events are also identified.

CONCLUSION

We reveal the significance of multiple introgressions from distant wild populations and alien species in shaping the genetic components of bread wheat, and provide important resources and new perspectives for future wheat breeding.

The wheat ribosomal DNA spacer region: Its structure and variation in populations and among species

The wheat rDNA clone pTA250 was examined in detail to provide a restriction enzyme map and the nucleotide sequence of two of the eleven, 130 bp repeating units found within the spacer region. The 130 bp units showed some sequence heterogeneity. The sequence difference between the two 130 bp units analysed (130.6 and 130.8) was at 7 positions and could be detected as a 4 °C shift in Tm when heterologous and homologous hybrids were compared. This corresponded to a 1.2% change in nucleotide sequence per ΔTm of 1 °C. The sensitivity of the Tm analysis using cloned sequences facilitated the analysis of small sequence variations in the spacer region of different Triticum aestivum cultivars and natural populations of T. turgidum ssp. dicoccoides (referred to as T. dicoccoides). In addition spacer length variation was assayed by restriction enzyme digestion and hybridization with spacer sequence probes.Extensive polymorphism was observed for the spacer region in various cultivars of T. aestivum, although within each cultivar the rDNA clusters were homogeneous and could be assigned to particular chromosomes. Within natural populations of T. dicoccoides polymorphism was also observed but, once again, within any one individual the rDNA clusters appeared to be homogeneous. The polymorphism, at the sequence level (assayed by Tm analysis), was not so great as to prevent the use of spacer sequence variation as a probe for evolutionary relationships. The length variation as assayed by restriction enzyme digestion did not appear to be as useful in this regard, since its range of variation was extensive even within populations of a species.

The use of cell culture for subchromosomal introgressions of barley yellow dwarf virus resistance from Thinopyrum intermedium to wheat

Barley yellow dwarf virus (BYDV) resistance has been transferred to wheat from a group 7 chromosome of Thinopyrum (Agropyron) intermedium. The source of the resistance gene was the L1 disomic addition line, which carries the 7Ai-1 chromosome. The resistance locus is on the long arm of this chromosome. BYDV resistant recombinant lines were identified after three or more generations of selection against a group 7 Th. intermedium short arm marker (red coleoptile) and selection for the presence of BYDV resistance. One recombinant line produced by ph. mutant induced homoeologous pairing and 14 recombinant lines induced by cell culture have been identified. Resistance in seven of the cell culture induced recombinants has been inherited via pollen according to Mendelian segregation ratios for up to eight generations. Meiotic analysis of heterozygotes indicates that the alien chromatin in the cell culture induced recombinants is small enough to allow regular meiotic behaviour. The ph-induced recombinant was less regular in meiosis. A probe, pEleAcc2, originally isolated from Th. elongatum and that hybridizes to dispersed repeated DNA sequences, was utilised to detect Th. intermedium chromatin, which confers resistance to BYDV, in wheat backgrounds. Quantification of these hybridization signals indicated that the translocations involved a portion of alien chromatin that was smaller than the complete long arm of 7Ai-1. Restriction fragment length polymorphism analysis confirmed the loss of the short arm of 7Ai-1 and indicated the retention of segments of the long arm of 7Ai-1. Two 7Ai-1L DNA markers always assorted with the BYDV resistance. A third 7Ai-IL DNA marker was also present in seven of eight recombinants. In all recombinants except TC7, the 7Ai-1L markers replaced the 7DL markers. None of the wheat group 7 markers was missing from TC7. It is concluded that all the resistant lines are the result of recombination with wheat chromosome 7D, except line TC7, which is the result of recombination with an unidentified nongroup 7 chromosome.

Proteome and phosphoproteome characterization reveals new response and defense mechanisms of Brachypodium distachyon leaves under salt stress

Salinity is a major abiotic stress affecting plant growth and development. Understanding the molecular mechanisms of salt response and defense in plants will help in efforts to improve the salt tolerance of crops. Brachypodium distachyon is a new model plant for wheat, barley, and several potential biofuel grasses. In the current study, proteome and phosphoproteome changes induced by salt stress were the focus. The Bd21 leaves were initially treated with salt in concentrations ranging from 80 to 320 mm and then underwent a recovery process prior to proteome analysis. A total of 80 differentially expressed protein spots corresponding to 60 unique proteins were identified. The sample treated with a median salt level of 240 mm and the control were selected for phosphopeptide purification using TiO2 microcolumns and LC-MS/MS for phosphoproteome analysis to identify the phosphorylation sites and phosphoproteins. A total of 1509 phosphoproteins and 2839 phosphorylation sites were identified. Among them, 468 phosphoproteins containing 496 phosphorylation sites demonstrated significant changes at the phosphorylation level. Nine phosphorylation motifs were extracted from the 496 phosphorylation sites. Of the 60 unique differentially expressed proteins, 14 were also identified as phosphoproteins. Many proteins and phosphoproteins, as well as potential signal pathways associated with salt response and defense, were found, including three 14-3-3s (GF14A, GF14B, and 14-3-3A) for signal transduction and several ABA signal-associated proteins such as ABF2, TRAB1, and SAPK8. Finally, a schematic salt response and defense mechanism in B. distachyon was proposed.

Trends in genetic and genome analyses in wheat: a review

The size and structure of the wheat genome makes it one of the most complex crop species for genetic analysis. The development of molecular techniques for genetic analysis, in particular the use of molecular markers to monitor DNA sequence variation between varieties, landraces, and wild relatives of wheat and related grass species, has led to a dramatic expansion in our understanding of wheat genetics and the structure and behaviour of the wheat genome. This review provides an overview of these developments, examines some of the special issues that have arisen in applying molecular techniques to genetic studies in wheat, and looks at the applications of these technologies to wheat breeding and to improving our understanding of the genetic basis of traits such as disease resistance and processing quality. The review also attempts to foreshadow some of the key molecular issues and developments that may occur in wheat genetics and breeding over the next few years.

Quantitative trait loci on chromosome 4B for coleoptile length and early vigour in wheat (Triticum aestivum L.)

The Norin-10 dwarfing genes, Rht-B1b (Rht1) and Rht-D1b (Rht2), are commonly used to reduce plant height and increase grain yield in wheat breeding programs. These dwarfing genes lower sensitivity of vegetative tissue to endogenous gibberellin to reduce cell and subsequent stem elongation. This reduction in cell elongation capacity reportedly results in a concomitant reduction in coleoptile length and early vigour (leaf area) thereby affecting seedling establishment and growth. A detailed genetic map from a cross between tall Halberd (Rht-B1a) and semidwarf Cranbrook (Rht-B1b) wheat cultivars was used to assess genetic factors affecting seedling growth. Parental and 150 doubled haploid progeny lines were characterised for seedling and height-related traits in controlled and field environments. Genotypic variation was large and predominantly under additive genetic control with evidence for transgressive segregation for some traits. Narrow-sense heritabilities were moderate to high (h2 = 0.31–0.91) indicating a strong genetic basis for differences between progeny. Molecular marker analyses identified a number of significant (P < 0.05) quantitative trait loci (QTL) for each trait. A major QTL, mapping directly to the Rht-B1 locus on chromosome arm 4BS, accounted for up to 49% of the genotypic variance in peduncle length and plant height, and 27–45% of the genotypic variance in coleoptile length across different temperatures. Another QTL, located close to the RFLP marker XksuC2 on the long arm of chromosome 4B, accounted for 15–27% of the genotypic variance in coleoptile length. The influence of the XksuC2-linked QTL on coleoptile length was greatest at 19˚C and decreased with cooler temperatures. The same QTL affected reductions in leaf size, and both coleoptile tiller size and presence to affect overall seedling vigour. There was also some evidence for epistatic interactions influencing coleoptile tiller growth. Reductions in plant size at the Rht-B1b and XksuC2 loci were associated with presence of the Cranbrook 4B allele. The negative genetic effect of the Rht-B1b dwarfing gene on early growth of wheat confirms phenotypic evidence of a pleiotropic effect of Rht-B1b on establishment and early vigour. Genetic increases in coleoptile length and early leaf area development are likely to be limited in wheat populations containing the Rht-B1b dwarfing gene.

Disomic Thinopyrum intermedium addition lines in wheat with barley yellow dwarf virus resistance and with rust resistances

Zhong 5 is a partial amphiploid (2n = 56) between Triticum aestivum (2n = 42) and Thinopyrum intermedium (2n = 42) carrying all the chromosomes of wheat and seven pairs of chromosomes from Th. intermedium. Following further backcrossing to wheat, six independent stable 2n = 44 lines were obtained representing 4 disomic chromosome addition lines. One chromosome confers barley yellow dwarf virus (BYDV) resistance, whereas two other chromosomes carry leaf and stem rust resistance; one of the latter also confers stripe rust resistance. Using RFLP and isozyme markers we have shown that the extra chromosome in the Zhong 5-derived BYDV resistant disomic addition lines (Z1, Z2, or Z6) belongs to the homoeologous group 2. It therefore carries a different locus to the BYDV resistant group 7 addition, L1, described previously. The leaf, stem, and stripe rust resistant line (Z4) carries an added group 7 chromosome. The line Z3 has neither BYDV nor rust resistance, is not a group 2 or group 7 addition, and is probably a group 1 addition. The line Z5 is leaf and stem rust resistant, is not stripe rust resistant, and its homoeology remains unknown.Key words: Agropyron, intermediate wheatgrass, leaf rust, stem rust, stripe rust, luteovirus.

Functional markers in wheat: current status and future prospects

Functional markers (FM) are developed from sequence polymorphisms present in allelic variants of a functional gene at a locus. FMs accurately discriminate alleles of a targeted gene, and are ideal molecular markers for marker-assisted selection in wheat breeding. In this paper, we summarize FMs developed and used in common wheat. To date, more than 30 wheat loci associated with processing quality, agronomic traits, and disease resistance, have been cloned, and 97 FMs were developed to identify 93 alleles based on the sequences of those genes. A general approach is described for isolation of wheat genes and development of FMs based on in silico cloning and comparative genomics. The divergence of DNA sequences of different alleles that affect gene function is summarized. In addition, 14 molecular markers specific for alien genes introduced from common wheat relatives were also described. This paper provides updated information on all FMs and gene-specific STS markers developed so far in wheat and should facilitate their application in wheat breeding programs.

The genetic analysis of D. melanogaster heterochromatin

Classical, or constitutive (Brown, 1966), heterochromatin is a nearly ubiquitous aspect of higher eukaryotic chromosome structure and its genetic significance has been debated since its first description (Heitz, 1928, 1929). This review describes genetic studies on the heterochromatin of Drosophila and the perspectives they provide on the possible biological role(s) of constitutive heterochromatin. In addition to documenting the existence of gene loci within heterochromatin, these studies provide evidence for long-range interactions within the nucleus and involving heterochromatin. These interactions may reflect more general properties of chromatin related to the three-dimensional structure of the interphase nucleus and its effects on gene expression.

Map-based cloning of a gene sequence encoding a nucleotide-binding domain and a leucine-rich region at the Cre3 nematode resistance locus of wheat

The Cre3 gene confers a high level of resistance to the root endoparasitic nematode Heterodera avenae in wheat. A DNA marker cosegregating with H. avenae resistance was used as an entry point for map-based cloning of a disease resistance gene family at the Cre3 locus. Two related gene sequences have been analysed at the Cre3 locus. One, identified as a cDNA clone, encodes a polypeptide with a nucleotide binding site (NBS) and a leucine-rich region; this member of the disease resistance gene family is expressed in roots. A second Cre3 gene sequence, cloned as genomic DNA, appears to be a pseudogene, with a frame shift caused by a deletion event. These two genes, related to members of the cytoplasmic NBS-leucine rich repeat class of plant disease resistance genes were physically mapped to the distal 0.06 fragment of the long arm of wheat chromosome 2D and cosegregated with nematode resistance.

Genome mapping of seed-borne allergens and immunoresponsive proteins in wheat

Wheat is an important staple grain for humankind globally because of its end-use quality and nutritional properties and its adaptability to diverse climates. For a small proportion of the population, specific wheat proteins can trigger adverse immune responses and clinical manifestations such as celiac disease, wheat allergy, baker's asthma, and wheat-dependent exercise-induced anaphylaxis (WDEIA). Establishing the content and distribution of the immunostimulatory regions in wheat has been hampered by the complexity of the wheat genome and the lack of complete genome sequence information. We provide novel insights into the wheat grain proteins based on a comprehensive analysis and annotation of the wheat prolamin Pfam clan grain proteins and other non-prolamin allergens implicated in these disorders using the new International Wheat Genome Sequencing Consortium bread wheat reference genome sequence, RefSeq v1.0. Celiac disease and WDEIA genes are primarily expressed in the starchy endosperm and show wide variation in protein- and transcript-level expression in response to temperature stress. Nonspecific lipid transfer proteins and α-amylase trypsin inhibitor gene families, implicated in baker's asthma, are primarily expressed in the aleurone layer and transfer cells of grains and are more sensitive to cold temperature. The study establishes a new reference map for immunostimulatory wheat proteins and provides a fresh basis for selecting wheat lines and developing diagnostics for products with more favorable consumer attributes.

The molecular–genetic analysis of Triticum tauschii, the D-genome donor to hexaploid wheat

DNA from Triticum tauschii (the D-genome donor to hexaploid wheat, Triticum aestivum) has been cloned using the restriction endonuclease PstI to generate fragments for insertion into the plasmid pBR322 or pUC118. A total of 143 clones were studied and demonstrated to contain one of the following sequence types: (i) a D-genome amplified repetitive sequence, (ii) polymorphic repetitive sequences ("fingerprint-type" sequences), (iii) polymorphic low to moderately repetitive sequences (PLR sequences), (iv) polymorphic low copy number sequences (PL sequences), and (v) invariant sequences. The D-genome amplified sequence was found to be located on all seven chromosomes. A genetic linkage map using PL and PLR sequences was produced by analysing F2 segregating progeny from crosses between two different taxa of T. tauschii. In addition to using DNA clones from T. tauschii, clones from other laboratories containing either anonymous sequences or genes coding for known products (e.g., 7S globulin, dehydrin, germin, storage protein) were used in the genetic linkage map. Multiple locations were mapped for the PLR sequences and were often clustered on single chromosomes. The restriction fragment length polymorphism markers and isozymes analysed were generally distributed over all the linkage groups that were identified and, when used in conjunction with published markers, provided the basis for a genetic map of T. tauschii.Key words: D genome, genetic linkage, restriction fragment length polymorphism, isozymes, chromosomal location.