Dynamics and differential proliferation of transposable elements during the evolution of the B and A genomes of wheat

Transposable elements (TEs) constitute >80% of the wheat genome but their dynamics and contribution to size variation and evolution of wheat genomes (Triticum and Aegilops species) remain unexplored. In this study, 10 genomic regions have been sequenced from wheat chromosome 3B and used to constitute, along with all publicly available genomic sequences of wheat, 1.98 Mb of sequence (from 13 BAC clones) of the wheat B genome and 3.63 Mb of sequence (from 19 BAC clones) of the wheat A genome. Analysis of TE sequence proportions (as percentages), ratios of complete to truncated copies, and estimation of insertion dates of class I retrotransposons showed that specific types of TEs have undergone waves of differential proliferation in the B and A genomes of wheat. While both genomes show similar rates and relatively ancient proliferation periods for the Athila retrotransposons, the Copia retrotransposons proliferated more recently in the A genome whereas Gypsy retrotransposon proliferation is more recent in the B genome. It was possible to estimate for the first time the proliferation periods of the abundant CACTA class II DNA transposons, relative to that of the three main retrotransposon superfamilies. Proliferation of these TEs started prior to and overlapped with that of the Athila retrotransposons in both genomes. However, they also proliferated during the same periods as Gypsy and Copia retrotransposons in the A genome, but not in the B genome. As estimated from their insertion dates and confirmed by PCR-based tracing analysis, the majority of differential proliferation of TEs in B and A genomes of wheat (87 and 83%, respectively), leading to rapid sequence divergence, occurred prior to the allotetraploidization event that brought them together in Triticum turgidum and Triticum aestivum, <0.5 million years ago. More importantly, the allotetraploidization event appears to have neither enhanced nor repressed retrotranspositions. We discuss the apparent proliferation of TEs as resulting from their insertion, removal, and/or combinations of both evolutionary forces.     

Rht-1 and Ppd-D1 associations with height, GA sensitivity, and days to heading in a worldwide bread wheat collection

Reduced height (Rht)-1 and Photoperiod (Ppd) have major effects on the adaptability of bread wheat (Triticum aestivum) to specific environments. Ppd–D1a is a photoperiod insensitive allele that reduces time to flowering. The gibberellin (GA) insensitive alleles Rht–B1b and Rht–D1b shorten plant stature and were important components of the ‘green revolution’. Two additional Rht–B1 alleles were recently identified that contain a 160 or 197 bp insertion upstream of the coding region and may affect plant height or GA sensitivity Wilhelm et al. (Theor Appl Gen doi:10.1007/s00122-013-2088-7, 2013b). We determined the frequency of the five alleles in a worldwide core collection of 372 wheat accessions (372CC) and estimated their effects on height, days to heading, and GA sensitivity when the collection was grown in pots outdoors or in the glasshouse. This revealed that each allele was widespread geographically with frequencies ranging from 0.12 to 0.25. Ppd-D1a was associated with significant (p ≤ 0.05) reductions in days to heading and height relative to photoperiod sensitive Ppd-D1b. Relative to wild type, Rht-B1b and Rht-D1b each resulted in significant reductions in height (approximately 30 %) and GA sensitivity. The 160 and 197 bp alleles were associated with significant height reductions of 18 and 12 %, respectively, and with non-significant reductions in GA sensitivity relative to wild type. Two statistical methods were developed and used to estimate GA sensitivity of the 372CC accessions, but novel GA insensitive alleles were not identified. Further characterization of the Rht-B1 insertion alleles is required, but our results suggest these may enable fine adjustments in plant height.     

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

Key message

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement.

Abstract

To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1U-6U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.     

Isolation of a gibberellin-insensitive dwarfing gene, Rht-B1e, and development of an allele-specific PCR marker

Gibberellins (GAs) are important phytohormones in plants. GAs promote plant growth by inducing the degradation of DELLA proteins, which serve as GA signal repressors. The semi-dwarfing genes Rht-B1b and Rht-D1b, derived from the Japanese variety Norin 10, are gain-of-function mutant alleles of the reduced height-1 genes (Rht-B1 and Rht-D1) encoding wheat DELLA proteins. Wheat varieties carrying these Rht alleles are shorter and insensitive to the GA response. At the Rht-B1 loci, an alternative GA-insensitive dwarfing gene, Rht-B1e, was found in the Russian mutant of Bezostaya1, or Krasznodari 1, by breeders, but its molecular mechanism for causing dwarfism remains unknown. In this study, the Rht-B1e allele was isolated using homology-based cloning. Sequence comparison between Rht-B1e and the wild-type Rht-B1a revealed an A-to-T substitution at nucleotide position 181 in Rht-B1e, which introduced a stop codon into the DELLA domain. Alignment of deduced amino acid sequences of Rht-B1e and Rht-B1b showed that the stop codon position in Rht-B1e was earlier than that of Rht-B1b by three amino acid residues, and it was also followed closely by several methionines, which may permit translational re-initiation, as seen in Rht-B1b. Yeast two-hybrid analysis revealed that the predicted Rht-B1e proteins did not interact with the GA receptor GID1 in the presence of GA, suggesting that the stop codon mutation in the DELLA domain is the molecular cause of GA insensitivity and dwarfism conferred by Rht-B1e in wheat. Meanwhile, we developed an allele-specific PCR marker for Rht-B1e, which may facilitate the use of the Rht-B1e dwarfing gene in wheat breeding programs.     

Comparative RFLP mapping of the chlorotoluron resistance gene (Su1) in cultivated wheat (Triticum aestivum) and wild wheat (Triticum dicoccoides)

Chlorotoluron is a selective phenylurea herbicide widely used for broad-leaved and annual grass weed control in cereals. Variation in the response to chlorotoluron (CT) was found in both hexaploid bread wheat (Triticum aestivum L.) and wild tetraploid wheat (Triticum dicoccoides KöRN.). Here, we describe the comparative mapping of the CT resistance gene (Su1) on chromosome 6B in bread and wild wheat using RFLP markers. In bread wheat, mapping was based on 58 F₄ single-seed descent (SSD) plants of the cross between a genotype sensitive to chlorotoluron, ‘Chinese Spring’ (CS), and a resistant derivative, the single chromosome substitution line, CS (‘Cappele-Desprez’ 6B) [CS (CAP6B). In T dicoccoides, mapping was based on 37 F₂ plants obtained from the cross between the CT-susceptible accession B-7 and the resistant accession B-35. Nine RFLP probes spanning the centromere were chosen for mapping. In bread wheat Su1 was found to be linked to α-Amy-1 (9.84 cM) and Xpsr371 (5.2 cM), both on the long arm of 6B, and Nor2 (2.74 cM) on the short arm. In wild wheat the most probable linkage map was Nor2-Xpsr312-Su1-Pgk2, and the genetic distances between the genes were 24.8cM, 5.3cM, and 6.8cM, respectively. These results along with other published map data indicate that the linear order of the genes is similar to that found in T. aestivum. The results of this study also show that the Su1 gene for differential response to chlorotoluron has evolved prior to the domestication of cultivated wheat and not in response to the development and use of chemicals.     

Cell organisation,sulphur metabolism and ion transport-related genes are differentially expressed in Paracoccidioides brasiliensis mycelium and yeast cells

Background

Bread wheat (Triticum aestivum) is an important staple food. However, wheat gluten proteins cause celiac disease (CD) in 0.5 to 1% of the general population. Among these proteins, the α-gliadins contain several peptides that are associated to the disease.

Results

We obtained 230 distinct α-gliadin gene sequences from severaldiploid wheat species representing the ancestral A, B, and D genomes of the hexaploid bread wheat. The large majority of these sequences (87%) contained an internal stop codon. All α-gliadin sequences could be distinguished according to the genome of origin on the basis of sequence similarity, of the average length of the polyglutamine repeats, and of the differences in the presence of four peptides that have been identified as T cell stimulatory epitopes in CD patients through binding to HLA-DQ2/8. By sequence similarity, α-gliadins from the public database of hexaploid T. aestivum could be assigned directly to chromosome 6A, 6B, or 6D. T. monococcum (A genome) sequences, as well as those from chromosome 6A of bread wheat, almost invariably contained epitope glia-α9 and glia-α20, but never the intact epitopes glia-α and glia-α2. A number of sequences from T. speltoides, as well as a number of sequences fromchromosome 6B of bread wheat, did not contain any of the four T cell epitopes screened for. The sequences from T. tauschii (D genome), as well as those from chromosome 6D of bread wheat, were found to contain all of these T cell epitopes in variable combinations per gene. The differences in epitope composition resulted mainly from point mutations. These substitutions appeared to be genome specific.

Conclusion

Our analysis shows that α-gliadin sequences from the three genomes of bread wheat form distinct groups. The four known T cell stimulatory epitopes are distributed non-randomly across the sequences, indicating that the three genomes contribute differently to epitope content. A systematic analysis of all known epitopes in gliadins and glutenins will lead to better understanding of the differences in toxiCity among wheat varieties. On the basis of such insight, breeding strategies can be designed to generate less toxic varieties of wheat which may be tolerated by at least part of the CD patient population.     

Genetic and epigenetic alteration among three homoeologous genes of a class E MADS box gene in hexaploid wheat

Bread wheat (Triticum aestivum) is a hexaploid species with A, B, and D ancestral genomes. Most bread wheat genes are present in the genome as triplicated homoeologous genes (homoeologs) derived from the ancestral species. Here, we report that both genetic and epigenetic alterations have occurred in the homoeologs of a wheat class E MADS box gene. Two class E genes are identified in wheat, wheat SEPALLATA (WSEP) and wheat LEAFY HULL STERILE1 (WLHS1), which are homologs of Os MADS45 and Os MADS1 in rice (Oryza sativa), respectively. The three wheat homoeologs of WSEP showed similar genomic structures and expression profiles. By contrast, the three homoeologs of WLHS1 showed genetic and epigenetic alterations. The A genome WLHS1 homoeolog (WLHS1-A) had a structural alteration that contained a large novel sequence in place of the K domain sequence. A yeast two-hybrid analysis and a transgenic experiment indicated that the WLHS1-A protein had no apparent function. The B and D genome homoeologs, WLHS1-B and WLHS1-D, respectively, had an intact MADS box gene structure, but WLHS1-B was predominantly silenced by cytosine methylation. Consequently, of the three WLHS1 homoeologs, only WLHS1-D functions in hexaploid wheat. This is a situation where three homoeologs are differentially regulated by genetic and epigenetic mechanisms.     

Resistance of bread wheat (Triticum aestivum L.) to preharvest sprouting: An association analysis

A statistical analysis of the data about 1422 bread wheat accessions with estimated preharvest sprouting was carried out. Close associations of preharvest sprouting resistance with the grain color and with resistance to Fusarium head blight were revealed, as well as weak, but statistically significant, associations with the habit, awnedness, and reduced height genes Rht-B1 and Rht-D1 (insensitive to gibberellin GA3). The pedigree analysis showed that the cluster structures of the gene pools of the North American red-grained and white-grained varieties are practically identical. In both groups, varieties that are resistant to preharvest sprouting differ from susceptible ones in the percentage of the contributions of the Crimean and Mediterranean landraces. Resistance is associated with a high contribution by the Crimean landrace and susceptibility is associated with a high contribution by the Mediterranean landrace.     

Wild-type alleles of Rht-B1 and Rht-D1 as independent determinants of thousand-grain weight and kernel number per spike in wheat

The utilization of dwarfing genes Rht-B1b and Rht-D1b in wheat significantly increased grain yield and contributed to the “green revolution”. However, the benefit of Rht-B1b and Rht-D1b in drought environments has been debated. Although quantitative trait loci (QTL) for kernel number per spike (KN) and thousand-grain weight (TGW) have been found to be associated with Rht-B1 and Rht-D1, the confounding effect of environmental variation has made a direct association difficult to find. In this study, we used a doubled haploid population (225 lines) of Westonia × Kauz, in which both Rht-B1b (Kauz) and Rht-D1b (Westonia) segregated. The purpose of the study was to determine the interaction of Rht-B1 and Rht-D1 with grain yield components, namely KN and TGW, and to investigate genotype-by-environment interactions in glasshouse and field trials conducted in 2010 and 2011 in Western Australia. A genetic map of 1,156 loci was constructed using 195 microsatellite markers, two gene-based markers for Rht-B1 and Rht-D1, and 959 single nucleotide polymorphisms. The major QTL for TGW and KN were strongly linked to Rht-B1 and Rht-D1 loci and the positive effects were associated with the wild-type alleles, Rht-B1a and Rht-D1a. The major QTL of TGW were on chromosome 2D and 4B. The significant genetic effects (14.6–22.9 %) of TGW indicated that marker-assisted selection for TGW is possible, and markers gwm192a (206 bp) or gwm192b (236 bp) can be used as indicators of high TGW. For KN, one major QTL was detected on chromosome 4D in the analysis across three environments. The association of the wild-type alleles Rht-B1a and Rht-D1a in drought environments is discussed.     

An Efficient Approach for the Development of Locus Specific Primers in Bread Wheat (Triticum aestivum L.) and Its Application to Re-Sequencing of Genes Involved in Frost Tolerance

Recent declines in costs accelerated sequencing of many species with large genomes, including hexaploid wheat (Triticum aestivum L.). Although the draft sequence of bread wheat is known, it is still one of the major challenges to developlocus specific primers suitable to be used in marker assisted selection procedures, due to the high homology of the three genomes. In this study we describe an efficient approach for the development of locus specific primers comprising four steps, i.e. (i) identification of genomic and coding sequences (CDS) of candidate genes, (ii) intron- and exon-structure reconstruction, (iii) identification of wheat A, B and D sub-genome sequences and primer development based on sequence differences between the three sub-genomes, and (iv); testing of primers for functionality, correct size and localisation. This approach was applied to single, low and high copy genes involved in frost tolerance in wheat. In summary for 27 of these genes for which sequences were derived from Triticum aestivum, Triticum monococcum and Hordeum vulgare, a set of 119 primer pairs was developed and after testing on Nulli-tetrasomic (NT) lines, a set of 65 primer pairs (54.6%), corresponding to 19 candidate genes, turned out to be specific. Out of these a set of 35 fragments was selected for validation via Sanger''s amplicon re-sequencing. All fragments, with the exception of one, could be assigned to the original reference sequence. The approach presented here showed a much higher specificity in primer development in comparison to techniques used so far in bread wheat and can be applied to other polyploid species with a known draft sequence.     

Association mapping for Fusarium head blight resistance in European soft winter wheat

Resistance to Fusarium head blight (FHB) is of great importance in wheat breeding programs in the northern hemisphere. In Europe, breeders prefer adapted germplasm as resistance donor because of high grain yield and quality demands. Our objective was to identify chromosomal regions affecting FHB resistance among 455 European soft winter wheat (Triticum aestivum L.) lines using a genome-wide association mapping approach and to analyze the importance of epistatic interactions. All entries were evaluated for FHB resistance by inoculation in two environments and several ratings. Wheat was genotyped by 115 simple sequence repeat markers randomly distributed across the genome and two allele-specific markers for Rht-B1 and Rht-D1 genes. The genome-wide scan revealed nine significant (P < 0.05) marker–phenotype associations on seven chromosomes including dwarfing gene Rht-D1. Using a Bonferroni–Holm correction, three significant associations remained on chromosomes 1B, 1D, and 2D. The proportion of the genotypic variance explained simultaneously by individual markers was 36% and increased to 50% when two digenic epistatic interactions were considered, one of them associated with Rht-B1. In conclusion, new genomic regions on chromosomes 1D and 3A could be found for FHB resistance in European wheat and the effect of epistatic interactions was substantial.     

Molecular characterization of field resistance to Fusarium head blight in two US soft red winter wheat cultivars

In the soft red winter wheat (Triticum aestivum L.) regions of the US, Fusarium head blight (FHB, caused by Fusarium spp.) resistance derived from locally adapted germplasm has been used predominantly. Two soft red winter wheat cultivars, Massey and Ernie, have moderate resistance to FHB. Mapping populations derived from Becker/Massey (B/M) and Ernie/MO 94-317 (E/MO) were evaluated for FHB resistance and other traits in multiple environments. Eight QTL in B/M and five QTL in E/MO were associated with FHB variables including incidence, severity (SEV), index (IND), Fusarium damaged kernels (FDK), deoxynivalenol (DON), and morphological traits flowering time and plant height. Four QTL were common to both populations. Three of them were located at or near known genes: Ppd-D1 on chromosome 2DS, Rht-B1 on 4BS, and Rht-D1 on 4DS. Alleles for dwarf plant height (Rht-B1b and Rht-D1b) and photoperiod insensitivity (Ppd-D1a) had pleiotropic effects in reducing height and increasing FHB susceptibility. The other QTL detected for FHB variables were on 3BL in both populations, 1AS, 1DS, 2BL, and 4DL in B/M, and 5AL (B1) and 6AL in E/MO. The additive effects of FHB variables ranged from 0.4 mg kg^?1 of DON to 6.2 % for greenhouse (GH) SEV in B/M and ranged from 0.3 mg kg^?1 of DON to 8.3 % for GH SEV in E/MO. The 4DS QTL had epistasis with Ppd-D1, Qdon.umc-6AL, and Qht.umc-4BS, and additive × additive × environment interactions with the 4BS QTL for SEV, IND, and FDK in E/MO. Marker-assisted selection might be used to enhance FHB resistance through selection of favorable alleles of significant QTL, taking into account genotypes at Rht-B1b, Rht-D1a and Ppd-D1a.     

Molecular characterization and genomic organization of low molecular weight glutenin subunit genes at the Glu-3 loci in hexaploid wheat (Triticum aestivum L.)

In this study, we report on the molecular characterization and genomic organization of the low molecular weight glutenin subunit (LMW-GS) gene family in hexaploid wheat (Triticum aestivum L.). Eighty-two positive BAC clones were identified to contain LMW-GS genes from the hexaploid wheat ‘Glenlea’ BAC library via filter hybridization and PCR validation. Twelve unique LMW glutenin genes and seven pseudogenes were isolated from these positive BAC clones by primer-template mismatch PCR and subsequent primer walking using hemi-nested touchdown PCR. These genes were sequenced and each consisted of a single-open reading frame (ORF) and untranslated 5′ and 3′ flanking regions. All 12 LMW glutenin subunits contained eight cysteine residues. The LMW-m-type subunits are the most abundant in hexaploid wheat. Of the 12 LMW-GS, 1, 2 and 9 are i-type, s-type and m-type, respectively. The phylogenetic analysis suggested that the LMW-i type gene showed greater differences to LMW-s and LMW-m-type genes, which, in turn, were more closely related to one another. On the basis of their N-terminal sequences, they were classified into nine groups. Fingerprinting of the 82 BAC clones indicated 30 BAC clones assembled into eight contigs, while the remaining clones were singletons. BAC end sequencing of the 82 clones revealed that long terminal repeat (LTR) retrotransposons were abundant in the Glu-3 regions. The average physical distance between two adjacent LMW-GS genes was estimated to be 81 kb. Most of LMW-GS genes are located in the d-genome, suggesting that the Glu-D3 locus is much larger than the Glu-B3 locus and Glu-A3 locus. Alignments of sequences indicated that the same type (starting with the same N-terminal sequence) LMW-GS genes were highly conserved in the homologous genomes between hexaploid wheat and its donors such as durum wheat and T. tauschii. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.