Comparative sequence analysis of the region harboring the hardness locus in barley and its colinear region in rice

The gene conferring susceptibility to spot blotch caused by Cochliobolus sativus is located at the Mla locus in barley cultivar Bowman

We identified, fine mapped, and physically anchored a dominant spot blotch susceptibility gene Scs6 to a 125 kb genomic region containing the Mla locus on barley chromosome 1H. Spot blotch caused by Cochliobolus sativus is an important disease of barley, but the molecular mechanisms underlying resistance and susceptibility to the disease are not well understood. In this study, we identified and mapped a gene conferring susceptibility to spot blotch caused by the pathotype 2 isolate (ND90Pr) of C. sativus in barley cultivar Bowman. Genetic analysis of F₁ and F₂ progeny as well as F₃ families from a cross between Bowman and ND 5883 indicated that a single dominant gene (designated as Scs6) conferred spot blotch susceptibility in Bowman. Using a doubled haploid (DH) population derived from a cross between Calicuchima-sib (resistant) and Bowman-BC (susceptible), we confirmed that Scs6, contributed by Bowman-BC, was localized at the same locus as the previously identified spot blotch resistance allele Rcs6, which was contributed by Calicuchima-sib and mapped on the short arm of chromosome 1H. Using a genome-wide putative linear gene index of barley (Genome Zipper), 13 cleaved amplified polymorphism markers were developed from 11 flcDNA and two EST sequences and mapped to the Scs6/Rcs6 region on a linkage map constructed with the DH population. Further fine mapping with markers developed from barley genome sequences and F₂ recombinants derived from Bowman × ND 5883 and Bowman × ND B112 crosses delimited Scs6 in a 125 kb genomic interval harboring the Mla locus on the reference genome of barley cv. Morex. This study provides a foundational step for further cloning of Scs6 using a map-based approach.

Conservation and divergence of small RNA pathways and microRNAs in land plants

Background

As key regulators of gene expression in eukaryotes, small RNAs have been characterized in many seed plants, and pathways for their biogenesis, degradation, and action have been defined in model angiosperms. However, both small RNAs themselves and small RNA pathways are not well characterized in other land plants such as lycophytes and ferns, preventing a comprehensive evolutionary perspective on small RNAs in land plants.

Results

Using 25 representatives from major lineages of lycophytes and ferns, most of which lack sequenced genomes, we characterized small RNAs and small RNA pathways in these plants. We identified homologs of DICER-LIKE (DCL), ARGONAUTE (AGO), and other genes involved in small RNA pathways, predicted over 2600 conserved microRNA (miRNA) candidates, and performed phylogenetic analyses on small RNA pathways as well as miRNAs. Pathways underlying miRNA biogenesis, degradation, and activity were established in the common ancestor of land plants, but the 24-nucleotide siRNA pathway that guides DNA methylation is incomplete in sister species of seed plants, especially lycophytes. We show that the functional diversification of key gene families such as DCL and AGO as observed in angiosperms occurred early in land plants followed by parallel expansion of the AGO family in ferns and angiosperms. We uncovered a conserved AGO subfamily absent in angiosperms.

Conclusions

Our phylogenetic analyses of miRNAs in bryophytes, lycophytes, ferns, and angiosperms refine the time-of-origin for conserved miRNA families as well as small RNA machinery in land plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1291-2) contains supplementary material, which is available to authorized users.

Assessing the Barley Genome Zipper and Genomic Resources for Breeding Purposes

The aim of this study was to estimate the accuracy and convergence of newly developed barley (Hordeum vulgare L.) genomic resources, primarily genome zipper (GZ) and population sequencing (POPSEQ), at the genome-wide level and to assess their usefulness in applied barley breeding by analyzing seven known loci. Comparison of barley GZ and POPSEQ maps to a newly developed consensus genetic map constructed with data from 13 individual linkage maps yielded an accuracy of 97.8% (GZ) and 99.3% (POPSEQ), respectively, regarding the chromosome assignment. The percentage of agreement in marker position indicates that on average only 3.7% GZ and 0.7% POPSEQ positions are not in accordance with their centimorgan coordinates in the consensus map. The fine-scale comparison involved seven genetic regions on chromosomes 1H, 2H, 4H, 6H, and 7H, harboring major genes and quantitative trait loci (QTL) for disease resistance. In total, 179 GZ loci were analyzed and 64 polymorphic markers were developed. Entirely, 89.1% of these were allocated within the targeted intervals and 84.2% followed the predicted order. Forty-four markers showed a match to a POPSEQ-anchored contig, the percentage of collinearity being 93.2%, on average. Forty-four markers allowed the identification of twenty-five fingerprinted contigs (FPCs) and a more clear delimitation of the physical regions containing the traits of interest. Our results demonstrate that an increase in marker density of barley maps by using new genomic data significantly improves the accuracy of GZ. In addition, the combination of different barley genomic resources can be considered as a powerful tool to accelerate barley breeding.

The barley Frost resistance-H2 locus

Frost resistance-H2 (Fr-H2) is a major QTL affecting freezing tolerance in barley, yet its molecular basis is still not clearly understood. To gain a better insight into the structural characterization of the locus, a high-resolution linkage map developed from the Nure × Tremois cross was initially implemented to map 13 loci which divided the 0.602 cM total genetic distance into ten recombination segments. A PCR-based screening was then applied to identify positive bacterial artificial chromosome (BAC) clones from two genomic libraries of the reference genotype Morex. Twenty-six overlapping BACs from the integrated physical-genetic map were 454 sequenced. Reads assembled in contigs were subsequently ordered, aligned and manually curated in 42 scaffolds. In a total of 1.47 Mbp, 58 protein-coding sequences were identified, 33 of which classified according to similarity with sequences in public databases. As three complete barley C-repeat Binding Factors (HvCBF) genes were newly identified, the locus contained13 full-length HvCBFs, four Related to AP2 Triticeae (RAPT) genes, and at least five CBF pseudogenes. The final overall assembly of Fr-H2 includes more than 90 % of target region: all genes were identified along the locus, and a general survey of Repetitive Elements obtained. We believe that this gold-standard sequence for the Morex Fr-H2 will be a useful genomic tool for structural and evolutionary comparisons with Fr-H2 in winter-hardy cultivars along with Fr-2 of other Triticeae crops.

The detection of QTLs in barley associated with endosperm hardness, grain density, grain size and malting quality using rapid phenotyping tools

Using a barley mapping population, 'Vlamingh' × 'Buloke' (V × B), whole grain analyses were undertaken for physical seed traits and malting quality. Grain density and size were predicted by digital image analysis (DIA), while malt extract and protein content were predicted using near infrared (NIR) analysis. Validation of DIA and NIR algorithms confirmed that data for QTL analysis was highly correlated (R (2) > 0.82), with high RPD values (the ratio of the standard error of prediction to the standard deviation, 2.31-9.06). Endosperm hardness was measured on this mapping population using the single kernel characterisation system. Grain density and endosperm hardness were significantly inter-correlated in all three environments (r > 0.22, P < 0.001); however, other grain components were found to interact with the traits. QTL for these traits were also found on different genomic regions, for example, grain density QTLs were found on chromosomes 2H and 6H, whereas endosperm hardness QTLs were found on 1H, 5H, and 7H. In this study, the majority of the genomic regions associated with grain texture were also coincident with QTLs for grain size, yield, flowering date and/or plant development genes. This study highlights the complexity of genomic regions associated with the variation of endosperm hardness and grain density, and their relationships with grain size traits, agronomic-related traits, and plant development loci.

Identification and distribution of Puroindoline b-2 variant gene homologs in Hordeum

The barley hordoindoline genes (Hina and Hinb) are homologous to the wheat puroindoline genes (Pina and Pinb). These genes are involved in grain hardness, which is an important quality for barley processing. We identified novel variants of Hina and Hinb in 10 wild Hordeum species (H. bogdanii, H. brachyantherum, H. bulbosum, H. chilense, H. comosum, H. marinum, H. murinum, H. patagonicum, H. pusillum, and H. roshevitzii) covering all Hordeum genomes and preliminarily named them Hinc. These nucleotide sequences were highly similar to those of Puroindoline b-2 variant genes (Pinb-2v) and were located on chromosome 7I in H. chilense. The Hinc genes in H. bogdanii, H. bulbosum, H. patagonicum, and H. roshevitzii were pseudogenes possessing in-frame stop codons. We also found a partial Hinc sequence in H. murinum. This gene was not found in cultivated barley and H. vulgare subsp. spontaneum. The phylogenetic tree of Gsp-1, Hin, and Pin genes demonstrates that Hinc and Pinb-2v genes formed one cluster. Therefore, we considered that Hinc and Pinb-2v genes shared a common ancestral gene and were homologous to each other. We also studied the evolutional process of Gsp-1, Hin, and Pin genes. Our results suggested that Gsp-1 might be the most closely related to a putative ancestral gene on Ha locus.

Towards positional isolation of three quantitative trait loci conferring resistance to powdery mildew in two Spanish barley landraces

Three quantitative trait loci (QTL) conferring broad spectrum resistance to powdery mildew, caused by the fungus Blumeria graminis f. sp. hordei, were previously identified on chromosomes 7HS, 7HL and 6HL in the Spanish barley landrace-derived lines SBCC097 and SBCC145. In the present work, a genome-wide putative linear gene index of barley (Genome Zipper) and the first draft of the physical, genetic and functional sequence of the barley genome were used to go one step further in the shortening and explicit demarcation on the barley genome of these regions conferring resistance to powdery mildew as well as in the identification of candidate genes. First, a comparative analysis of the target regions to the barley Genome Zippers of chromosomes 7H and 6H allowed the development of 25 new gene-based molecular markers, which slightly better delimit the QTL intervals. These new markers provided the framework for anchoring of genetic and physical maps, figuring out the outline of the barley genome at the target regions in SBCC097 and SBCC145. The outermost flanking markers of QTLs on 7HS, 7HL and 6HL defined a physical area of 4 Mb, 3.7 Mb and 3.2 Mb, respectively. In total, 21, 10 and 16 genes on 7HS, 7HL and 6HL, respectively, could be interpreted as potential candidates to explain the resistance to powdery mildew, as they encode proteins of related functions with respect to the known pathogen defense-related processes. The majority of these were annotated as belonging to the NBS-LRR class or protein kinase family.

New algorithm improves fine structure of the barley consensus SNP map

BACKGROUND

The need to integrate information from multiple linkage maps is a long-standing problem in genetics. One way to visualize the complex ordinal relationships is with a directed graph, where each vertex in the graph is a bin of markers. When there are no ordering conflicts between the linkage maps, the result is a directed acyclic graph, or DAG, which can then be linearized to produce a consensus map.

RESULTS

New algorithms for the simplification and linearization of consensus graphs have been implemented as a package for the R computing environment called DAGGER. The simplified consensus graphs produced by DAGGER exactly capture the ordinal relationships present in a series of linkage maps. Using either linear or quadratic programming, DAGGER generates a consensus map with minimum error relative to the linkage maps while remaining ordinally consistent with them. Both linearization methods produce consensus maps that are compressed relative to the mean of the linkage maps. After rescaling, however, the consensus maps had higher accuracy (and higher marker density) than the individual linkage maps in genetic simulations. When applied to four barley linkage maps genotyped at nearly 3000 SNP markers, DAGGER produced a consensus map with improved fine structure compared to the existing barley consensus SNP map. The root-mean-squared error between the linkage maps and the DAGGER map was 0.82 cM per marker interval compared to 2.28 cM for the existing consensus map. Examination of the barley hardness locus at the 5HS telomere, for which there is a physical map, confirmed that the DAGGER output was more accurate for fine structure analysis.

CONCLUSIONS

The R package DAGGER is an effective, freely available resource for integrating the information from a set of consistent linkage maps.

Chromosomal loci associated with endosperm hardness in a malting barley cross

A breeding objective for the malting barley industry is to produce lines with softer, plumper grain containing moderate protein content (9-12%) as they are more likely to imbibe water readily and contain more starch per grain, which in turn produces higher levels of malt extract. In a malting barley mapping population, 'Arapiles' × 'Franklin', the most significant and robust quantitative trait locus (QTL) for endosperm hardness was observed on the short arm of chromosome 1H, across three environments over two growing seasons. This accounted for 22.6% (Horsham 2000), 26.8% (Esperance 2001), and 12.0% (Tarranyurk 2001) of the genetic variance and significantly increased endosperm hardness by 2.06-3.03 SKCS hardness units. Interestingly, Arapiles and Franklin do not vary in Ha locus alleles. Therefore, this region, near the centromere on chromosome 1H, may be of great importance when aiming to manipulate endosperm hardness and malting quality. Interestingly, this region, close to the centromere on chromosome 1H, in our study, aligns with the region of the genome that includes the HvCslF9 and the HvGlb1 genes. Potentially, one or both of these genes could be considered to be candidate genes that influence endosperm hardness in the barley grain. Additional QTLs for endosperm hardness were detected on chromosomes 2H, 3H, 6H and 7H, confirming that the hardness trait in barley is complex and multigenic, similar to many malting quality traits of interest.

Genetic analysis and ecological association of Hina genes based on single nucleotide polymorphisms (SNPs) in wild barley, Hordeum spontaneum

Specific primers were designed to amplify the sequences of Hina genes from 121 wild barley (Hordeum spontaneum) accessions belonging to 18 populations from Iran, Israel and Turkey. Forty-nine single nucleotide polymorphisms (SNPs), nine indels, and 26 haplotypes were determined by sequence analysis. The genetic polymorphism (P), genetic diversity (He), and Shannon's information index (I) in the 18 populations were 0.486, 0.181 and 0.269, respectively. Approximately 2/3 genetic variations of Hina genes were presented within populations, while approximately 1/3 genetic variations were observed between populations. Broad gene flow (Nm= 3.31) and low genetic variation (Gst= 0.0702) were detected. However, the genetic differentiation between populations was independent of geographical distances according to the Mantel test (p = 0.478). The result of Spearman rank correlations (r(s)) showed that the genetic indices (P, He and I) of Hina were not significantly correlated with ecological factors. Only eight SNP positions correlated significantly with ecological factors. Of the eight SNP positions that positively correlated with ecological factors, only one SNP (769, T-C) was located in the coding region; however, it was not responsible for the amino acid change.

Tracing the evolution of the floral homeotic B- and C-function genes through genome synteny

The evolution of the floral homeotic genes has been characterized using phylogenetic and functional studies. It is possible to enhance these studies by comparing gene content and order between species to determine the evolutionary history of the regulatory genes. Here, we use a synteny-based approach to trace the evolution of the floral B- and C-function genes that are required for specification of the reproductive organs. Consistent with previous phylogenetic studies, we show that the euAP3–TM6 split occurred after the monocots and dicots diverged. The Arabidopsis TM6 and papaya euAP3 genes are absent from the respective genomes, and we have detected loci from which these genes were lost. These data indicate that either the TM6 or the euAP3 lineage genes can be lost without detriment to flower development. In contrast, PI is essential for male reproductive organ development; yet, contrary to predictions, complex genomic rearrangements have resulted in almost complete breakdown of synteny at the PI locus. In addition to showing the evolution of B-function genes through the prediction of ancestral loci, similar reconstructions reveal the origins of the C-function AG and PLE lineages in dicots, and show the shared ancestry with the monocot C-function genes. During our studies, we found that transposable elements (TEs) present in sequenced Antirrhinum genomic clones limited comparative studies. A pilot survey of the Antirrhinum data revealed that gene-rich regions contain an unusually high degree of TEs of very varied types, which will be an important consideration for future genome sequencing efforts.

A highly conserved gene island of three genes on chromosome 3B of hexaploid wheat: diverse gene function and genomic structure maintained in a tightly linked block

BACKGROUND

The complexity of the wheat genome has resulted from waves of retrotransposable element insertions. Gene deletions and disruptions generated by the fast replacement of repetitive elements in wheat have resulted in disruption of colinearity at a micro (sub-megabase) level among the cereals. In view of genomic changes that are possible within a given time span, conservation of genes between species tends to imply an important functional or regional constraint that does not permit a change in genomic structure. The ctg1034 contig completed in this paper was initially studied because it was assigned to the Sr2 resistance locus region, but detailed mapping studies subsequently assigned it to the long arm of 3B and revealed its unusual features.

RESULTS

BAC shotgun sequencing of the hexaploid wheat (Triticum aestivum cv. Chinese Spring) genome has been used to assemble a group of 15 wheat BACs from the chromosome 3B physical map FPC contig ctg1034 into a 783,553 bp genomic sequence. This ctg1034 sequence was annotated for biological features such as genes and transposable elements. A three-gene island was identified among >80% repetitive DNA sequence. Using bioinformatics analysis there were no observable similarity in their gene functions. The ctg1034 gene island also displayed complete conservation of gene order and orientation with syntenic gene islands found in publicly available genome sequences of Brachypodium distachyon, Oryza sativa, Sorghum bicolor and Zea mays, even though the intergenic space and introns were divergent.

CONCLUSION

We propose that ctg1034 is located within the heterochromatic C-band region of deletion bin 3BL7 based on the identification of heterochromatic tandem repeats and presence of significant matches to chromodomain-containing gypsy LTR retrotransposable elements. We also speculate that this location, among other highly repetitive sequences, may account for the relative stability in gene order and orientation within the gene island.Sequence data from this article have been deposited with the GenBank Data Libraries under accession no. GQ422824.

Promoters of the barley germin-like GER4 gene cluster enable strong transgene expression in response to pathogen attack

Immunity of plants triggered by pathogen-associated molecular patterns (PAMPs) is based on the execution of an evolutionarily conserved defense response that includes the accumulation of pathogenesis-related (PR) proteins as well as multiple other defenses. The most abundant PR transcript of barley (Hordeum vulgare) leaf epidermis attacked by the powdery mildew fungus Blumeria graminis f. sp hordei encodes the germin-like protein GER4, which has superoxide dismutase activity and functions in PAMP-triggered immunity. Here, we show that barley GER4 is encoded by a dense cluster of tandemly duplicated genes (GER4a-h) that underwent several cycles of duplication. The genomic organization of the GER4 locus also provides evidence for repeated gene birth and death cycles. The GER4 promoters contain multiple WRKY factor binding sites (W-boxes) preferentially located in promoter fragments that were exchanged between subfamily members by gene conversion. Mutational analysis of TATA-box proximal W-boxes used GER4c promoter-beta-glucuronidase fusions to reveal their enhancing effects and functional redundancy on pathogen-induced promoter activity. The data suggest enhanced transcript dosage as an evolutionary driving force for the local expansion and functional redundancy of the GER4 locus. In addition, the GER4c promoter provides a tool to study signal transduction of PAMP-triggered immunity and to engineer strictly localized and pathogen-regulated disease resistance in transgenic cereal crops.

The determinants of grain texture in cereals

Kernel hardness is an important agronomic trait that influences end-product properties. In wheat cultivars, this trait is determined by the Puroindoline a (Pina) and Puroindoline b (Pinb) genes, located in the Hardness locus (Ha) on chromosome 5DS of the D genome. Wild type alleles code puroindoline a (PINA) and puroindoline b (PINB) proteins, which form a 15-kDa friabilin present on the surface of water-washed starch granules. Both the proteins are accumulated in the starch endosperm cells and aleurone of the mature kernels. Puroindoline-like genes coding puroindoline-like proteins in the starch endosperm occur in some of the genomes of Triticeae and Aveneae cereals. Orthologs are present in barley, rye and oats. However, some genomes of these diploid and polyploid cereals, like that of Triticum turgidum var. durum (AABB) lack the puroindoline genes, having a very hard kernel texture. The two wild type alleles in opposition (dominant loci) control the soft phenotype. Mutation either in Pina or Pinb or in both leads to a medium-hard or hard kernel texture. The most frequent types of Pin mutations are point mutations within the coding sequence resulting in the substitution of a single amino acid or a null allele. The latter is the result of a frame shift determined by base deletion or insertion or a one-point mutation to the stop codon. The lipid-binding properties of the puroindolines affect not only the dough quality but also the plants' resistance to pathogens. Genetic modification of cereals with Puroindoline genes and/or their promoters enable more detailed functional analyses and the production of plants with the desired characteristics.

Fine mapping, physical mapping and development of diagnostic markers for the Rrs2 scald resistance gene in barley

The Rrs2 gene confers resistance to the fungal pathogen Rhynchosporium secalis which causes leaf scald, a major barley disease. The Rrs2 gene was fine mapped to an interval of 0.08 cM between markers 693M6_6 and P1D23R on the distal end of barley chromosome 7HS using an Atlas (resistant) x Steffi (susceptible) mapping population of 9,179 F(2)-plants. The establishment of a physical map of the Rrs2 locus led to the discovery that Rrs2 is located in an area of suppressed recombination within this mapping population. The analysis of 58 barley genotypes revealed a large linkage block at the Rrs2 locus extending over several hundred kb which is present only in Rrs2 carrying cultivars. Due to the lack of recombination in the mapping population and the presence of a Rrs2-specific linkage block, we assume a local chromosomal rearrangement (alien introgression or inversion) in Rrs2 carrying varieties. The variety analysis led to the discovery of eight SNPs which were diagnostic for the Rrs2 phenotype. Based on these SNPs diagnostic molecular markers (CAPS and pyrosequencing markers) were developed which are highly useful for marker-assisted selection in resistance gene pyramiding programmes for Rhynchosporium secalis resistance in barley.

Fine mapping and marker development for the crossability gene SKr on chromosome 5BS of hexaploid wheat (Triticum aestivum L.)

Most elite wheat varieties cannot be crossed with related species thereby restricting greatly the germplasm that can be used for alien introgression in breeding programs. Inhibition to crossability is controlled genetically and a number of QTL have been identified to date, including the major gene Kr1 on 5BL and SKr, a strong QTL affecting crossability between wheat and rye on chromosome 5BS. In this study, we used a recombinant SSD population originating from a cross between the poorly crossable cultivar Courtot (Ct) and the crossable line MP98 to characterize the major dominant effect of SKr and map the gene at the distal end of the chromosome near the 5B homeologous GSP locus. Colinearity with barley and rice was used to saturate the SKr region with new markers and establish orthologous relationships with a 54-kb region on rice chromosome 12. In total, five markers were mapped within a genetic interval of 0.3 cM and 400 kb of BAC contigs were established on both sides of the gene to lay the foundation for map-based cloning of SKr. Two SSR markers completely linked to SKr were used to evaluate a collection of crossable wheat progenies originating from primary triticale breeding programs. The results confirm the major effect of SKr on crossability and the usefulness of the two markers for the efficient introgression of crossability in elite wheat varieties.

Sixty million years in evolution of soft grain trait in grasses: emergence of the softness locus in the common ancestor of Pooideae and Ehrhartoideae, after their divergence from Panicoideae

Together maize, Sorghum, rice, and wheat grass (Poaceae) species are the most important cereal crops in the world and exhibit different "grain endosperm texture." This trait has been studied extensively in wheat because of its pivotal role in determining quality of products obtained from wheat grain. Grain softness protein-1 and Puroindolines A and B (grain storage proteins), encoded by Ha-like genes: Gsp-1, Pina, and Pinb, of the Hardness (Ha) locus, are the main determinants of the grain softness/hardness trait in wheat. The origin and evolution of grain endosperm texture in grasses was addressed by comparing genomic sequences of the Ha orthologous region of wheat, Brachypodium, rice, and Sorghum. Results show that the Ha-like genes are present in wheat and Brachypodium but are absent from Sorghum bicolor. A truncated remnant of an Ha-like gene is present in rice. Synteny analysis of the genomes of these grass species shows that only one of the paralogous Ha regions, created 70 My by whole-genome duplication, contained Ha-like genes. The comparative genome analysis and evolutionary comparison with genes encoding grain reserve proteins of grasses suggest that an ancestral Ha-like gene emerged, as a new member of the prolamin gene family, in a common ancestor of the Pooideae (Triticeae and Brachypoidieae tribes) and Ehrhartoideae (rice), between 60 and 50 My, after their divergence from Panicoideae (Sorghum). It was subsequently lost in Ehrhartoideae. Recurring duplications, deletions, and/or truncations occurred independently and appear to characterize Ha-like gene evolution in the grass species. The Ha-like genes gained a new function in Triticeae, such as wheat, underlying the soft grain phenotype. Loss of these genes in some wheat species leads, in turn, to hard endosperm seeds.

Low-pass shotgun sequencing of the barley genome facilitates rapid identification of genes, conserved non-coding sequences and novel repeats

Background

Barley has one of the largest and most complex genomes of all economically important food crops. The rise of new short read sequencing technologies such as Illumina/Solexa permits such large genomes to be effectively sampled at relatively low cost. Based on the corresponding sequence reads a Mathematically Defined Repeat (MDR) index can be generated to map repetitive regions in genomic sequences.

Results

We have generated 574 Mbp of Illumina/Solexa sequences from barley total genomic DNA, representing about 10% of a genome equivalent. From these sequences we generated an MDR index which was then used to identify and mark repetitive regions in the barley genome. Comparison of the MDR plots with expert repeat annotation drawing on the information already available for known repetitive elements revealed a significant correspondence between the two methods. MDR-based annotation allowed for the identification of dozens of novel repeat sequences, though, which were not recognised by hand-annotation. The MDR data was also used to identify gene-containing regions by masking of repetitive sequences in eight de-novo sequenced bacterial artificial chromosome (BAC) clones. For half of the identified candidate gene islands indeed gene sequences could be identified. MDR data were only of limited use, when mapped on genomic sequences from the closely related species Triticum monococcum as only a fraction of the repetitive sequences was recognised.

Conclusion

An MDR index for barley, which was obtained by whole-genome Illumina/Solexa sequencing, proved as efficient in repeat identification as manual expert annotation. Circumventing the labour-intensive step of producing a specific repeat library for expert annotation, an MDR index provides an elegant and efficient resource for the identification of repetitive and low-copy (i.e. potentially gene-containing sequences) regions in uncharacterised genomic sequences. The restriction that a particular MDR index can not be used across species is outweighed by the low costs of Illumina/Solexa sequencing which makes any chosen genome accessible for whole-genome sequence sampling.

Structure-function analysis of the barley genome: the gene-rich region of chromosome 2HL

A major gene-rich region on the end of the long arm of Triticeae group 2 chromosomes exhibits high recombination frequencies, making it an attractive region for positional cloning. Traits known to be controlled by this region include chasmogamy/cleistogamy, frost tolerance at flowering, grain yield, head architecture, and resistance to Fusarium head blight and rusts. To assist these cloning efforts, we constructed detailed genetic maps of barley chromosome 2H, including 61 polymerase chain reaction markers. Colinearity with rice occurred in eight distinct blocks, including five blocks in the terminal gene-rich region. Alignment of rice sequences from the junctions of colinear chromosome segments provided no evidence for the involvement of long (>2.5 kb) inverted repeats in generating inversions. However, reuse of some junction sequences in two or three separate evolutionary breakage/fusion events was implicated, suggesting the presence of fragile sites. Sequencing across 91 gene fragments totaling 107 kb from four barley genotypes revealed the highest single nucleotide substitution and insertion-deletion polymorphism levels in the terminal regions of the chromosome arms. The maps will assist in the isolation of genes from the chromosome 2L gene-rich region in barley and wheat by providing markers and accelerating the identification of the corresponding points in the rice genome sequence.

The 172-kb genomic DNA region of the O. rufipogon yld1.1 locus: comparative sequence analysis with O. sativa ssp. japonica and O. sativa ssp. indica

Common wild rice (Oryza rufipogon) plays an important role by contributing to modern rice breeding. In this paper, we report the sequence and analysis of a 172-kb genomic DNA region of wild rice around the RM5 locus, which is associated with the yield QTL yld1.1. Comparative sequence analysis between orthologous RM5 regions from Oryza sativa ssp. japonica, O. sativa ssp. indica and O. rufipogon revealed a high level of conserved synteny in the content, homology, structure, orientation, and physical distance of all 14 predicted genes. Twelve of the putative genes were supported by matches to proteins with known function, whereas two were predicted by homology to rice and other plant expressed sequence tags or complementary DNAs. The remarkably high level of conservation found in coding, intronic and intergenic regions may indicate high evolutionary selection on the RM5 region. Although our analysis has not defined which gene(s) determine the yld1.1 phenotype, allelic variation and the insertion of transposable elements, among other nucleotide changes, represent potential variation responsible for the yield QTL. However, as suggested previously, two putative receptor-like protein kinase genes remain the key suspects for yld1.1.

Contrasted microcolinearity and gene evolution within a homoeologous region of wheat and barley species

We study here the evolution of genes located in the same physical locus using the recently sequenced Ha locus in seven wheat genomes in diploid, tetraploid, and hexaploid species and compared them with barley and rice orthologous regions. We investigated both the conservation of microcolinearity and the molecular evolution of genes, including coding and noncoding sequences. Microcolinearity is restricted to two groups of genes (Unknown gene-2, VAMP, BGGP, Gsp-1, and Unknown gene-8 surrounded by several copies of ATPase), almost conserved in rice and barley, but in a different relative position. Highly conserved genes between wheat and rice run along with genes harboring different copy numbers and highly variable sequences between close wheat genomes. The coding sequence evolution appeared to be submitted to heterogeneous selective pressure and intronic sequences analysis revealed that the molecular clock hypothesis is violated in most cases.

Molecular genetics of puroindolines and related genes: allelic diversity in wheat and other grasses

The hardness or texture of cereal grains is a primary determinant of their technological and processing quality. Among members of the Triticeae, most notably wheat, much of the variation in texture is controlled by a single locus comprised of the Puroindoline a, Puroindoline b and Grain Softness Protein-1 (Gsp-1) genes. Puroindolines confer the three major texture classes of soft and hard common wheat and the very hard durum wheat. The protein products of these genes interact with lipids and are associated with the surface of isolated starch (as a protein fraction known as 'friabilin'). During the past ten years a great diversity of alleles of both Puroindoline genes have been discovered and significant advances made in understanding the relationship between the gene presence/absence, sequence polymorphism and texture of cereal grains. Efforts have also focussed on Puroindoline and Gsp-1 genes in diploid progenitors, other Triticeae grasses and synthetic wheats in order to understand the evolution of this gene family and find potentially useful variants. The puroindoline homologues in other cereals such as rye and barley are also receiving attention. This work summarises new developments in molecular genetics of puroindolines in wheat and related Triticeae grasses, and the related genes in other cereals.

Gene expression quantitative trait locus analysis of 16 000 barley genes reveals a complex pattern of genome-wide transcriptional regulation

Transcript abundance from cRNA hybridizations to Affymetrix microarrays can be used for simultaneous marker development and genome-wide gene expression quantitative trait locus (eQTL) analysis of crops. We have previously shown that it is easily possible to use Affymetrix expression arrays to profile individuals from a segregating population to accurately identify robust polymorphic molecular genetic markers. We applied the method to identify more than 2000 genetic polymorphisms (transcript derived markers, TDMs) from an experiment involving two commercial varieties of barley (Hordeum vulgare; Steptoe and Morex) and their doubled-haploid progeny. With this set of TDMs, we constructed a genetic map and used it for the genome-wide eQTL analysis of about 16 000 genes in a relatively large population (n = 139). We identified 23 738 significant eQTLs at a genome-wide significance (P </= 0.05), affecting the expression of 12 987 genes. Over a third of these genes with expression variation have only one identified eQTL while the rest have two to six. A large proportion of the quantitatively controlled transcripts appear to be regulated by both cis and trans effects. More than half of the quantitatively controlled transcripts appear to be primarily regulated by cis eQTLs in this population. We show that although there appear to be eQTL hotspots many of these are in chromosomal regions of low recombination, such as genetic centromeres, and so have a high gene density per centimorgan. Some chromosomal regions have a significant excess of eQTL over the number expected from gene density, and many of these are biased towards eQTL for which the allele from one particular parent increases the expression level.

TEnest: automated chronological annotation and visualization of nested plant transposable elements

Organisms with a high density of transposable elements (TEs) exhibit nesting, with subsequent repeats found inside previously inserted elements. Nesting splits the sequence structure of TEs and makes annotation of repetitive areas challenging. We present TEnest, a repeat identification and display tool made specifically for highly repetitive genomes. TEnest identifies repetitive sequences and reconstructs separated sections to provide full-length repeats and, for long-terminal repeat (LTR) retrotransposons, calculates age since insertion based on LTR divergence. TEnest provides a chronological insertion display to give an accurate visual representation of TE integration history showing timeline, location, and families of each TE identified, thus creating a framework from which evolutionary comparisons can be made among various regions of the genome. A database of repeats has been developed for maize (Zea mays), rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare) to illustrate the potential of TEnest software. All currently finished maize bacterial artificial chromosomes totaling 29.3 Mb were analyzed with TEnest to provide a characterization of the repeat insertions. Sixty-seven percent of the maize genome was found to be made up of TEs; of these, 95% are LTR retrotransposons. The rate of solo LTR formation is shown to be dissimilar across retrotransposon families. Phylogenetic analysis of TE families reveals specific events of extreme TE proliferation, which may explain the high quantities of certain TE families found throughout the maize genome. The TEnest software package is available for use on PlantGDB under the tools section (http://www.plantgdb.org/prj/TE_nest/TE_nest.html); the source code is available from (http://wiselab.org).

Fine mapping and targeted SNP survey using rice-wheat gene colinearity in the region of the Bo1 boron toxicity tolerance locus of bread wheat

Toxicity due to high levels of soil boron (B) represents a significant limitation to cereal production in some regions, and the Bo1 gene provides a major source of B toxicity tolerance in bread wheat (Triticum aestivum L.). A novel approach was used to develop primers to amplify and sequence gene fragments specifically from the Bo1 region of the hexaploid wheat genome. Single-nucleotide polymorphisms (SNPs) identified were then used to generate markers close to Bo1 on the distal end of chromosome 7BL. In the 16 gene fragments totaling 19.6 kb, SNPs were observed between the two cultivars Cranbrook and Halberd at a low frequency (one every 613 bp). Furthermore, SNPs were distributed unevenly, being limited to only two genes. In contrast, RFLP provided a much greater number of genetic markers, with every tested gene identifying polymorphism. Bo1 previously known only as a QTL was located as a discrete Mendelian locus. In total, 28 new RFLP, PCR and SSR markers were added to the existing map. The 1.8 cM Bo1 interval of wheat corresponds to a 227 kb section of rice chromosome 6L encoding 21 predicted proteins with no homology to any known B transporters. The co-dominant PCR marker AWW5L7 co-segregated with Bo1 and was highly predictive of B tolerance status within a set of 94 Australian bread wheat cultivars and breeding lines. The markers and rice colinearity described here represent tools that will assist B tolerance breeding and the positional cloning of Bo1.

Analysis of the barley chromosome 2 region containing the six-rowed spike gene vrs1 reveals a breakdown of rice-barley micro collinearity by a transposition

In cultivated barley (Hordeum vulgare ssp. vulgare), six-rowed spikes produce three times as many seeds per spike as do two-rowed spikes. The determinant of this trait is the Mendelian gene vrs1, located on chromosome 2H, which is syntenous with rice (Oryza sativa) chromosomes 4 and 7. We exploited barley-rice micro-synteny to increase marker density in the vrs1 region as a prelude to its map-based cloning. The rice genomic sequence, covering a 980 kb contig, identified barley ESTs linked to vrs1. A high level of conservation of gene sequence was obtained between barley chromosome 2H and rice chromosome 4. A total of 22 EST-based STS markers were placed within the target region, and the linear order of these markers in barley and rice was identical. The genetic window containing vrs1 was narrowed from 0.5 to 0.06 cM, which facilitated covering the vrs1 region by a 518 kb barley BAC contig. An analysis of the contig sequence revealed that a rice Vrs1 orthologue is present on chromosome 7, suggesting a transposition of the chromosomal segment containing Vrs1 within barley chromosome 2H. The breakdown of micro-collinearity illustrates the limitations of synteny cloning, and stresses the importance of implementing genomic studies directly in the target species.

Analysis of the barley chromosome 2 region containing the six-rowed spike gene vrs1 reveals a breakdown of rice-barley micro collinearity by a transposition

In cultivated barley (Hordeum vulgare ssp. vulgare), six-rowed spikes produce three times as many seeds per spike as do two-rowed spikes. The determinant of this trait is the Mendelian gene vrs1, located on chromosome 2H, which is syntenous with rice (Oryza sativa) chromosomes 4 and 7. We exploited barley-rice micro-synteny to increase marker density in the vrs1 region as a prelude to its map-based cloning. The rice genomic sequence, covering a 980 kb contig, identified barley ESTs linked to vrs1. A high level of conservation of gene sequence was obtained between barley chromosome 2H and rice chromosome 4. A total of 22 EST-based STS markers were placed within the target region, and the linear order of these markers in barley and rice was identical. The genetic window containing vrs1 was narrowed from 0.5 to 0.06 cM, which facilitated covering the vrs1 region by a 518 kb barley BAC contig. An analysis of the contig sequence revealed that a rice Vrs1 orthologue is present on chromosome 7, suggesting a transposition of the chromosomal segment containing Vrs1 within barley chromosome 2H. The breakdown of micro-collinearity illustrates the limitations of synteny cloning, and stresses the importance of implementing genomic studies directly in the target species.

Triticeae genomics: advances in sequence analysis of large genome cereal crops

Whole genome sequencing provides direct access to all genes of an organism and represents an essential step towards a systematic understanding of (crop) plant biology. Wheat and barley, two of the most important crop species worldwide, have two- to five-fold larger genomes than human - too large to be completely sequenced at current costs. Nevertheless, significant progress has been made to unlock the gene contents of these species by sequencing expressed sequence tags (EST) for high-density mapping and as a basis for elucidating gene function on a large scale. Several megabases of genomic (BAC) sequences have been obtained providing a first insight into the complexity of these huge cereal genomes. However, to fully exploit the information of the wheat and barley genomes for crop improvement, sequence analysis of a significantly larger portion of the Triticeae genomes is needed. In this review an overview of the current status of Triticeae genome sequencing and a perspective concerning future developments in cereal structural genomics is provided.

The promoter of the wheat puroindoline-a gene (PinA) exhibits a more complex pattern of activity than that of the PinB gene and is induced by wounding and pathogen attack in rice

Puroindolines form the molecular basis of wheat grain hardness. However, little is known about puroindoline gene regulation. We previously reported that the Triticum aestivum puroindoline-b gene (PinB) promoter directs beta-glucuronidase gene (uidA) seed-specific expression in transgenic rice. In this study, we isolated a puroindoline-a gene (PinA), analyzed PinA promoter activity by 5' deletions and compared PinA and PinB promoters in transgenic rice. Seeds of PinA-1214 and PinB-1063 transgenic plants strongly expressed uidA in endosperm, in the aleurone layer and in epidermis cells in a developmentally regulated manner. The GUS activity was also observed in PinA-1214 embryos. Whereas the PinB promoter is seed specific, the PinA promoter also directed, but to a lower level, uidA expression in roots of seedlings and in the vascular tissues of palea and pollen grains of dehiscent anthers during flower development. In addition, the PinA promoter was induced by wounding and by Magnaporthe grisea. By deletion analysis, we showed that the "390-bp" PinA promoter drives the same expression pattern as the "1214-bp" promoter. Moreover, the "214-bp" PinA promoter drives uidA expression solely in pollen grains of dehiscent anthers. The presence of putative cis-regulatory elements that may be related to PinA expression is discussed from an evolutionary point of view. By electrophoretic mobility shift assay, we showed that putative cis-elements (WUN-box, TCA motifs and as-1-like binding sites) whose presence in the PinA promoter may be related to wounding and/or the pathogen response form complexes with nuclear extracts isolated from wounded wheat leaves.

Recent history of artificial outcrossing facilitates whole-genome association mapping in elite inbred crop varieties

Genomewide association studies depend on the extent of linkage disequilibrium (LD), the number and distribution of markers, and the underlying structure in populations under study. Outbreeding species generally exhibit limited LD, and consequently, a very large number of markers are required for effective whole-genome association genetic scans. In contrast, several of the world's major food crops are self-fertilizing inbreeding species with narrow genetic bases and theoretically extensive LD. Together these are predicted to result in a combination of low resolution and a high frequency of spurious associations in LD-based studies. However, inbred elite plant varieties represent a unique human-induced pseudo-outbreeding population that has been subjected to strong selection for advantageous alleles. By assaying 1,524 genomewide SNPs we demonstrate that, after accounting for population substructure, the level of LD exhibited in elite northwest European barley, a typical inbred cereal crop, can be effectively exploited to map traits by using whole-genome association scans with several hundred to thousands of biallelic SNPs.

High-frequency Ds remobilization over multiple generations in barley facilitates gene tagging in large genome cereals

Transposable elements have certain advantages over other approaches for identifying and determining gene function in large genome cereals. Different strategies have been used to exploit the maize Activator/dissociation (Ac/Ds) transposon system for functional genomics in heterologous species. Either large numbers of independent Ds insertion lines or transposants (TNPs) are generated and screened phenotypically, or smaller numbers of TNPs are produced, Ds locations mapped and remobilized for localized gene targeting. It is imperative to characterize key features of the system in order to utilize the latter strategy, which is more feasible in large genome cereals like barley and wheat. In barley, we generated greater than 100 single-copy Ds TNPs and determined remobilization frequencies of primary, secondary, and tertiary TNPs with intact terminal inverted repeats (TIRs); frequencies ranged from 11.8 to 17.1%. In 16% of TNPs that had damaged TIRs no transposition was detected among progeny of crosses using those TNPs as parental lines. In half of the greater than 100 TNP lines, the nature of flanking sequences and status of the 11 bp TIRs and 8-bp direct repeats were determined. BLAST searches using a gene prediction program revealed that 86% of TNP flanking sequences matched either known or putative genes, indicating preferential Ds insertion into genic regions, critical in large genome species. Observed remobilization frequencies of primary, secondary, tertiary, and quaternary TNPs, coupled with the tendency for localized Ds transposition, validates a saturation mutagenesis approach using Ds to tag and characterize genes linked to Ds in large genome cereals like barley and wheat.

Molecular evolution of the puroindoline-a, puroindoline-b, and grain softness protein-1 genes in the tribe Triticeae

The genome organization of the Hardness locus in the tribe Triticeae constitutes an excellent model for studying the mechanisms of evolution that played a role in the preservation and potential functional innovations of duplicate genes. Here we applied the nonsynonymous-synonymous rate ratio (d ( N )/d ( S ) or omega) to measure the selective pressures at the paralogous puroindoline-a (Pina), puroindoline-b (Pinb), and grain softness protein-1 (Gsp-1) genes located at this locus. Puroindolines represent the molecular-genetic basis of grain texture. In addition, the puroindoline gene products have antimicrobial properties with potential role in plant defense. We document the complete coding sequences from the Triticum/Aegilops taxa, rye and barley including the A, D, C, H, M, N, R, S, and U genomes of the Triticeae. Maximum likelihood analyses performed on Bayesian phylogenetic trees showed distinct evolutionary patterns among Pina, Pinb, and Gsp-1. Positive diversifying selection appeared to drive the evolution of at least one of the three genes examined, suggesting that adaptive forces have operated at this locus. Results evidenced positive selection (omega > 4) at Pina and detected amino acid residues along the mature PIN-a protein with a high probability (>95%) of having evolved under adaptation. We hypothesized that positive selection at the Pina region is congruent with its role as a plant defense gene.

Genome-wide SNP discovery and linkage analysis in barley based on genes responsive to abiotic stress

More than 2,000 genome-wide barley single nucleotide polymorphisms (SNPs) were developed by resequencing unigene fragments from eight diverse accessions. The average genome-wide SNP frequency observed in 877 unigenes was 1 SNP per 200 bp. However, SNP frequency was highly variable with the least number of SNP and SNP haplotypes observed within European cultivated germplasm reflecting effects of breeding history on genetic diversity. More than 300 SNP loci were mapped genetically in three experimental mapping populations which allowed the construction of an integrated SNP map incorporating a large number of RFLP, AFLP and SSR markers (1,237 loci in total). The genes used for SNP discovery were selected based on their transcriptional response to a variety of abiotic stresses. A set of known barley abiotic stress QTL was positioned on the linkage map, while the available sequence and gene expression information facilitated the identification of genes potentially associated with these traits. Comparison of the sequenced SNP loci to the rice genome sequence identified several regions of highly conserved gene order providing a framework for marker saturation in barley genomic regions of interest. The integration of genome-wide SNP and expression data with available genetic and phenotypic information will facilitate the identification of gene function in barley and other non-model organisms.

Genome-wide SNP discovery and linkage analysis in barley based on genes responsive to abiotic stress

More than 2,000 genome-wide barley single nucleotide polymorphisms (SNPs) were developed by resequencing unigene fragments from eight diverse accessions. The average genome-wide SNP frequency observed in 877 unigenes was 1 SNP per 200 bp. However, SNP frequency was highly variable with the least number of SNP and SNP haplotypes observed within European cultivated germplasm reflecting effects of breeding history on genetic diversity. More than 300 SNP loci were mapped genetically in three experimental mapping populations which allowed the construction of an integrated SNP map incorporating a large number of RFLP, AFLP and SSR markers (1,237 loci in total). The genes used for SNP discovery were selected based on their transcriptional response to a variety of abiotic stresses. A set of known barley abiotic stress QTL was positioned on the linkage map, while the available sequence and gene expression information facilitated the identification of genes potentially associated with these traits. Comparison of the sequenced SNP loci to the rice genome sequence identified several regions of highly conserved gene order providing a framework for marker saturation in barley genomic regions of interest. The integration of genome-wide SNP and expression data with available genetic and phenotypic information will facilitate the identification of gene function in barley and other non-model organisms.

Extreme population-dependent linkage disequilibrium detected in an inbreeding plant species, Hordeum vulgare

In human genetics a detailed knowledge of linkage disequilibrium (LD) is considered a prerequisite for effective population-based, high-resolution gene mapping and cloning. Similar opportunities exist for plants; however, differences in breeding system and population history need to be considered. Here we report a detailed study of localized LD in different populations of an inbreeding crop species. We measured LD between and within four gene loci within the region surrounding the hardness locus in three different gene pools of barley (Hordeum vulgare). We demonstrate that LD extends to at least 212 kb in elite barley cultivars but is rapidly eroded in related inbreeding ancestral populations. Our results indicate that haplotype-based sequence analysis in multiple populations will provide new opportunities to adjust the resolution of association studies in inbreeding crop species.

Novel puroindoline and grain softness protein alleles in Aegilops species with the C, D, S, M and U genomes

DNA from six hexaploid, tetraploid and diploid species of Aegilops with the C, D, S, M and U genomes was amplified with specific PCR primers to identify sequences encoding puroindolines (Pins) a and b and grain softness protein (GSP), all of which are encoded by genes at the Ha (hardness) locus, with Ae. tauschii (DD) and bread wheat (T. aestivum) (AABBDD) cv Hiline being studied as controls. Seven new allelic forms of Pin a and Pin b were identified, including forms with mutations within or close to the tryptophan motif. In addition, five new forms of GSP were detected. In all species both genomic DNA from leaves and cDNA from developing grain were analysed. This revealed the presence of both silent genes (with premature stop codons) and multiple genes, with the latter being confirmed by Southern blot analysis. Freeze fracture analysis demonstrated that all except one accession (Ae. sharonensis) were soft textured. However, this difference cannot be accounted for by the sequences of the Pin alleles present in this line.

Single-feature polymorphism discovery in the barley transcriptome

A probe-level model for analysis of GeneChip gene-expression data is presented which identified more than 10,000 single-feature polymorphisms (SFP) between two barley genotypes. The method has good sensitivity, as 67% of known single-nucleotide polymorphisms (SNP) were called as SFPs. This method is applicable to all oligonucleotide microarray data, accounts for SNP effects in gene-expression data and represents an efficient and versatile approach for highly parallel marker identification in large genomes.

Single-feature polymorphism discovery in the barley transcriptome

A probe level model for analysis of GeneChip gene expression data is presented which identified more than 10,000 single-feature polymorphisms between two barley genotypes, with a high sensitivity. This method is applicable to all oligonucleotide microarray data.

A probe-level model for analysis of GeneChip gene-expression data is presented which identified more than 10,000 single-feature polymorphisms (SFP) between two barley genotypes. The method has good sensitivity, as 67% of known single-nucleotide polymorphisms (SNP) were called as SFPs. This method is applicable to all oligonucleotide microarray data, accounts for SNP effects in gene-expression data and represents an efficient and versatile approach for highly parallel marker identification in large genomes.

Molecular basis of evolutionary events that shaped the hardness locus in diploid and polyploid wheat species (Triticum and Aegilops)

The Hardness (Ha) locus controls grain hardness in hexaploid wheat (Triticum aestivum) and its relatives (Triticum and Aegilops species) and represents a classical example of a trait whose variation arose from gene loss after polyploidization. In this study, we investigated the molecular basis of the evolutionary events observed at this locus by comparing corresponding sequences of diploid, tertraploid, and hexaploid wheat species (Triticum and Aegilops). Genomic rearrangements, such as transposable element insertions, genomic deletions, duplications, and inversions, were shown to constitute the major differences when the same genomes (i.e., the A, B, or D genomes) were compared between species of different ploidy levels. The comparative analysis allowed us to determine the extent and sequences of the rearranged regions as well as rearrangement breakpoints and sequence motifs at their boundaries, which suggest rearrangement by illegitimate recombination. Among these genomic rearrangements, the previously reported Pina and Pinb genes loss from the Ha locus of polyploid wheat species was caused by a large genomic deletion that probably occurred independently in the A and B genomes. Moreover, the Ha locus in the D genome of hexaploid wheat (T. aestivum) is 29 kb smaller than in the D genome of its diploid progenitor Ae. tauschii, principally because of transposable element insertions and two large deletions caused by illegitimate recombination. Our data suggest that illegitimate DNA recombination, leading to various genomic rearrangements, constitutes one of the major evolutionary mechanisms in wheat species.

Specific patterns of changes in wheat gene expression after treatment with three antifungal compounds

The two fungicides azoxystrobin and fenpropimorph are used against powdery mildew and rust diseases in wheat (Triticum aestivumL). Azoxystrobin, a strobilurin, inhibits fungal mitochondrial respiration and fenpropimorph, a morpholin, represses biosynthesis of ergosterol, the major sterol of fungal membranes. Although the fungitoxic activity of these compounds is well understood, their effects on plant metabolism remain unclear. In contrast to the fungicides which directly affect pathogen metabolism, benzo(1,2,3) thiadiazole-7-carbothioic acid S-methylester (BTH) induces resistance against wheat pathogens by the activation of systemic acquired resistance in the host plant. In this study, we monitored gene expression in spring wheat after treatment with each of these agrochemicals in a greenhouse trial using a microarray containing 600 barley cDNA clones. Defence-related genes were strongly induced after treatment with BTH, confirming the activation of a similar set of genes as in dicot plants following salicylic acid treatment. A similar gene expression pattern was observed after treatment with fenpropimorph and some defence-related genes were induced by azoxystrobin, demonstrating that these fungicides also activate a defence reaction. However, less intense responses were triggered than with BTH. The same experiments performed under field conditions gave dramatically different results. No gene showed differential expression after treatment and defence genes were already expressed at a high level before application of the agrochemicals. These differences in the expression patterns between the two environments demonstrate the importance of plant growth conditions for testing the impact of agrochemicals on plant metabolism.