Genome-wide SNPs and re-sequencing of growth habit and inflorescence genes in barley: implications for association mapping in germplasm arrays varying in size and structure

BACKGROUND

Considerations in applying association mapping (AM) to plant breeding are population structure and size: not accounting for structure and/or using small populations can lead to elevated false-positive rates. The principal determinants of population structure in cultivated barley are growth habit and inflorescence type. Both are under complex genetic control: growth habit is controlled by the epistatic interactions of several genes. For inflorescence type, multiple loss-of-function alleles in one gene lead to the same phenotype. We used these two traits as models for assessing the effectiveness of AM. This research was initiated using the CAP Core germplasm array (n = 102) assembled at the start of the Barley Coordinated Agricultural Project (CAP). This array was genotyped with 4,608 SNPs and we re-sequenced genes involved in morphology, growth and development. Larger arrays of breeding germplasm were subsequently genotyped and phenotyped under the auspices of the CAP project. This provided sets of 247 accessions phenotyped for growth habit and 2,473 accessions phenotyped for inflorescence type. Each of the larger populations was genotyped with 3,072 SNPs derived from the original set of 4,608.

RESULTS

Significant associations with SNPs located in the vicinity of the loci involved in growth habit and inflorescence type were found in the CAP Core. Differentiation of true and spurious associations was not possible without a priori knowledge of the candidate genes, based on re-sequencing. The re-sequencing data were used to define allele types of the determinant genes based on functional polymorphisms. In a second round of association mapping, these synthetic markers based on allele types gave the most significant associations. When the synthetic markers were used as anchor points for analysis of interactions, we detected other known-function genes and candidate loci involved in the control of growth habit and inflorescence type. We then conducted association analyses--with SNP data only--in the larger germplasm arrays. For both vernalization sensitivity and inflorescence type, the most significant associations in the larger data sets were found with SNPs coincident with the synthetic markers used in the CAP Core and with SNPs detected via interaction analysis in the CAP Core.

CONCLUSIONS

Small and highly structured collections of germplasm, such as the CAP Core, are cost-effectively phenotyped and genotyped with high-throughput markers. They are also useful for characterizing allelic diversity at loci in germplasm of interest. Our results suggest that discovery-oriented exercises in AM in such small arrays may generate a large number of false-positives. However, if haplotypes in candidate genes are available, they may be used as anchors in an analysis of interactions to identify other candidate regions harboring genes determining target traits. Using larger germplasm arrays, genome regions where the principal genes determining vernalization sensitivity and row type are located were identified.

Genome-wide SNPs and re-sequencing of growth habit and inflorescence genes in barley: implications for association mapping in germplasm arrays varying in size and structure

Background

Considerations in applying association mapping (AM) to plant breeding are population structure and size: not accounting for structure and/or using small populations can lead to elevated false-positive rates. The principal determinants of population structure in cultivated barley are growth habit and inflorescence type. Both are under complex genetic control: growth habit is controlled by the epistatic interactions of several genes. For inflorescence type, multiple loss-of-function alleles in one gene lead to the same phenotype. We used these two traits as models for assessing the effectiveness of AM. This research was initiated using the CAP Core germplasm array (n = 102) assembled at the start of the Barley Coordinated Agricultural Project (CAP). This array was genotyped with 4,608 SNPs and we re-sequenced genes involved in morphology, growth and development. Larger arrays of breeding germplasm were subsequently genotyped and phenotyped under the auspices of the CAP project. This provided sets of 247 accessions phenotyped for growth habit and 2,473 accessions phenotyped for inflorescence type. Each of the larger populations was genotyped with 3,072 SNPs derived from the original set of 4,608.

Results

Significant associations with SNPs located in the vicinity of the loci involved in growth habit and inflorescence type were found in the CAP Core. Differentiation of true and spurious associations was not possible without a priori knowledge of the candidate genes, based on re-sequencing. The re-sequencing data were used to define allele types of the determinant genes based on functional polymorphisms. In a second round of association mapping, these synthetic markers based on allele types gave the most significant associations. When the synthetic markers were used as anchor points for analysis of interactions, we detected other known-function genes and candidate loci involved in the control of growth habit and inflorescence type. We then conducted association analyses - with SNP data only - in the larger germplasm arrays. For both vernalization sensitivity and inflorescence type, the most significant associations in the larger data sets were found with SNPs coincident with the synthetic markers used in the CAP Core and with SNPs detected via interaction analysis in the CAP Core.

Conclusions

Small and highly structured collections of germplasm, such as the CAP Core, are cost-effectively phenotyped and genotyped with high-throughput markers. They are also useful for characterizing allelic diversity at loci in germplasm of interest. Our results suggest that discovery-oriented exercises in AM in such small arrays may generate a large number of false-positives. However, if haplotypes in candidate genes are available, they may be used as anchors in an analysis of interactions to identify other candidate regions harboring genes determining target traits. Using larger germplasm arrays, genome regions where the principal genes determining vernalization sensitivity and row type are located were identified.

Transcriptional profiling of root-knot nematode induced feeding sites in cowpea (Vigna unguiculata L. Walp.) using a soybean genome array

BACKGROUND

The locus Rk confers resistance against several species of root-knot nematodes (Meloidogyne spp., RKN) in cowpea (Vigna unguiculata). Based on histological and reactive oxygen species (ROS) profiles, Rk confers a delayed but strong resistance mechanism without a hypersensitive reaction-mediated cell death process, which allows nematode development but blocks reproduction.

RESULTS

Responses to M. incognita infection in roots of resistant genotype CB46 and a susceptible near-isogenic line (null-Rk) were investigated using a soybean Affymetrix GeneChip expression array at 3 and 9 days post-inoculation (dpi). At 9 dpi 552 genes were differentially expressed in incompatible interactions (infected resistant tissue compared with non-infected resistant tissue) and 1,060 genes were differentially expressed in compatible interactions (infected susceptible tissue compared with non-infected susceptible tissue). At 3 dpi the differentially expressed genes were 746 for the incompatible and 623 for the compatible interactions. When expression between infected resistant and susceptible genotypes was compared, 638 and 197 genes were differentially expressed at 9 and 3 dpi, respectively.

CONCLUSIONS

In comparing the differentially expressed genes in response to nematode infection, a greater number and proportion of genes were down-regulated in the resistant than in the susceptible genotype, whereas more genes were up-regulated in the susceptible than in the resistant genotype. Gene ontology based functional categorization revealed that the typical defense response was partially suppressed in resistant roots, even at 9 dpi, allowing nematode juvenile development. Differences in ROS concentrations, induction of toxins and other defense related genes seem to play a role in this unique resistance mechanism.

Association mapping of spot blotch resistance in wild barley

Spot blotch, caused by Cochliobolus sativus, is an important foliar disease of barley. The disease has been controlled for over 40 years through the deployment of cultivars with durable resistance derived from the line NDB112. Pathotypes of C. sativus with virulence for the NDB112 resistance have been detected in Canada; thus, many commercial cultivars are vulnerable to spot blotch epidemics. To increase the diversity of spot blotch resistance in cultivated barley, we evaluated 318 diverse wild barley accessions comprising the Wild Barley Diversity Collection (WBDC) for reaction to C. sativus at the seedling stage and utilized an association mapping (AM) approach to identify and map resistance loci. A high frequency of resistance was found in the WBDC as 95% (302/318) of the accessions exhibited low infection responses. The WBDC was genotyped with 558 Diversity Array Technology (DArT((R))) and 2,878 single nucleotide polymorphism (SNP) markers and subjected to structure analysis before running the AM procedure. Thirteen QTL for spot blotch resistance were identified with DArT and SNP markers. These QTL were found on chromosomes 1H, 2H, 3H, 5H, and 7H and explained from 2.3 to 3.9% of the phenotypic variance. Nearly half of the identified QTL mapped to chromosome bins where spot blotch resistance loci were previously reported, offering some validation for the AM approach. The other QTL mapped to unique genomic regions and may represent new spot blotch resistance loci. This study demonstrates that AM is an effective technique for identifying and mapping QTL for disease resistance in a wild crop progenitor. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-010-9402-8) contains supplementary material, which is available to authorized users.

Response to zinc deficiency of two rice lines with contrasting tolerance is determined by root growth maintenance and organic acid exudation rates, and not by zinc-transporter activity

*Zinc (Zn)-deficient soils constrain rice (Oryza sativa) production and cause Zn malnutrition. The identification of Zn-deficiency-tolerant rice lines indicates that breeding might overcome these constraints. Here, we seek to identify processes underlying Zn-deficiency tolerance in rice at the physiological and transcriptional levels. *A Zn-deficiency-tolerant line RIL46 acquires Zn more efficiently and produces more biomass than its nontolerant maternal line (IR74) at low [Zn](ext) under field conditions. We tested if this was the result of increased expression of Zn(2+) transporters; increased root exudation of deoxymugineic acid (DMA) or low-molecular-weight organic acids (LMWOAs); and/or increased root production. Experiments were performed in field and controlled environment conditions. *There was little genotypic variation in transcript abundance of Zn-responsive root Zn(2+)-transporters between the RIL46 and IR74. However, root exudation of DMA and LMWOA was greater in RIL46, coinciding with increased root expression of putative ligand-efflux genes. Adventitious root production was maintained in RIL46 at low [Zn](ext), correlating with altered expression of root-specific auxin-responsive genes. *Zinc-deficiency tolerance in RIL46 is most likely the result of maintenance of root growth, increased efflux of Zn ligands, and increased uptake of Zn-ligand complexes at low [Zn](ext); these traits are potential breeding targets.

Development and implementation of high-throughput SNP genotyping in barley

BACKGROUND

High density genetic maps of plants have, nearly without exception, made use of marker datasets containing missing or questionable genotype calls derived from a variety of genic and non-genic or anonymous markers, and been presented as a single linear order of genetic loci for each linkage group. The consequences of missing or erroneous data include falsely separated markers, expansion of cM distances and incorrect marker order. These imperfections are amplified in consensus maps and problematic when fine resolution is critical including comparative genome analyses and map-based cloning. Here we provide a new paradigm, a high-density consensus genetic map of barley based only on complete and error-free datasets and genic markers, represented accurately by graphs and approximately by a best-fit linear order, and supported by a readily available SNP genotyping resource.

RESULTS

Approximately 22,000 SNPs were identified from barley ESTs and sequenced amplicons; 4,596 of them were tested for performance in three pilot phase Illumina GoldenGate assays. Data from three barley doubled haploid mapping populations supported the production of an initial consensus map. Over 200 germplasm selections, principally European and US breeding material, were used to estimate minor allele frequency (MAF) for each SNP. We selected 3,072 of these tested SNPs based on technical performance, map location, MAF and biological interest to fill two 1536-SNP "production" assays (BOPA1 and BOPA2), which were made available to the barley genetics community. Data were added using BOPA1 from a fourth mapping population to yield a consensus map containing 2,943 SNP loci in 975 marker bins covering a genetic distance of 1099 cM.

CONCLUSION

The unprecedented density of genic markers and marker bins enabled a high resolution comparison of the genomes of barley and rice. Low recombination in pericentric regions is evident from bins containing many more than the average number of markers, meaning that a large number of genes are recombinationally locked into the genetic centromeric regions of several barley chromosomes. Examination of US breeding germplasm illustrated the usefulness of BOPA1 and BOPA2 in that they provide excellent marker density and sensitivity for detection of minor alleles in this genetically narrow material.

Development and implementation of high-throughput SNP genotyping in barley

BACKGROUND

High density genetic maps of plants have, nearly without exception, made use of marker datasets containing missing or questionable genotype calls derived from a variety of genic and non-genic or anonymous markers, and been presented as a single linear order of genetic loci for each linkage group. The consequences of missing or erroneous data include falsely separated markers, expansion of cM distances and incorrect marker order. These imperfections are amplified in consensus maps and problematic when fine resolution is critical including comparative genome analyses and map-based cloning. Here we provide a new paradigm, a high-density consensus genetic map of barley based only on complete and error-free datasets and genic markers, represented accurately by graphs and approximately by a best-fit linear order, and supported by a readily available SNP genotyping resource.

RESULTS

Approximately 22,000 SNPs were identified from barley ESTs and sequenced amplicons; 4,596 of them were tested for performance in three pilot phase Illumina GoldenGate assays. Data from three barley doubled haploid mapping populations supported the production of an initial consensus map. Over 200 germplasm selections, principally European and US breeding material, were used to estimate minor allele frequency (MAF) for each SNP. We selected 3,072 of these tested SNPs based on technical performance, map location, MAF and biological interest to fill two 1536-SNP "production" assays (BOPA1 and BOPA2), which were made available to the barley genetics community. Data were added using BOPA1 from a fourth mapping population to yield a consensus map containing 2,943 SNP loci in 975 marker bins covering a genetic distance of 1099 cM.

CONCLUSION

The unprecedented density of genic markers and marker bins enabled a high resolution comparison of the genomes of barley and rice. Low recombination in pericentric regions is evident from bins containing many more than the average number of markers, meaning that a large number of genes are recombinationally locked into the genetic centromeric regions of several barley chromosomes. Examination of US breeding germplasm illustrated the usefulness of BOPA1 and BOPA2 in that they provide excellent marker density and sensitivity for detection of minor alleles in this genetically narrow material.

A consensus genetic map of cowpea [Vigna unguiculata (L) Walp.] and synteny based on EST-derived SNPs

Consensus genetic linkage maps provide a genomic framework for quantitative trait loci identification, map-based cloning, assessment of genetic diversity, association mapping, and applied breeding in marker-assisted selection schemes. Among "orphan crops" with limited genomic resources such as cowpea [Vigna unguiculata (L.) Walp.] (2n = 2x = 22), the use of transcript-derived SNPs in genetic maps provides opportunities for automated genotyping and estimation of genome structure based on synteny analysis. Here, we report the development and validation of a high-throughput EST-derived SNP assay for cowpea, its application in consensus map building, and determination of synteny to reference genomes. SNP mining from 183,118 ESTs sequenced from 17 cDNA libraries yielded approximately 10,000 high-confidence SNPs from which an Illumina 1,536-SNP GoldenGate genotyping array was developed and applied to 741 recombinant inbred lines from six mapping populations. Approximately 90% of the SNPs were technically successful, providing 1,375 dependable markers. Of these, 928 were incorporated into a consensus genetic map spanning 680 cM with 11 linkage groups and an average marker distance of 0.73 cM. Comparison of this cowpea genetic map to reference legumes, soybean (Glycine max) and Medicago truncatula, revealed extensive macrosynteny encompassing 85 and 82%, respectively, of the cowpea map. Regions of soybean genome duplication were evident relative to the simpler diploid cowpea. Comparison with Arabidopsis revealed extensive genomic rearrangement with some conserved microsynteny. These results support evolutionary closeness between cowpea and soybean and identify regions for synteny-based functional genomics studies in legumes.

Integrated genetic map and genetic analysis of a region associated with root traits on the short arm of rye chromosome 1 in bread wheat

A rye-wheat centric chromosome translocation 1RS.1BL has been widely used in wheat breeding programs around the world. Increased yield of translocation lines was probably a consequence of increased root biomass. In an effort to map loci-controlling root characteristics, homoeologous recombinants of 1RS with 1BS were used to generate a consensus genetic map comprised of 20 phenotypic and molecular markers, with an average spacing of 2.5 cM. Physically, all recombination events were located in the distal 40% of the arms. A total of 68 recombinants was used and recombination breakpoints were aligned and ordered over map intervals with all the markers, integrated together in a genetic map. This approach enabled dissection of genetic components of quantitative traits, such as root traits, present on 1S. To validate our hypothesis, phenotyping of 45-day-old wheat roots was performed in five lines including three recombinants representative of the entire short arm along with bread wheat parents 'Pavon 76' and Pavon 1RS.1BL. Individual root characteristics were ranked and the genotypic rank sums were subjected to Quade analysis to compare the overall rooting ability of the genotypes. It appears that the terminal 15% of the rye 1RS arm carries gene(s) for greater rooting ability in wheat.

Comparing genomic expression patterns across plant species reveals highly diverged transcriptional dynamics in response to salt stress

Background

Rice and barley are both members of Poaceae (grass family) but have a marked difference in salt tolerance. The molecular mechanism underlying this difference was previously unexplored. This study employs a comparative genomics approach to identify analogous and contrasting gene expression patterns between rice and barley.

Results

A hierarchical clustering approach identified several interesting expression trajectories among rice and barley genotypes. There were no major conserved expression patterns between the two species in response to salt stress. A wheat salt-stress dataset was queried for comparison with rice and barley. Roughly one-third of the salt-stress responses of barley were conserved with wheat while overlap between wheat and rice was minimal. These results demonstrate that, at transcriptome level, rice is strikingly different compared to the more closely related barley and wheat. This apparent lack of analogous transcriptional programs in response to salt stress is further highlighted through close examination of genes associated with root growth and development.

Conclusion

The analysis provides support for the hypothesis that conservation of transcriptional signatures in response to environmental cues depends on the genetic similarity among the genotypes within a species, and on the phylogenetic distance between the species.

Evidence and evolutionary analysis of ancient whole-genome duplication in barley predating the divergence from rice

BACKGROUND

Well preserved genomic colinearity among agronomically important grass species such as rice, maize, Sorghum, wheat and barley provides access to whole-genome structure information even in species lacking a reference genome sequence. We investigated footprints of whole-genome duplication (WGD) in barley that shaped the cereal ancestor genome by analyzing shared synteny with rice using a approximately 2000 gene-based barley genetic map and the rice genome reference sequence.

RESULTS

Based on a recent annotation of the rice genome, we reviewed the WGD in rice and identified 24 pairs of duplicated genomic segments involving 70% of the rice genome. Using 968 putative orthologous gene pairs, synteny covered 89% of the barley genetic map and 63% of the rice genome. We found strong evidence for seven shared segmental genome duplications, corresponding to more than 50% of the segmental genome duplications previously determined in rice. Analysis of synonymous substitution rates (Ks) suggested that shared duplications originated before the divergence of these two species. While major genome rearrangements affected the ancestral genome of both species, small paracentric inversions were found to be species specific.

CONCLUSION

We provide a thorough analysis of comparative genome evolution between barley and rice. A barley genetic map of approximately 2000 non-redundant EST sequences provided sufficient density to allow a detailed view of shared synteny with the rice genome. Using an indirect approach that included the localization of WGD-derived duplicated genome segments in the rice genome, we determined the current extent of shared WGD-derived genome duplications that occurred prior to species divergence.

A compartmentalized approach to the assembly of physical maps

BACKGROUND

Physical maps have been historically one of the cornerstones of genome sequencing and map-based cloning strategies. They also support marker assisted breeding and EST mapping. The problem of building a high quality physical map is computationally challenging due to unavoidable noise in the input fingerprint data.

RESULTS

We propose a novel compartmentalized method for the assembly of high quality physical maps from fingerprinted clones. The knowledge of genetic markers enables us to group clones into clusters so that clones in the same cluster are more likely to overlap. For each cluster of clones, a local physical map is first constructed using FingerPrinted Contigs (FPC). Then, all the individual maps are carefully merged into the final physical map. Experimental results on the genomes of rice and barley demonstrate that the compartmentalized assembly produces significantly more accurate maps, and that it can detect and isolate clones that would induce "chimeric" contigs if used in the final assembly.

CONCLUSION

The software is available for download at http://www.cs.ucr.edu/~sbozdag/assembler/

The K-segment of maize DHN1 mediates binding to anionic phospholipid vesicles and concomitant structural changes

Dehydrins (DHNs; late embryogenesis abundant D11 family) are a family of intrinsically unstructured plant proteins that accumulate in the late stages of seed development and in vegetative tissues subjected to water deficit, salinity, low temperature, or abscisic acid treatment. We demonstrated previously that maize (Zea mays) DHNs bind preferentially to anionic phospholipid vesicles; this binding is accompanied by an increase in alpha-helicity of the protein, and adoption of alpha-helicity can be induced by sodium dodecyl sulfate. All DHNs contain at least one "K-segment," a lysine-rich 15-amino acid consensus sequence. The K-segment is predicted to form a class A2 amphipathic alpha-helix, a structural element known to interact with membranes and proteins. Here, three K-segment deletion proteins of maize DHN1 were produced. Lipid vesicle-binding assays revealed that the K-segment is required for binding to anionic phospholipid vesicles, and adoption of alpha-helicity of the K-segment accounts for most of the conformational change of DHNs upon binding to anionic phospholipid vesicles or sodium dodecyl sulfate. The adoption of structure may help stabilize cellular components, including membranes, under stress conditions.

The K-segment of maize DHN1 mediates binding to anionic phospholipid vesicles and concomitant structural changes

Dehydrins (DHNs; late embryogenesis abundant D11 family) are a family of intrinsically unstructured plant proteins that accumulate in the late stages of seed development and in vegetative tissues subjected to water deficit, salinity, low temperature, or abscisic acid treatment. We demonstrated previously that maize (Zea mays) DHNs bind preferentially to anionic phospholipid vesicles; this binding is accompanied by an increase in alpha-helicity of the protein, and adoption of alpha-helicity can be induced by sodium dodecyl sulfate. All DHNs contain at least one "K-segment," a lysine-rich 15-amino acid consensus sequence. The K-segment is predicted to form a class A2 amphipathic alpha-helix, a structural element known to interact with membranes and proteins. Here, three K-segment deletion proteins of maize DHN1 were produced. Lipid vesicle-binding assays revealed that the K-segment is required for binding to anionic phospholipid vesicles, and adoption of alpha-helicity of the K-segment accounts for most of the conformational change of DHNs upon binding to anionic phospholipid vesicles or sodium dodecyl sulfate. The adoption of structure may help stabilize cellular components, including membranes, under stress conditions.

Detection and validation of single feature polymorphisms using RNA expression data from a rice genome array

BACKGROUND

A large number of genetic variations have been identified in rice. Such variations must in many cases control phenotypic differences in abiotic stress tolerance and other traits. A single feature polymorphism (SFP) is an oligonucleotide array-based polymorphism which can be used for identification of SNPs or insertion/deletions (INDELs) for high throughput genotyping and high density mapping. Here we applied SFP markers to a lingering question about the source of salt tolerance in a particular rice recombinant inbred line (RIL) derived from a salt tolerant and salt sensitive parent.

RESULTS

Expression data obtained by hybridizing RNA to an oligonucleotide array were analyzed using a statistical method called robustified projection pursuit (RPP). By applying the RPP method, a total of 1208 SFP probes were detected between two presumed parental genotypes (Pokkali and IR29) of a RIL population segregating for salt tolerance. We focused on the Saltol region, a major salt tolerance QTL. Analysis of FL478, a salt tolerant RIL, revealed a small (< 1 Mb) region carrying alleles from the presumed salt tolerant parent, flanked by alleles matching the salt sensitive parent IR29. Sequencing of putative SFP-containing amplicons from this region and other positions in the genome yielded a validation rate more than 95%.

CONCLUSION

Recombinant inbred line FL478 contains a small (< 1 Mb) segment from the salt tolerant parent in the Saltol region. The Affymetrix rice genome array provides a satisfactory platform for high resolution mapping in rice using RNA hybridization and the RPP method of SFP analysis.

Orphan legume crops enter the genomics era!

Many of the world's most important food legumes are grown in arid and semi-arid regions of Africa and Asia, where crop productivity is hampered by biotic and abiotic stresses. Until recently, these crops have also suffered from a dearth of genomic and molecular-genetic resources and thus were 'orphans' of the genome revolution. However, the community of legume researchers has begun a concerted effort to change this situation. The driving force is a series of international collaborations that benefit from recent advances in genome sequencing and genotyping technologies. The focus of these activities is the development of genome-scale data sets that can be used in high-throughput approaches to facilitate genomics-assisted breeding in these legumes.

Mapping QTL for drought stress-induced premature senescence and maturity in cowpea [Vigna unguiculata (L.) Walp.]

Cowpea is an important crop for subsistence farmers in arid regions of Africa, Asia, and South America. Efforts to develop cultivars with improved productivity under drought conditions are constrained by lack of molecular markers associated with drought tolerance. Here, we report the mapping of 12 quantitative trait loci (QTL) associated with seedling drought tolerance and maturity in a cowpea recombinant inbred (RIL) population. One hundred and twenty-seven F(8) RILs developed from a cross between IT93K503-1 and CB46 were screened with 62 EcoR1 and Mse1 primer combinations to generate 306 amplified fragment length polymorphisms for use in genetic linkage mapping. The same population was phenotyped for maintenance of stem greenness (stg) and recovery dry weight (rdw) after drought stress in six greenhouse experiments. In field experiments conducted over 3 years, visual ratings and dry weights were used to phenotype drought stress-induced premature senescence in the RIL population. Kruskall-Wallis and multiple-QTL model mapping analysis were used to identify QTL associated with drought response phenotypes. Observed QTL were highly reproducible between stg and rdw under greenhouse conditions. Field studies confirmed all ten drought-response QTL observed under greenhouse conditions. Regions harboring drought-related QTL were observed on linkage groups 1, 2, 3, 5, 6, 7, 9, and 10 accounting for between 4.7 and 24.2% of the phenotypic variance (R(2)). Further, two QTL for maturity (R(2) = 14.4-28.9% and R(2) = 11.7-25.2%) mapped on linkage groups 7 and 8 separately from drought-related QTL. These results provide a platform for identification of genetic determinants of seedling drought tolerance in cowpea.

Mapping QTL for drought stress-induced premature senescence and maturity in cowpea [Vigna unguiculata (L.) Walp.]

Cowpea is an important crop for subsistence farmers in arid regions of Africa, Asia, and South America. Efforts to develop cultivars with improved productivity under drought conditions are constrained by lack of molecular markers associated with drought tolerance. Here, we report the mapping of 12 quantitative trait loci (QTL) associated with seedling drought tolerance and maturity in a cowpea recombinant inbred (RIL) population. One hundred and twenty-seven F(8) RILs developed from a cross between IT93K503-1 and CB46 were screened with 62 EcoR1 and Mse1 primer combinations to generate 306 amplified fragment length polymorphisms for use in genetic linkage mapping. The same population was phenotyped for maintenance of stem greenness (stg) and recovery dry weight (rdw) after drought stress in six greenhouse experiments. In field experiments conducted over 3 years, visual ratings and dry weights were used to phenotype drought stress-induced premature senescence in the RIL population. Kruskall-Wallis and multiple-QTL model mapping analysis were used to identify QTL associated with drought response phenotypes. Observed QTL were highly reproducible between stg and rdw under greenhouse conditions. Field studies confirmed all ten drought-response QTL observed under greenhouse conditions. Regions harboring drought-related QTL were observed on linkage groups 1, 2, 3, 5, 6, 7, 9, and 10 accounting for between 4.7 and 24.2% of the phenotypic variance (R(2)). Further, two QTL for maturity (R(2) = 14.4-28.9% and R(2) = 11.7-25.2%) mapped on linkage groups 7 and 8 separately from drought-related QTL. These results provide a platform for identification of genetic determinants of seedling drought tolerance in cowpea.

Towards systems genetic analyses in barley: Integration of phenotypic, expression and genotype data into GeneNetwork

BACKGROUND

A typical genetical genomics experiment results in four separate data sets; genotype, gene expression, higher-order phenotypic data and metadata that describe the protocols, processing and the array platform. Used in concert, these data sets provide the opportunity to perform genetic analysis at a systems level. Their predictive power is largely determined by the gene expression dataset where tens of millions of data points can be generated using currently available mRNA profiling technologies. Such large, multidimensional data sets often have value beyond that extracted during their initial analysis and interpretation, particularly if conducted on widely distributed reference genetic materials. Besides quality and scale, access to the data is of primary importance as accessibility potentially allows the extraction of considerable added value from the same primary dataset by the wider research community. Although the number of genetical genomics experiments in different plant species is rapidly increasing, none to date has been presented in a form that allows quick and efficient on-line testing for possible associations between genes, loci and traits of interest by an entire research community.

DESCRIPTION

Using a reference population of 150 recombinant doubled haploid barley lines we generated novel phenotypic, mRNA abundance and SNP-based genotyping data sets, added them to a considerable volume of legacy trait data and entered them into the GeneNetwork http://www.genenetwork.org. GeneNetwork is a unified on-line analytical environment that enables the user to test genetic hypotheses about how component traits, such as mRNA abundance, may interact to condition more complex biological phenotypes (higher-order traits). Here we describe these barley data sets and demonstrate some of the functionalities GeneNetwork provides as an easily accessible and integrated analytical environment for exploring them.

CONCLUSION

By integrating barley genotypic, phenotypic and mRNA abundance data sets directly within GeneNetwork's analytical environment we provide simple web access to the data for the research community. In this environment, a combination of correlation analysis and linkage mapping provides the potential to identify and substantiate gene targets for saturation mapping and positional cloning. By integrating datasets from an unsequenced crop plant (barley) in a database that has been designed for an animal model species (mouse) with a well established genome sequence, we prove the importance of the concept and practice of modular development and interoperability of software engineering for biological data sets.

Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph

Genetic linkage maps are cornerstones of a wide spectrum of biotechnology applications, including map-assisted breeding, association genetics, and map-assisted gene cloning. During the past several years, the adoption of high-throughput genotyping technologies has been paralleled by a substantial increase in the density and diversity of genetic markers. New genetic mapping algorithms are needed in order to efficiently process these large datasets and accurately construct high-density genetic maps. In this paper, we introduce a novel algorithm to order markers on a genetic linkage map. Our method is based on a simple yet fundamental mathematical property that we prove under rather general assumptions. The validity of this property allows one to determine efficiently the correct order of markers by computing the minimum spanning tree of an associated graph. Our empirical studies obtained on genotyping data for three mapping populations of barley (Hordeum vulgare), as well as extensive simulations on synthetic data, show that our algorithm consistently outperforms the best available methods in the literature, particularly when the input data are noisy or incomplete. The software implementing our algorithm is available in the public domain as a web tool under the name MSTmap.

Genetic linkage maps are cornerstones of a wide spectrum of biotechnology applications. In recent years, new high-throughput genotyping technologies have substantially increased the density and diversity of genetic markers, creating new algorithmic challenges for computational biologists. In this paper, we present a novel algorithmic method to construct genetic maps based on a new theoretical insight. Our approach outperforms the best methods available in the scientific literature, particularly when the input data are noisy or incomplete.

Deconvoluting BAC-gene relationships using a physical map

Deconvolution of relationships between bacterial artificial chromosome (BAC) clones and genes is a crucial step in the selective sequencing of regions of interest in a genome. It often includes combinatorial pooling of unique probes obtained from the genes (unigenes), and screening of the BAC library using the pools in a hybridization experiment. Since several probes can hybridize to the same BAC, in order for the deconvolution to be achievable the pooling design has to be able to handle a large number of positives. As a consequence, smaller pools need to be designed, which in turn increases the number of hybridization experiments, possibly making the entire protocol unfeasible. We propose a new algorithm that is capable of producing high-accuracy deconvolution even in the presence of a weak pooling design, i.e. when pools are rather large. The algorithm compensates for the decrease of information in the hybridization data by taking advantage of a physical map of the BAC clones. We show that the right combination of combinatorial pooling and our algorithm not only dramatically reduces the number of pools required, but also successfully deconvolutes the BAC-gene relationships with almost perfect accuracy. Software is available on request from the first author.

Exploiting regulatory variation to identify genes underlying quantitative resistance to the wheat stem rust pathogen Puccinia graminis f. sp. tritici in barley

We previously mapped mRNA transcript abundance traits (expression-QTL or eQTL) using the Barley1 Affymetrix array and 'whole plant' tissue from 139 progeny of the Steptoe x Morex (St/Mx) reference barley mapping population. Of the 22,840 probesets (genes) on the array, 15,987 reported transcript abundance signals that were suitable for eQTL analysis, and this revealed a genome-wide distribution of 23,738 significant eQTLs. Here we have explored the potential of using these mRNA abundance eQTL traits as surrogates for the identification of candidate genes underlying the interaction between barley and the wheat stem rust fungus Puccinia graminis f. sp. tritici. We re-analysed quantitative 'resistance phenotype' data collected on this population in 1990/1991 and identified six loci associated with barley's reaction to stem rust. One of these coincided with the major stem rust resistance locus Rpg1, that we had previously positionally cloned using this population. Correlation analysis between phenotype values for rust infection and mRNA abundance values reported by the 22,840 GeneChip probe sets placed Rpg1, which is on the Barley1 GeneChip, in the top five candidate genes for the major QTL on chromosome 7H corresponding to the location of Rpg1. A second co-located with the rpg4/Rpg5 stem rust resistance locus that has been mapped in a different population and the remaining four were novel. Correlation analyses identified candidate genes for the rpg4/Rpg5 locus on chromosome 5H. By combining our data with additional published mRNA profiling data sets, we identify a putative sensory transduction histidine kinase as a strong candidate for a novel resistance locus on chromosome 2H and compile candidate gene lists for the other three loci.

Coupling amplified DNA from flow-sorted chromosomes to high-density SNP mapping in barley

Background

Flow cytometry facilitates sorting of single chromosomes and chromosome arms which can be used for targeted genome analysis. However, the recovery of microgram amounts of DNA needed for some assays requires sorting of millions of chromosomes which is laborious and time consuming. Yet, many genomic applications such as development of genetic maps or physical mapping do not require large DNA fragments. In such cases time-consuming de novo sorting can be minimized by utilizing whole-genome amplification.

Results

Here we report a protocol optimized in barley including amplification of DNA from only ten thousand chromosomes, which can be isolated in less than one hour. Flow-sorted chromosomes were treated with proteinase K and amplified using Phi29 multiple displacement amplification (MDA). Overnight amplification in a 20-microlitre reaction produced 3.7 – 5.7 micrograms DNA with a majority of products between 5 and 30 kb. To determine the purity of sorted fractions and potential amplification bias we used quantitative PCR for specific genes on each chromosome. To extend the analysis to a whole genome level we performed an oligonucleotide pool assay (OPA) for interrogation of 1524 loci, of which 1153 loci had known genetic map positions. Analysis of unamplified genomic DNA of barley cv. Akcent using this OPA resulted in 1426 markers with present calls. Comparison with three replicates of amplified genomic DNA revealed >99% concordance. DNA samples from amplified chromosome 1H and a fraction containing chromosomes 2H – 7H were examined. In addition to loci with known map positions, 349 loci with unknown map positions were included. Based on this analysis 40 new loci were mapped to 1H.

Conclusion

The results indicate a significant potential of using this approach for physical mapping. Moreover, the study showed that multiple displacement amplification of flow-sorted chromosomes is highly efficient and representative which considerably expands the potential of chromosome flow sorting in plant genomics.

Coupling amplified DNA from flow-sorted chromosomes to high-density SNP mapping in barley

BACKGROUND

Flow cytometry facilitates sorting of single chromosomes and chromosome arms which can be used for targeted genome analysis. However, the recovery of microgram amounts of DNA needed for some assays requires sorting of millions of chromosomes which is laborious and time consuming. Yet, many genomic applications such as development of genetic maps or physical mapping do not require large DNA fragments. In such cases time-consuming de novo sorting can be minimized by utilizing whole-genome amplification.

RESULTS

Here we report a protocol optimized in barley including amplification of DNA from only ten thousand chromosomes, which can be isolated in less than one hour. Flow-sorted chromosomes were treated with proteinase K and amplified using Phi29 multiple displacement amplification (MDA). Overnight amplification in a 20-microlitre reaction produced 3.7 - 5.7 micrograms DNA with a majority of products between 5 and 30 kb. To determine the purity of sorted fractions and potential amplification bias we used quantitative PCR for specific genes on each chromosome. To extend the analysis to a whole genome level we performed an oligonucleotide pool assay (OPA) for interrogation of 1524 loci, of which 1153 loci had known genetic map positions. Analysis of unamplified genomic DNA of barley cv. Akcent using this OPA resulted in 1426 markers with present calls. Comparison with three replicates of amplified genomic DNA revealed >99% concordance. DNA samples from amplified chromosome 1H and a fraction containing chromosomes 2H - 7H were examined. In addition to loci with known map positions, 349 loci with unknown map positions were included. Based on this analysis 40 new loci were mapped to 1H.

CONCLUSION

The results indicate a significant potential of using this approach for physical mapping. Moreover, the study showed that multiple displacement amplification of flow-sorted chromosomes is highly efficient and representative which considerably expands the potential of chromosome flow sorting in plant genomics.

Dehydrin gene expression provides an indicator of low temperature and drought stress: transcriptome-based analysis of barley (Hordeum vulgare L.)

Low temperature and drought have major influences on plant growth and productivity. To identify barley genes involved in responses to these stresses and to specifically test the hypothesis that the dehydrin (Dhn) multigene family can serve as an indicator of the entire transcriptome response, we investigated the response of barley cv. Morex to: (1) gradual drought over 21 days and (2) low temperature including chilling, freeze-thaw cycles, and deacclimation over 33 days. We found 4,153 genes that responded to at least one component of these two stress regimes, about one fourth of all genes called "present" under any condition. About 44% (1,822 of 4,153) responded specifically to drought, whereas only 3.8% (158 of 4,153) were chilling specific and 2.8% (119 of 4,153) freeze-thaw specific, with 34.1% responsive to freeze-thaw and drought. The intersection between chilling and drought (31.9%) was somewhat smaller than the intersection between freeze-thaw and drought, implying an element of osmotic stress response to freeze-thaw. About 82.4% of the responsive genes were similar to Arabidopsis genes. The expression of 13 barley Dhn genes mirrored the global clustering of all transcripts, with specific combinations of Dhn genes providing an excellent indicator of each stress response. Data from these studies provide a robust reference data set for abiotic stress.

Barley grain maturation and germination: metabolic pathway and regulatory network commonalities and differences highlighted by new MapMan/PageMan profiling tools

Plant seeds prepare for germination already during seed maturation. We performed a detailed transcriptome analysis of barley (Hordeum vulgare) grain maturation, desiccation, and germination in two tissue fractions (starchy endosperm/aleurone and embryo/scutellum) using the Affymetrix Barley1 GeneChip. To aid data evaluation, Arabidopsis thaliana MapMan and PageMan tools were adapted to barley. The analyses allow a number of conclusions: (1) Cluster analysis revealed a smooth transition in transcription programs between late seed maturation and germination within embryo tissues, but not in the endosperm/aleurone. (2) More than 12,000 transcripts are stored in the embryo of dry barley grains, many of which are presumably activated during germination. (3) Transcriptional activation of storage reserve mobilization events occurs at an early stage of germination, well before radicle protrusion. (4) Key genes of gibberellin (GA) biosynthesis are already active during grain maturation at a time when abscisic acid peaks suggesting the formation of an endogenous store of GA in the aleurone. This GA probably acts later during germination in addition to newly synthesized GA. (5) Beside the well-known role of GA in gene activation during germination spatiotemporal expression data and cis-element searches in homologous rice promoters confirm an equally important gene-activating role of abscisic acid during this developmental period. The respective regulatory webs are linked to auxin and ethylene controlled networks. In summary, new bioinformatics PageMan and MapMan tools developed in barley have been successfully used to investigate in detail the transcriptome relationships between seed maturation and germination in an important crop plant.

Detection and validation of single feature polymorphisms in cowpea (Vigna unguiculata L. Walp) using a soybean genome array

BACKGROUND

Cowpea (Vigna unguiculata L. Walp) is an important food and fodder legume of the semiarid tropics and subtropics worldwide, especially in sub-Saharan Africa. High density genetic linkage maps are needed for marker assisted breeding but are not available for cowpea. A single feature polymorphism (SFP) is a microarray-based marker which can be used for high throughput genotyping and high density mapping.

RESULTS

Here we report detection and validation of SFPs in cowpea using a readily available soybean (Glycine max) genome array. Robustified projection pursuit (RPP) was used for statistical analysis using RNA as a surrogate for DNA. Using a 15% outlying score cut-off, 1058 potential SFPs were enumerated between two parents of a recombinant inbred line (RIL) population segregating for several important traits including drought tolerance, Fusarium and brown blotch resistance, grain size and photoperiod sensitivity. Sequencing of 25 putative polymorphism-containing amplicons yielded a SFP probe set validation rate of 68%.

CONCLUSION

We conclude that the Affymetrix soybean genome array is a satisfactory platform for identification of some 1000's of SFPs for cowpea. This study provides an example of extension of genomic resources from a well supported species to an orphan crop. Presumably, other legume systems are similarly tractable to SFP marker development using existing legume array resources.

On the accurate construction of consensus genetic maps

We study the problem of merging genetic maps, when the individual genetic maps are given as directed acyclic graphs. The problem is to build a consensus map, which includes and is consistent with all (or, the vast majority of) the markers in the individual maps. When markers in the input maps have ordering conflicts, the resulting consensus map will contain cycles. We formulate the problem of resolving cycles in a combinatorial optimization framework, which in turn is expressed as an integer linear program. A faster approximation algorithm is proposed, and an additional speed-up heuristic is developed. According to an extensive set of experimental results, our tool is consistently better than JOINMAP, both in terms of accuracy and running time.

Histological characterization of root-knot nematode resistance in cowpea and its relation to reactive oxygen species modulation

Root-knot nematodes (Meloidogyne spp.) are sedentary endoparasites with a broad host range which includes economically important crop species. Cowpea (Vigna unguiculata L. Walp) is an important food and fodder legume grown in many regions where root-knot nematodes are a major problem in production fields. Several sources of resistance to root-knot nematode have been identified in cowpea, including the widely used Rk gene. As part of a study to elucidate the mechanism of Rk-mediated resistance, the histological response to avirulent M. incognita feeding of a resistant cowpea cultivar CB46 was compared with a susceptible near-isogenic line (in CB46 background). Most root-knot nematode resistance mechanisms in host plants that have been examined induced a hypersensitive response (HR). However, there was no typical HR in resistant cowpea roots and nematodes were able to develop normal feeding sites similar to those in susceptible roots up to 9-14 d post inoculation (dpi). From 14-21 dpi giant cell deterioration was observed and the female nematodes showed arrested development and deterioration. Nematodes failed to reach maturity and did not initiate egg laying in resistant roots. These results confirmed that the induction of resistance is relatively late in this system. Typically in pathogen resistance HR is closely associated with an oxidative burst (OB) in infected tissue. The level of reactive oxygen species release in both compatible and incompatible reactions during early and late stages of infection was also quantified. Following a basal OB during early infection in both susceptible and resistant roots, which was also observed in mechanically wounded root tissues, no significant OB was detected up to 14 dpi, a profile consistent with the histological observations of a delayed resistance response. These results will be useful to design gene expression experiments to dissect Rk-mediated resistance at the molecular level.

Mapping translocation breakpoints using a wheat microarray

We report mapping of translocation breakpoints using a microarray. We used complex RNA to compare normal hexaploid wheat (17,000 Mb genome) to a ditelosomic stock missing the short arm of chromosome 1B (1BS) and wheat-rye translocations that replace portions of 1BS with rye 1RS. Transcripts detected by a probe set can come from all three Triticeae genomes in ABD hexaploid wheat, and sequences of homoeologous genes on 1AS, 1BS and 1DS often differ from each other. Absence or replacement of 1BS therefore must sometimes result in patterns within a probe set that deviate from hexaploid wheat. We termed these 'high variance probe sets' (HVPs) and examined the extent to which HVPs associated with 1BS aneuploidy are related to rice genes on syntenic rice chromosome 5 short arm (5S). We observed an enrichment of such probe sets to 15-20% of all HVPs, while 1BS represents approximately 2% of the total genome. In total 257 HVPs constitute wheat 1BS markers. Two wheat-rye translocations subdivided 1BS HVPs into three groups, allocating translocation breakpoints to narrow intervals defined by rice 5S coordinates. This approach could be extended to the entire wheat genome or any organism with suitable aneuploid or translocation stocks.

Large-scale expression profiling and physiological characterization of jasmonic acid-mediated adaptation of barley to salinity stress

Barley (Hordeum vulgare L.) is a salt-tolerant member of the Triticeae. Recent transcriptome studies on salinity stress response in barley revealed regulation of jasmonic acid (JA) biosynthesis and JA-responsive genes by salt stress. From that observation and several other physiological reports, it was hypothesized that JA is involved in the adaptation of barley to salt stress. Here we tested that hypothesis by applying JA to barley plants and observing the physiological responses and transcriptome changes. Photosynthetic and sodium ion accumulation responses were compared after (1) salinity stress, (2) JA treatment and (3) JA pre-treatment followed by salinity stress. The JA-pre-treated salt-stressed plants accumulated strikingly low levels of Na(+) in the shoot tissue compared with untreated salt-stressed plants after several days of exposure to stress. In addition, pre-treatment with JA partially alleviated photosynthetic inhibition caused by salinity stress. Expression profiling after a short-term exposure to salinity stress indicated a considerable overlap between genes regulated by salinity stress and JA application. Three JA-regulated genes, arginine decarboxylase, ribulose 1.5-bisphosphate carboxylase/oxygenase (Rubisco) activase and apoplastic invertase are possibly involved in salinity tolerance mediated by JA. This work provides a reference data set for further study of the role of JA in salinity tolerance in barley and other plants species.

Large-scale expression profiling and physiological characterization of jasmonic acid-mediated adaptation of barley to salinity stress

Barley (Hordeum vulgare L.) is a salt-tolerant member of the Triticeae. Recent transcriptome studies on salinity stress response in barley revealed regulation of jasmonic acid (JA) biosynthesis and JA-responsive genes by salt stress. From that observation and several other physiological reports, it was hypothesized that JA is involved in the adaptation of barley to salt stress. Here we tested that hypothesis by applying JA to barley plants and observing the physiological responses and transcriptome changes. Photosynthetic and sodium ion accumulation responses were compared after (1) salinity stress, (2) JA treatment and (3) JA pre-treatment followed by salinity stress. The JA-pre-treated salt-stressed plants accumulated strikingly low levels of Na(+) in the shoot tissue compared with untreated salt-stressed plants after several days of exposure to stress. In addition, pre-treatment with JA partially alleviated photosynthetic inhibition caused by salinity stress. Expression profiling after a short-term exposure to salinity stress indicated a considerable overlap between genes regulated by salinity stress and JA application. Three JA-regulated genes, arginine decarboxylase, ribulose 1.5-bisphosphate carboxylase/oxygenase (Rubisco) activase and apoplastic invertase are possibly involved in salinity tolerance mediated by JA. This work provides a reference data set for further study of the role of JA in salinity tolerance in barley and other plants species.

Array-based genotyping and expression analysis of barley cv. Maythorpe and Golden Promise

BACKGROUND

Golden Promise is a salt-tolerant spring barley closely related to Maythorpe. Salt tolerance in Golden Promise has been attributed to a single mutation at the Ari-e locus (on 5H) resulting from irradiation of Maythorpe. Golden Promise accumulates lower shoot Na+ compared to Maythorpe when growing under saline conditions. This study focused on elucidating the genetic basis and mechanisms involved in this difference.

RESULTS

The level of polymorphism between the two genotypes was explored using the Barley1 GeneChip for single feature polymorphisms (SFPs) and an oligonucleotide pool assay for single nucleotide polymorphisms (SNPs). Polymorphism analyses revealed three haplotype blocks spanning 6.4 cM on chromosome 1H, 23.7 cM on chromosome 4H and 3.0 cM on 5H. The Barley1 GeneChip was used to examine transcript abundance in different tissues and stages during development. Several genes within the polymorphic haplotype blocks were differentially regulated. Additionally, a more global difference in the jasmonic acid pathway regulation was detected between the two genotypes.

CONCLUSION

The results confirm that Golden Promise and Maythorpe are genetically very closely related but establish that they are not isogenic, as previously reported, due to three polymorphic haplotype blocks. Transcriptome analysis indicates that the response of the two genotypes to salinity stress is quite different. Additionally, the response to salinity stress in the roots and shoot tissue is strikingly different.

Array-based genotyping and expression analysis of barley cv. Maythorpe and Golden Promise

BACKGROUND

Golden Promise is a salt-tolerant spring barley closely related to Maythorpe. Salt tolerance in Golden Promise has been attributed to a single mutation at the Ari-e locus (on 5H) resulting from irradiation of Maythorpe. Golden Promise accumulates lower shoot Na+ compared to Maythorpe when growing under saline conditions. This study focused on elucidating the genetic basis and mechanisms involved in this difference.

RESULTS

The level of polymorphism between the two genotypes was explored using the Barley1 GeneChip for single feature polymorphisms (SFPs) and an oligonucleotide pool assay for single nucleotide polymorphisms (SNPs). Polymorphism analyses revealed three haplotype blocks spanning 6.4 cM on chromosome 1H, 23.7 cM on chromosome 4H and 3.0 cM on 5H. The Barley1 GeneChip was used to examine transcript abundance in different tissues and stages during development. Several genes within the polymorphic haplotype blocks were differentially regulated. Additionally, a more global difference in the jasmonic acid pathway regulation was detected between the two genotypes.

CONCLUSION

The results confirm that Golden Promise and Maythorpe are genetically very closely related but establish that they are not isogenic, as previously reported, due to three polymorphic haplotype blocks. Transcriptome analysis indicates that the response of the two genotypes to salinity stress is quite different. Additionally, the response to salinity stress in the roots and shoot tissue is strikingly different.

Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage

Rice yield is most sensitive to salinity stress imposed during the panicle initiation (PI) stage. In this study, we have focused on physiological and transcriptional responses of four rice genotypes exposed to salinity stress during PI. The genotypes selected included a pair of indicas (IR63731 and IR29) and a pair of japonica (Agami and M103) rice subspecies with contrasting salt tolerance. Physiological characterization showed that tolerant genotypes maintained a much lower shoot Na+ concentration relative to sensitive genotypes under salinity stress. Global gene expression analysis revealed a strikingly large number of genes which are induced by salinity stress in sensitive genotypes, IR29 and M103 relative to tolerant lines. We found 19 probe sets to be commonly induced in all four genotypes. We found several salinity modulated, ion homeostasis related genes from our analysis. We also studied the expression of SKC1, a cation transporter reported by others as a major source of variation in salt tolerance in rice. The transcript abundance of SKC1 did not change in response to salinity stress at PI stage in the shoot tissue of all four genotypes. However, we found the transcript abundance of SKC1 to be significantly higher in tolerant japonica Agami relative to sensitive japonica M103 under control and stressed conditions during PI stage.

Deconvoluting the BAC-gene relationships using a physical map

MOTIVATION

The deconvolution of the relationships between BAC clones and genes is a crucial step in the selective sequencing of the regions of interest in a genome. It usually requires combinatorial pooling of unique probes obtained from the genes (unigenes), and the screening of the BAC library using the pools in a hybridization experiment. Since several probes can hybridize to the same BAC, in order for the deconvolution to be achievable the pooling design has to be able to handle a large number of positives. As a consequence, smaller pools need to be designed which in turn increases the number of hybridization experiments possibly making the entire protocol unfeasible.

RESULTS

We propose a new algorithm that is capable of producing high accuracy deconvolution even in the presence of a weak pooling design, i.e., when pools are rather large. The algorithm compensates for the decrease of information in the hybridization data by taking advantage of a physical map of the BAC clones. We show that the right combination of combinatorial pooling and our algorithm not only dramatically reduces the number of pools required, but also successfully deconvolutes the BAC-gene relationships with almost perfect accuracy.

An improved method to identify BAC clones using pooled overgos

Hybridization using overgo probes is an established approach for screening arrayed bacterial artificial chromosome (BAC) libraries. We have improved the use of overgos by increasing the yield of positive clones using reduced levels of radioisotopes and enzyme. The strategy involves labeling with all four radiolabeled nucleotides in a hot pulse followed by a cold nucleotide chase and then extending the exposure time to compensate for reduced specific activity of the probes. The resulting cost savings and reduced human exposure to radiation make the use of highly pooled overgo probes a more attractive approach for screening of BAC libraries from organisms with large genomes.

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.

454 sequencing put to the test using the complex genome of barley

Background

During the past decade, Sanger sequencing has been used to completely sequence hundreds of microbial and a few higher eukaryote genomes. In recent years, a number of alternative technologies became available, among them adaptations of the pyrosequencing procedure (i.e. "454 sequencing"), promising a ~100-fold increase in throughput over Sanger technology – an advancement which is needed to make large and complex genomes more amenable to full genome sequencing at affordable costs. Although several studies have demonstrated its potential usefulness for sequencing small and compact microbial genomes, it was unclear how the new technology would perform in large and highly repetitive genomes such as those of wheat or barley.

Results

To study its performance in complex genomes, we used 454 technology to sequence four barley Bacterial Artificial Chromosome (BAC) clones and compared the results to those from ABI-Sanger sequencing. All gene containing regions were covered efficiently and at high quality with 454 sequencing whereas repetitive sequences were more problematic with 454 sequencing than with ABI-Sanger sequencing. 454 sequencing provided a much more even coverage of the BAC clones than ABI-Sanger sequencing, resulting in almost complete assembly of all genic sequences even at only 9 to 10-fold coverage. To obtain highly advanced working draft sequences for the BACs, we developed a strategy to assemble large parts of the BAC sequences by combining comparative genomics, detailed repeat analysis and use of low-quality reads from 454 sequencing. Additionally, we describe an approach of including small numbers of ABI-Sanger sequences to produce hybrid assemblies to partly compensate the short read length of 454 sequences.

Conclusion

Our data indicate that 454 pyrosequencing allows rapid and cost-effective sequencing of the gene-containing portions of large and complex genomes and that its combination with ABI-Sanger sequencing and targeted sequence analysis can result in large regions of high-quality finished genomic sequences.

Use of a large-scale Triticeae expressed sequence tag resource to reveal gene expression profiles in hexaploid wheat (Triticum aestivum L.)

The US Wheat Genome Project, funded by the National Science Foundation, developed the first large public Triticeae expressed sequence tag (EST) resource. Altogether, 116,272 ESTs were produced, comprising 100,674 5' ESTs and 15 598 3' ESTs. These ESTs were derived from 42 cDNA libraries, which were created from hexaploid bread wheat (Triticum aestivum L.) and its close relatives, including diploid wheat (T. monococcum L. and Aegilops speltoides L.), tetraploid wheat (T. turgidum L.), and rye (Secale cereale L.), using tissues collected from various stages of plant growth and development and under diverse regimes of abiotic and biotic stress treatments. ESTs were assembled into 18,876 contigs and 23,034 singletons, or 41,910 wheat unigenes. Over 90% of the contigs contained fewer than 10 EST members, implying that the ESTs represented a diverse selection of genes and that genes expressed at low and moderate to high levels were well sampled. Statistical methods were used to study the correlation of gene expression patterns, based on the ESTs clustered in the 1536 contigs that contained at least 10 5' EST members and thus representing the most abundant genes expressed in wheat. Analysis further identified genes in wheat that were significantly upregulated (p < 0.05) in tissues under various abiotic stresses when compared with control tissues. Though the function annotation cannot be assigned for many of these genes, it is likely that they play a role associated with the stress response. This study predicted the possible functionality for 4% of total wheat unigenes, which leaves the remaining 96% with their functional roles and expression patterns largely unknown. Nonetheless, the EST data generated in this project provide a diverse and rich source for gene discovery in wheat.

A branched-chain amino acid aminotransferase gene isolated from Hordeum vulgare is differentially regulated by drought stress

Differential display was used to isolate cDNA clones showing differential expression in response to ABA, drought and cold in barley seedling shoots. One drought-regulated cDNA clone (DD12) was further analyzed and found to encode a branched-chain amino acid aminotransferase (HvBCAT-1). A genomic clone was isolated by probing the Morex BAC library with the cDNA clone DD12 and the structure of Hvbcat-1 was elucidated. The coding region is interrupted by six introns and contains a predicted mitochondrial transit peptide. Hvbcat1 was mapped to chromosome 4H. A comparison was made to rice and Arabidopsis genes to identify conserved structural patterns. Complementation of a yeast (Saccharomyces cerevisiae) double knockout strain revealed that HvBCAT-1 can function as the mitochondrial (catabolic) BCATs in vivo. Transcript levels of Hvbcat-1, increased in response to drought stress. As the first enzyme in the branched-chain amino acid (BCAA) catabolic pathway, HvBCAT-1 might have a role in the degradation of BCAA. Degradation of BCAA could serve as a detoxification mechanism that maintains the pool of free branched-chain amino acids at low and non toxic levels, under drought stress conditions.

Transcriptome analysis of cold acclimation in barley albina and xantha mutants

Previously, we have shown that barley (Hordeum vulgare) plants carrying a mutation preventing chloroplast development are completely frost susceptible as well as impaired in the expression of several cold-regulated genes. Here we investigated the transcriptome of barley albina and xantha mutants and the corresponding wild type to assess the effect of the chloroplast on expression of cold-regulated genes. First, by comparing control wild type against cold-hardened wild-type plants 2,735 probe sets with statistically significant changes (P = 0.05; > or = 2-fold change) were identified. Expression of these wild-type cold-regulated genes was then analyzed in control and cold-hardened mutants. Only about 11% of the genes cold regulated in wild type were regulated to a similar extent in all genotypes (chloroplast-independent cold-regulated genes); this class includes many genes known to be under C-repeat binding factor control. C-repeat binding factor genes were also equally induced in mutants and wild-type plants. About 67% of wild-type cold-regulated genes were not regulated by cold in any mutant (chloroplast-dependent cold-regulated genes). We found that the lack of cold regulation in the mutants is due to the presence of signaling pathway(s) normally cold activated in wild type but constitutively active in the mutants, as well as to the disruption of low-temperature signaling pathway(s) due to the absence of active chloroplasts. We also found that photooxidative stress signaling pathway is constitutively active in the mutants. These results demonstrate the major role of the chloroplast in the control of the molecular adaptation to cold.

Expression analysis of barley (Hordeum vulgare L.) during salinity stress

Barley (Hordeum vulgare L.) is a salt-tolerant crop species with considerable economic importance in salinity-affected arid and semiarid regions of the world. In this work, barley cultivar Morex was used for transcriptional profiling during salinity stress using a microarray containing approximately 22,750 probe sets. The experiment was designed to target the early responses of genes to a salinity stress at seedling stage. We found a comparable number of probe sets up-regulated and down-regulated in response to salinity. The differentially expressed genes were broadly characterized using gene ontology and through expression-based hierarchical clustering to identify interesting features in the data. A prominent feature of the response to salinity was the induction of genes involved in jasmonic acid biosynthesis and genes known to respond to jasmonic acid treatment. A large number of abiotic stress (heat, drought, and low temperature) related genes were also found to be responsive to salinity stress. Our results also indicate osmoprotection to be an early response of barley under salinity stress. Additionally, we compared the results of our studies with two other reports characterizing gene expression of barley under salinity stress and found very few genes in common.

OligoSpawn: a software tool for the design of overgo probes from large unigene datasets

BACKGROUND

Expressed sequence tag (EST) datasets represent perhaps the largest collection of genetic information. ESTs can be exploited in a variety of biological experiments and analysis. Here we are interested in the design of overlapping oligonucleotide (overgo) probes from large unigene (EST-contigs) datasets.

RESULTS

OLIGOSPAWN is a suite of software tools that offers two complementary services, namely (1) the selection of "unique" oligos each of which appears in one unigene but does not occur (exactly or approximately) in any other and (2) the selection of "popular" oligos each of which occurs (exactly or approximately) in as many unigenes as possible. In this paper, we describe the functionalities of OLIGOSPAWN and the computational methods it employs, and we report on experimental results for the overgo probes designed with it.

CONCLUSION

The algorithms we designed are highly efficient and capable of processing unigene datasets of sizes on the order of several tens of Mb in a few hours on a regular PC. The software has been used to design overgo probes employed to screen a barley BAC library (Hordeum vulgare). OLIGOSPAWN is freely available at http://oligospawn.ucr.edu/.

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.

Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage

Rice (Oryza sativa), a salt-sensitive species, has considerable genetic variation for salt tolerance within the cultivated gene pool. Two indica rice genotypes, FL478, a recombinant inbred line derived from a population developed for salinity tolerance studies, and IR29, the sensitive parent of the population, were selected for this study. We used the Affymetrix rice genome array containing 55,515 probe sets to explore the transcriptome of the salt-tolerant and salt-sensitive genotypes under control and salinity-stressed conditions during vegetative growth. Response of the sensitive genotype IR29 is characterized by induction of a relatively large number of probe sets compared to tolerant FL478. Salinity stress induced a number of genes involved in the flavonoid biosynthesis pathway in IR29 but not in FL478. Cell wall-related genes were responsive in both genotypes, suggesting cell wall restructuring is a general adaptive mechanism during salinity stress, although the two genotypes also had some differences. Additionally, the expression of genes mapping to the Saltol region of chromosome 1 were examined in both genotypes. Single-feature polymorphism analysis of expression data revealed that IR29 was the source of the Saltol region in FL478, contrary to expectation. This study provides a genome-wide transcriptional analysis of two well-characterized, genetically related rice genotypes differing in salinity tolerance during a gradually imposed salinity stress under greenhouse conditions.

Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage

Rice (Oryza sativa), a salt-sensitive species, has considerable genetic variation for salt tolerance within the cultivated gene pool. Two indica rice genotypes, FL478, a recombinant inbred line derived from a population developed for salinity tolerance studies, and IR29, the sensitive parent of the population, were selected for this study. We used the Affymetrix rice genome array containing 55,515 probe sets to explore the transcriptome of the salt-tolerant and salt-sensitive genotypes under control and salinity-stressed conditions during vegetative growth. Response of the sensitive genotype IR29 is characterized by induction of a relatively large number of probe sets compared to tolerant FL478. Salinity stress induced a number of genes involved in the flavonoid biosynthesis pathway in IR29 but not in FL478. Cell wall-related genes were responsive in both genotypes, suggesting cell wall restructuring is a general adaptive mechanism during salinity stress, although the two genotypes also had some differences. Additionally, the expression of genes mapping to the Saltol region of chromosome 1 were examined in both genotypes. Single-feature polymorphism analysis of expression data revealed that IR29 was the source of the Saltol region in FL478, contrary to expectation. This study provides a genome-wide transcriptional analysis of two well-characterized, genetically related rice genotypes differing in salinity tolerance during a gradually imposed salinity stress under greenhouse conditions.