Identification and analysis of expressed resistance gene sequences in wheat

Forty-eight resistance (R) genes conferring resistance to various types of pests have been cloned from 12 plant species. Irrespective of the host or the pest type, most R genes share a strong protein sequence similarity especially for domains and motifs. The objective of this study was to identify expressed R genes of wheat, the fraction of which is expected to be very low in the genome. Using modified RNA fingerprinting and data mining approaches we identified 220 expressed R-gene candidates. Of these, 125 sequences structurally resembled known R genes. In addition to 25-87% protein sequence similarity with the known R genes, the sequence, order, and distribution of the domains and motifs were also the same. Among the remaining 95, 17 were probable R-related, 21 were a new class of nucleotide-binding kinases, 21 were probable kinases, and 36 were p-loop-containing unknown sequences. About 76% were rare including 73 novel sequences. Three new R-gene specific motifs were also identified. Physical mapping of the 164 best R-gene candidates on 339 deletion lines localized 121 mappable R-gene candidates to 26 small chromosomal regions encompassing about 16% of the genome. About 90 of the 110 phenotypically characterized wheat R genes corresponding to 18 different pests also mapped in these regions.     

Arabidopsis thaliana: a source of candidate disease-resistance genes for Brassica napus.

Common structural and amino acid motifs among cloned plant disease-resistance genes (R genes), have made it possible to identify putative disease-resistance sequences based on DNA sequence identity. Mapping of such R-gene homologues will identify candidate disease-resistance loci to expedite map-based cloning strategies in complex crop genomes. Arabidopsis thaliana expressed sequence tags (ESTs) with homology to cloned plant R genes (R-ESTs), were mapped in both A. thaliana and Brassica napus to identify candidate R-gene loci and investigate intergenomic collinearity. Brassica R-gene homologous sequences were also mapped in B. napus. In total, 103 R-EST loci and 36 Brassica R-gene homologous loci were positioned on the N-fo-61-9 B. napus genetic map, and 48 R-EST loci positioned on the Columbia x Landsberg A. thaliana map. The mapped loci identified collinear regions between Arabidopsis and Brassica which had been observed in previous comparative mapping studies; the detection of syntenic genomic regions indicated that there was no apparent rapid divergence of the identified genomic regions housing the R-EST loci.     

Development of a new wheat microarray from a durum wheat totipotent cDNA library used for a powdery mildew resistance study

Totipotent cDNA libraries representative of all the potentially expressed sequences in a genome would be of great benefit to gene expression studies. Here, we report on an innovative method for creating such a library for durum wheat (Triticum turgidum L. var. durum) and its application for gene discovery. The use of suitable quantities of 5-azacytidine during the germination phase induced the demethylation of total DNA, and the resulting seedlings potentially express all of the genes present in the genome. A new wheat microarray consisting of 4925 unigenes was developed from the totipotent cDNA library and used to screen for genes that may contribute to differences in the disease resistance of two near-isogenic lines, the durum wheat cultivar Latino and the line 5BIL-42, which are respectively susceptible and resistant to powdery mildew. Fluorescently labeled cDNA was prepared from the RNA of seedlings of the two near-isogenic wheat lines after infection with a single powdery mildew isolate under controlled conditions in the greenhouse. Hybridization to the microarray identified six genes that were differently expressed in the two lines. Four of the sequences could be assigned putative functions based on their similarity to known genes in public databases. Physical mapping of the six genes localized them to two regions of the genome: the centromeric region of chromosome 5B, where the Pm36 resistance gene was previously localized, and chromosome 6B.     

A resistance-like gene identified by EST mapping and its association with a QTL controlling Fusarium head blight infection on wheat chromosome 3BS.

Fusarium head blight (FHB) is a major disease in the wheat growing regions of the world. A quantitative trait locus (QTL) on the short arm of chromosome 3B controls much of the variation for resistance. The cloning of candidate disease-resistance genes for FHB QTLs on chromosome 3B can provide further elucidation of the mechanisms that control resistance. However, rearrangements and divergence during plant genome evolution often hampers the identification of sequences with similarity to known disease-resistance genes. This study focuses on the use of wheat expressed sequence tags (ESTs) that map to the region on chromosome 3B containing the QTL for FHB resistance and low-stringency BLAST searching to identify sequences with similarity to known disease-resistance genes. One EST rich with leucine repeats and low similarity to a protein kinase domain of the barley Rpg1 gene was identified. Genetic mapping using a Ning894037 x Alondra recombinant inbred (RI) population showed that this EST mapped to the QTL on the short arm of chromosome 3B and may represent a portion of a newly diverged gene contributing to FHB resistance. The EST is a new marker suitable for marker-assisted selection and provides a starting point to begin map-based cloning for chromosome walking and investigate new diverged genes at this locus.     

Characterization and linkage mapping of R-gene analogous DNA sequences in pea (Pisum sativum L.)

Pea (Pisum sativum L.) sequences that are analogous to the conserved nucleotide binding site (NBS) domain found in a number of plant disease resistance genes (R-genes) were cloned. Using redundant oligonucleotide primers and the polymerase chain reaction (PCR), we amplified nine pea sequences and characterised their sequences. The pea R-gene analog (RGA)- deduced amino acid sequences demonstrated significant sequence similarity with known R-gene sequences lodged in public databases. The genomic locations of eight of the pea RGAs were determined by linkage mapping. The eight RGAs identified ten loci that mapped to six linkage groups. In addition, the genomic organization of the RGAs was inferred. Both single-copy and multicopy sequence families were present among the RGAs, and the multicopy families occurred most often as tightly linked clusters of related sequences. Intraspecific copy number variability was observed in three of the RGA sequence families, suggesting that these sequence families are evolving rapidly. The genomic locations of the pea RGAs were compared with the locations of known pea R-genes and sym genes involved in the pea-rhizobia symbiosis. Two pea RGAs mapped in the genomic region containing a pea R-gene, Fw, and four pea RGAs mapped in regions of the genome containing sym genes. Received: 4 August 1999 / Accepted: 11 November 1999     

Structure, function and evolution of plant disease resistance genes

Gene-for-gene plant disease resistance involves two basic processes: perception of pathogen attack, followed by responses to limit disease. Perception involves receptors with high degrees of specificity for pathogen strains, which are encoded by disease resistance genes. Large repertoires of distantly related resistance (R) genes with diverse recognitional specificities are found within a single plant species. The generation of R-gene polymorphism involves gene duplication, followed by DNA-sequence divergence by point mutation, and by deletion and duplication of intragenic DNA repeats encoding blocks of leucine-rich elements. Recombination between related genes reassorts this variation to further diversify gene sequences. Pathogen pressure selects functional resistance specificities and results in the maintenance of R-gene diversity. Recent genome-sequence data reveal that the NBS-LRR (i.e. nucleotide-binding site-leucine-rich repeat) class of R genes represents as much as 1% of the Arabidopsis genome. Experimental data have shown that the LRR has a role in determination of specificity. Mutation experiments, in which R-gene signaling has been dissociated from specificity in constitutive signal mutants, have provided the potential for non-specific resistance to be expressed from pathogen-infection-induced promoters in transgenic plants.     

Macro- and microcolinearity between the genomic region of wheat chromosome 5B containing the Tsn1 gene and the rice genome

The Tsn1 gene in wheat confers sensitivity to a proteinaceous host-selective toxin (Ptr ToxA) produced by the tan spot fungus (Pyrenophora tritici-repentis) and lies within a gene-rich region of chromosome 5B. To use the rice genome sequence information for the map-based cloning of Tsn1, colinearity between the wheat genomic region containing Tsn1 and the rice genome was determined at the macro- and microlevels. Macrocolinearity was determined by testing 28 expressed sequence markers (ESMs) spanning a 25.5-cM segment and encompassing Tsn1 for similarity to rice sequences. Twelve ESMs had no similarity to rice sequences, and 16 had similarity to sequences on seven different rice chromosomes. Segments of colinearity with rice chromosomes 3 and 9 were identified, but frequent rearrangements and disruptions occurred. Microcolinearity was determined by testing the sequences of 26 putative genes identified from BAC contigs of 205 and 548 kb in length and flanking Tsn1 for similarity to rice genomic sequences. Fourteen of the predicted genes detected orthologous sequences on six different rice chromosomes, whereas the remaining 12 had no similarity with rice sequences. Four genes were colinear on rice chromosome 9, but multiple disruptions, rearrangements, and duplications were observed in wheat relative to rice. The data reported provide a detailed analysis of a region of wheat chromosome 5B that is highly rearranged relative to rice.     

Comparative mapping of the barley Ppd-H1 photoperiod response gene region,which lies close to a junction between two rice linkage segments

Comparative mapping of cereals has shown that chromosomes of barley, wheat, and maize can be described in terms of rice "linkage segments." However, little is known about marker order in the junctions between linkage blocks or whether this will impair comparative analysis of major genes that lie in such regions. We used genetic and physical mapping to investigate the relationship between the distal part of rice chromosome 7L, which contains the Hd2 heading date gene, and the region of barley chromosome 2HS containing the Ppd-H1 photoperiod response gene, which lies near the junction between rice 7 and rice 4 linkage segments. RFLP markers were mapped in maize to identify regions that might contain Hd2 or Ppd-H1 orthologs. Rice provided useful markers for the Ppd-H1 region but comparative mapping was complicated by loss of colinearity and sequence duplications that predated the divergence of rice, maize, and barley. The sequences of cDNA markers were used to search for homologs in the Arabidopsis genome. Homologous sequences were found for 13 out of 16 markers but they were dispersed in Arabidopsis and did not identify any candidate equivalent region. The implications of the results for comparative trait mapping in junction regions are discussed.     

Functional features of a single chromosome arm in wheat (1AL) determined from its structure

Bread wheat (Triticum aestivum L.) is one of the most important crops globally and a high priority for genetic improvement, but its large and complex genome has been seen as intractable to whole genome sequencing. Isolation of individual wheat chromosome arms has facilitated large-scale sequence analyses. However, so far there is no such survey of sequences from the A genome of wheat. Greater understanding of an A chromosome could facilitate wheat improvement and future sequencing of the entire genome. We have constructed BAC library from the long arm of T. aestivum chromosome 1A (1AL) and obtained BAC end sequences from 7,470 clones encompassing the arm. We obtained 13,445 (89.99%) useful sequences with a cumulative length of 7.57 Mb, representing 1.43% of 1AL and about 0.14% of the entire A genome. The GC content of the sequences was 44.7%, and 90% of the chromosome was estimated to comprise repeat sequences, while just over 1% encoded expressed genes. From the sequence data, we identified a large number of sites suitable for development of molecular markers (362 SSR and 6,948 ISBP) which will have utility for mapping this chromosome and for marker assisted breeding. From 44 putative ISBP markers tested 23 (52.3%) were found to be useful. The BAC end sequence data also enabled the identification of genes and syntenic blocks specific to chromosome 1AL, suggesting regions of particular functional interest and targets for future research.     

Mapping of EST- and STS-markers in the human genome using a panel fo radiation hybrids

Ten DNA markers were localized in the human genome by a screening procedure against the radiation hybrid somatic cell panel (GeneBridge 4 RH Panel) using polymerase chain reaction (RH mapping method). DNA markers were developed to nucleotide sequences adjacent to NotI sites of human chromosome 3 (NotI-STS markers) and also to nucleotide sequences of human cDNA (EST markers). Three EST markers mapped (B10164, S16R and 18F5R) were localized in the human genome for the first time. Marker B10164 was found to be homologous to the nucleotide sequence of the BASP1 gene coding a major receptor protein. Markers S16R and 18F5R presumably tagged new genes, because no homologies were revealed among the nucleotide sequences presented in the databases. For four NotI-STS, more precise localization on human chromosome 3 was determined. On the basis of the data obtained, the NotI map may be integrated with other types of physical maps of human chromosome 3. RH mapping with a standard commercial panel of radiation hybrid somatic cells provided a chance to integrate the data obtained into international databases and existing integrated human chromosomal maps.     

Characterization of wheat DArT markers: genetic and functional features

Diversity array technology (DArT) markers are largely used for mapping, genetic diversity, and association mapping studies. For years, they have been used as anonymous genomic markers, as their sequences were not known. As the sequences of 2,000 wheat DArT clones are now available, this study was designed to analyze these sequences with bioinformatic approaches, and to study the genetic features of a subset of 291 markers positioned on the A and B genomes in three durum wheat genetic maps. A set of 1,757 non-redundant sequences was identified, and used as queries for similarity searches. Analysis of the genetic positions of markers corresponding to nearly identical sequences indicates that redundancy of sequences is one of the factors that explains the clustering of these markers in specific genomic regions. Of a total of 1,124 DArT clones (64?%) that represent putatively expressed sequences, putative functions are proposed for more than 700 of them. Of note, many clones correspond to genes that are related to disease resistance, as characterized by leucine-rich repeat domains, and 40 of these clones are positioned in the three genetic maps presented in this study. Finally, DArT markers have been used to find syntenic regions in the Brachypodium and rice genomes. In conclusion, the analyses herein presented contribute to explain the main features of DArT markers observed in genetic maps, as clustering in short chromosome regions. Moreover, the attribution of putative gene functions for more than 700 sequences makes these markers an optimal tool for collinearity studies or for the identification of candidate genes.     

Comparative genetics of nucleotide binding site-leucine rich repeat resistance gene homologues in the genomes of two dicotyledons: tomato and arabidopsis

The presence of a single resistance (R) gene allele can determine plant disease resistance. The protein products of such genes may act as receptors that specifically interact with pathogen-derived factors. Most functionally defined R-genes are of the nucleotide binding site-leucine rich repeat (NBS-LRR) supergene family and are present as large multigene families. The specificity of R-gene interactions together with the robustness of plant-pathogen interactions raises the question of their gene number and diversity in the genome. Genomic sequences from tomato showing significant homology to genes conferring race-specific resistance to pathogens were identified by systematically "scanning" the genome using a variety of primer pairs based on ubiquitous NBS motifs. Over 70 sequences were isolated and 10% are putative pseudogenes. Mapping of the amplified sequences on the tomato genetic map revealed their organization as mixed clusters of R-gene homologues that showed in many cases linkage to genetically characterized tomato resistance loci. Interspecific examination within Lycopersicon showed the existence of a null allele. Consideration of the tomato and potato comparative genetic maps unveiled conserved syntenic positions of R-gene homologues. Phylogenetic clustering of R-gene homologues within tomato and other Solanaceae family members was observed but not with R-gene homologues from Arabidopsis thaliana. Our data indicate remarkably rapid evolution of R-gene homologues during diversification of plant families.     

Detailed comparative gene map of rat chromosome 1 with mouse and human genomes and physical mapping of an evolutionary chromosomal breakpoint

We report the localization of 92 new gene-based markers assigned to rat chromosome 1 by linkage or radiation hybrid mapping. The markers were chosen to enrich gene mapping data in a region of the rat chromosome known to contain several of the principal quantitative trait loci in rodent models of human multifactorial disease. The composite map reported here provides map information on a total of 139 known genes, including 80 that have been localized in mouse and 109 that have been localized in human, and integrates the gene-based markers with anonymous microsatellites. The evolutionary breakpoints identifying 16 segments that are homologous regions in the human genome are defined. These data will facilitate genetic and comparative mapping studies and identification of novel candidate genes for the quantitative trait loci that have been localized to the region.     

Efficient targeting of plant disease resistance loci using NBS profiling

The conserved sequences in the nucleotide-binding sites of the nucleotide-binding site-leucine-rich repeat (NBS-LRR) class of disease resistance (R) genes have been used for PCR-based R-gene isolation and subsequent development of molecular markers. Here we present a PCR-based approach (NBS profiling) that efficiently targets R genes and R-gene analogs (RGAs) and, at the same time, produces polymorphic markers in these genes. In NBS profiling, genomic DNA is digested with a restriction enzyme, and an NBS-specific (degenerate) primer is used in a PCR reaction towards an adapter linked to the resulting DNA fragments. The NBS profiling protocol generates a reproducible polymorphic multilocus marker profile on a sequencing gel that is highly enriched for R genes and RGAs. NBS profiling was successfully used in potato with several restriction enzymes, and several primers targeted to different conserved motifs in the NBS. Across primers and enzymes, the NBS profiles contained 50–90% fragments that were significantly similar to known R-gene and RGA sequences. The protocol was similarly successful in other crops (including tomato, barley, and lettuce) without modifications. NBS profiling can thus be used to produce markers tightly linked to R genes and R-gene clusters for genomic mapping and positional cloning and to mine for new alleles and new sources of disease resistance in available germplasm.Communicated by H.F. Linskens     

High-density mapping and comparative analysis of agronomically important traits on wheat chromosome 3A

Bread wheat chromosome 3A has been shown to contain genes/QTLs controlling grain yield and other agronomic traits. The objectives of this study were to generate high-density physical and genetic-linkage maps of wheat homoeologous group 3 chromosomes and reveal the physical locations of genes/QTLs controlling yield and its component traits, as well as agronomic traits, to obtain a precise estimate of recombination for the corresponding regions and to enrich the QTL-containing regions with markers. Physical mapping was accomplished by 179 DNA markers mostly representing expressed genes using 41 single-break deletion lines. Polymorphism survey of cultivars Cheyenne (CNN) and Wichita (WI), and a substitution line of CNN carrying chromosome 3A from WI [CNN(WI3A)], with 142 RFLP probes and 55 SSR markers revealed that the extent of polymorphism is different among various group 3 chromosomal regions as well as among the homoeologs. A genetic-linkage map for chromosome 3A was developed by mapping 17 QTLs for seven agronomic traits relative to 26 RFLP and 15 SSR chromosome 3A-specific markers on 95 single-chromosome recombinant inbred lines. Comparison of the physical maps with the 3A genetic-linkage map localized the QTLs to gene-containing regions and accounted for only about 36% of the chromosome. Two chromosomal regions containing 9 of the 17 QTLs encompassed less than 10% of chromosome 3A but accounted for almost all of the arm recombination. To identify rice chromosomal regions corresponding to the particular QTL-containing wheat regions, 650 physically mapped wheat group 3 sequences were compared with rice genomic sequences. At an E value of E < or = 10(-5), 82% of the wheat group 3 sequences identified rice homologs, of which 54% were on rice chromosome 1. The rice chromosome 1 region collinear with the two wheat regions that contained 9 QTLs was about 6.5 Mb.     

Genome-wide linkage disequilibrium analysis in bread wheat and durum wheat.

Bread wheat and durum wheat were examined for linkage disequilibrium (LD) using microsatellite markers distributed across the genome. The allele database consisted of 189 bread wheat accessions genotyped at 370 loci and 93 durum wheat accessions genotyped at 245 loci. A significance level of p < 0.001 was set for all comparisons. The bread and durum wheat collections showed that 47.9% and 14.0% of all locus pairs were in LD, respectively. LD was more prevalent between loci on the same chromosome compared with loci on independent chromosomes and was highest between adjacent loci. Only a small fraction (bread wheat, 0.9%; durum wheat, 3.2%) of the locus pairs in LD showed R2 values > 0.2. The LD between adjacent locus pairs extended (R2 > 0.2) approximately 2-3 cM, on average, but some regions of the bread and durum wheat genomes showed high levels of LD (R2 = 0.7 and 1.0, respectively) extending 41.2 and 25.5 cM, respectively. The wheat collections were clustered by similarity into subpopulations using unlinked microsatellite data and the software Structure. Analysis within subpopulations showed 14- to 16-fold fewer locus pairs in LD, higher R2 values for those pairs in LD, and LD extending further along the chromosome. The data suggest that LD mapping of wheat can be performed with simple sequence repeats to a resolution of <5 cM.     

Characterization and Fine Localization of Two New Genes in Xq28 Using the Genomic Sequence/EST Database Screening Approach

Two new genes were identified and mapped by searching the EST databases with genomic sequences obtained from putative CpG islands of the rodent–human hybrid X3000. Previous mapping of these CpG islands in the proximity of the host cell factor (HCFC1) and GdX genes automatically localized these two new genes to Xq28 in the interval between the L1 cell adhesion molecule (L1CAM) and the glucose-6-phosphate dehydrogenase (G6PD) loci. Both genes are relatively short, contain an ORF of 261 and 105 amino acids, respectively, and are ubiquitously expressed. Combining sequencing of selected CpG islands, derived from hybrids containing small portions of the human genome, with an EST database search is an easy method of identifying and mapping new genes to specific regions of the genome.     

An EST survey of the sugarcane transcriptome

Its large genome and high polyploidy makes sugarcane (Saccharum spp.) a singularly challenging crop to study and improve using genetic approaches. To provide large numbers of functionally characterized candidate genes that might be tested for direct association (rather than distant linkage) with economically important traits, we sequenced the 5' ends of 9,216 clones from three cDNA libraries (apex, leaf and mature internode), representing 3,401 non-redundant sequences. About 57% of these sequences could be assigned a tentative function based on statistically significant similarity to previously characterized proteins or DNA sequences. Another 28% corresponded to previously identified, but uncharacterized, sequences. Some of the remaining unidentified sequences were predicted to be genes which could potentially be new to plants or unique to sugarcane. Comparisons of the sugarcane ESTs to a large sorghum EST database revealed similar compositions of expressed genes between some different tissues. Comparison to a detailed Arabidopsis protein database showed some highly conserved sequences, which might be useful DNA markers for pan-angiosperm comparative mapping. These EST sequences provide a foundation for many new studies to accelerate isolation of agronomically important genes from the cumbersome sugarcane genome.