Detecting single-feature polymorphisms using oligonucleotide arrays and robustified projection pursuit

MOTIVATION

Genomic DNA was hybridized to oligonucleotide microarrays to identify single-feature polymorphisms (SFP) for Arabidopsis, which has a genome size of approximately 130 Mb. However, that method does not work well for organisms such as barley, with a much larger 5200 Mb genome. In the present study, we demonstrate SFP detection using a small number of replicate datasets and complex RNA as a surrogate for barley DNA. To identify single probes defining SFPs in the data, we developed a method using robustified projection pursuit (RPP). This method first evaluates, for each probe set, the overall differentiation of signal intensities between two genotypes and then measures the contribution of the individual probes within the probe set to the overall differentiation.

RESULTS

RNA from whole seedlings with and without dehydration stress provided 'present' calls for approximately 75% of probe sets. Using triplicated data, among the 5% of 'present' probe sets identified as most likely to contain at least one SFP probe, at least 80% are correctly predicted. This was determined by direct sequencing of PCR amplicons derived from barley genomic DNA. Using a 5 percentile cutoff, we defined 2007 SFP probes contained in 1684 probe sets by combining three parental genotype comparisons: Steptoe versus Morex, Morex versus Barke and Oregon Wolfe Barley Dominant versus Recessive.

AVAILABILITY

The algorithm is available upon request from the corresponding author.

CONTACT

xinping.cui@ucr.edu

SUPPLEMENTARY INFORMATION

http://faculty.ucr.edu/~xpcui.

Deletion mapping of homoeologous group 6-specific wheat expressed sequence tags

To localize wheat (Triticum aestivum L.) ESTs on chromosomes, 882 homoeologous group 6-specific ESTs were identified by physically mapping 7965 singletons from 37 cDNA libraries on 146 chromosome, arm, and sub-arm aneuploid and deletion stocks. The 882 ESTs were physically mapped to 25 regions (bins) flanked by 23 deletion breakpoints. Of the 5154 restriction fragments detected by 882 ESTs, 2043 (loci) were localized to group 6 chromosomes and 806 were mapped on other chromosome groups. The number of loci mapped was greatest on chromosome 6B and least on 6D. The 264 ESTs that detected orthologous loci on all three homoeologs using one restriction enzyme were used to construct a consensus physical map. The physical distribution of ESTs was uneven on chromosomes with a tendency toward higher densities in the distal halves of chromosome arms. About 43% of the wheat group 6 ESTs identified rice homologs upon comparisons of genome sequences. Fifty-eight percent of these ESTs were present on rice chromosome 2 and the remaining were on other rice chromosomes. Even within the group 6 bins, rice chromosomal blocks identified by 1-6 wheat ESTs were homologous to up to 11 rice chromosomes. These rice-block contigs were used to resolve the order of wheat ESTs within each bin.

Construction and evaluation of cDNA libraries for large-scale expressed sequence tag sequencing in wheat (Triticum aestivum L.)

A total of 37 original cDNA libraries and 9 derivative libraries enriched for rare sequences were produced from Chinese Spring wheat (Triticum aestivum L.), five other hexaploid wheat genotypes (Cheyenne, Brevor, TAM W101, BH1146, Butte 86), tetraploid durum wheat (T. turgidum L.), diploid wheat (T. monococcum L.), and two other diploid members of the grass tribe Triticeae (Aegilops speltoides Tausch and Secale cereale L.). The emphasis in the choice of plant materials for library construction was reproductive development subjected to environmental factors that ultimately affect grain quality and yield, but roots and other tissues were also included. Partial cDNA expressed sequence tags (ESTs) were examined by various measures to assess the quality of these libraries. All ESTs were processed to remove cloning system sequences and contaminants and then assembled using CAP3. Following these processing steps, this assembly yielded 101,107 sequences derived from 89,043 clones, which defined 16,740 contigs and 33,213 singletons, a total of 49,953 "unigenes." Analysis of the distribution of these unigenes among the libraries led to the conclusion that the enrichment methods were effective in reducing the most abundant unigenes and to the observation that the most diverse libraries were from tissues exposed to environmental stresses including heat, drought, salinity, or low temperature.

A 2500-locus bin map of wheat homoeologous group 5 provides insights on gene distribution and colinearity with rice

We constructed high-density deletion bin maps of wheat chromosomes 5A, 5B, and 5D, including 2338 loci mapped with 1052 EST probes and 217 previously mapped loci (total 2555 loci). This information was combined to construct a consensus chromosome bin map of group 5 including 24 bins. A relatively higher number of loci were mapped on chromosome 5B (38%) compared to 5A (34%) and 5D (28%). Differences in the levels of polymorphism among the three chromosomes were partially responsible for these differences. A higher number of duplicated loci was found on chromosome 5B (42%). Three times more loci were mapped on the long arms than on the short arms, and a significantly higher number of probes, loci, and duplicated loci were mapped on the distal halves than on the proximal halves of the chromosome arms. Good overall colinearity was observed among the three homoeologous group 5 chromosomes, except for the previously known 5AL/4AL translocation and a putative small pericentric inversion in chromosome 5A. Statistically significant colinearity was observed between low-copy-number ESTs from wheat homoeologous group 5 and rice chromosomes 12 (88 ESTs), 9 (72 ESTs), and 3 (84 ESTs).

A chromosome bin map of 16,000 expressed sequence tag loci and distribution of genes among the three genomes of polyploid wheat

Because of the huge size of the common wheat (Triticum aestivum L., 2n = 6x = 42, AABBDD) genome of 17,300 Mb, sequencing and mapping of the expressed portion is a logical first step for gene discovery. Here we report mapping of 7104 expressed sequence tag (EST) unigenes by Southern hybridization into a chromosome bin map using a set of wheat aneuploids and deletion stocks. Each EST detected a mean of 4.8 restriction fragments and 2.8 loci. More loci were mapped in the B genome (5774) than in the A (5173) or D (5146) genomes. The EST density was significantly higher for the D genome than for the A or B. In general, EST density increased relative to the physical distance from the centromere. The majority of EST-dense regions are in the distal parts of chromosomes. Most of the agronomically important genes are located in EST-dense regions. The chromosome bin map of ESTs is a unique resource for SNP analysis, comparative mapping, structural and functional analysis, and polyploid evolution, as well as providing a framework for constructing a sequence-ready, BAC-contig map of the wheat genome.

Barley Dhn13 encodes a KS-type dehydrin with constitutive and stress responsive expression

Dehydrins (DHNs) compose a family of intrinsically unstructured proteins that have high water solubility and accumulate during late seed development, low temperature or water deficit conditions, and are thought to play a protective role in freezing and drought tolerance in plants. Twelve Dhn genes were previously described in the barley genome. Here, we report an additional member of this multigene family, Dhn13. The Dhn13 gene is located in chromosome 4 near marker MWG634 and encodes a 107-amino acid KS-type DHN. Semi-quantitative reverse transcriptase PCR data indicated that Dhn13 is constitutively expressed in seedling tissues and embryos of developing seeds. Microarray data were consistent with these results and showed a considerable increase of Dhn13 transcripts when plants were subjected to chilling and freezing temperatures. The highest transcript levels where observed in anthers. The presence of ABRE, MYC, DRE, and POLLEN1LELAT52 regulatory elements in the putative Dhn13 promoter region is in agreement with expression data.

A chromosome bin map of 2148 expressed sequence tag loci of wheat homoeologous group 7

The objectives of this study were to develop a high-density chromosome bin map of homoeologous group 7 in hexaploid wheat (Triticum aestivum L.), to identify gene distribution in these chromosomes, and to perform comparative studies of wheat with rice and barley. We mapped 2148 loci from 919 EST clones onto group 7 chromosomes of wheat. In the majority of cases the numbers of loci were significantly lower in the centromeric regions and tended to increase in the distal regions. The level of duplicated loci in this group was 24% with most of these loci being localized toward the distal regions. One hundred nineteen EST probes that hybridized to three fragments and mapped to the three group 7 chromosomes were designated landmark probes and were used to construct a consensus homoeologous group 7 map. An additional 49 probes that mapped to 7AS, 7DS, and the ancestral translocated segment involving 7BS also were designated landmarks. Landmark probe orders and comparative maps of wheat, rice, and barley were produced on the basis of corresponding rice BAC/PAC and genetic markers that mapped on chromosomes 6 and 8 of rice. Identification of landmark ESTs and development of consensus maps may provide a framework of conserved coding regions predating the evolution of wheat genomes.

Development of an expressed sequence tag (EST) resource for wheat (Triticum aestivum L.): EST generation, unigene analysis, probe selection and bioinformatics for a 16,000-locus bin-delineated map

This report describes the rationale, approaches, organization, and resource development leading to a large-scale deletion bin map of the hexaploid (2n = 6x = 42) wheat genome (Triticum aestivum L.). Accompanying reports in this issue detail results from chromosome bin-mapping of expressed sequence tags (ESTs) representing genes onto the seven homoeologous chromosome groups and a global analysis of the entire mapped wheat EST data set. Among the resources developed were the first extensive public wheat EST collection (113,220 ESTs). Described are protocols for sequencing, sequence processing, EST nomenclature, and the assembly of ESTs into contigs. These contigs plus singletons (unassembled ESTs) were used for selection of distinct sequence motif unigenes. Selected ESTs were rearrayed, validated by 5' and 3' sequencing, and amplified for probing a series of wheat aneuploid and deletion stocks. Images and data for all Southern hybridizations were deposited in databases and were used by the coordinators for each of the seven homoeologous chromosome groups to validate the mapping results. Results from this project have established the foundation for future developments in wheat genomics.

Analysis of expressed sequence tag loci on wheat chromosome group 4

A total of 1918 loci, detected by the hybridization of 938 expressed sequence tag unigenes (ESTs) from 26 Triticeae cDNA libraries, were mapped to wheat (Triticum aestivum L.) homoeologous group 4 chromosomes using a set of deletion, ditelosomic, and nulli-tetrasomic lines. The 1918 EST loci were not distributed uniformly among the three group 4 chromosomes; 41, 28, and 31% mapped to chromosomes 4A, 4B, and 4D, respectively. This pattern is in contrast to the cumulative results of EST mapping in all homoeologous groups, as reported elsewhere, that found the highest proportion of loci mapped to the B genome. Sixty-five percent of these 1918 loci mapped to the long arms of homoeologous group 4 chromosomes, while 35% mapped to the short arms. The distal regions of chromosome arms showed higher numbers of loci than the proximal regions, with the exception of 4DL. This study confirmed the complex structure of chromosome 4A that contains two reciprocal translocations and two inversions, previously identified. An additional inversion in the centromeric region of 4A was revealed. A consensus map for homoeologous group 4 was developed from 119 ESTs unique to group 4. Forty-nine percent of these ESTs were found to be homoeologous to sequences on rice chromosome 3, 12% had matches with sequences on other rice chromosomes, and 39% had no matches with rice sequences at all. Limited homology (only 26 of the 119 consensus ESTs) was found between wheat ESTs on homoeologous group 4 and the Arabidopsis genome. Forty-two percent of the homoeologous group 4 ESTs could be classified into functional categories on the basis of blastX searches against all protein databases.

A 2600-locus chromosome bin map of wheat homoeologous group 2 reveals interstitial gene-rich islands and colinearity with rice

The complex hexaploid wheat genome offers many challenges for genomics research. Expressed sequence tags facilitate the analysis of gene-coding regions and provide a rich source of molecular markers for mapping and comparison with model organisms. The objectives of this study were to construct a high-density EST chromosome bin map of wheat homoeologous group 2 chromosomes to determine the distribution of ESTs, construct a consensus map of group 2 ESTs, investigate synteny, examine patterns of duplication, and assess the colinearity with rice of ESTs assigned to the group 2 consensus bin map. A total of 2600 loci generated from 1110 ESTs were mapped to group 2 chromosomes by Southern hybridization onto wheat aneuploid chromosome and deletion stocks. A consensus map was constructed of 552 ESTs mapping to more than one group 2 chromosome. Regions of high gene density in distal bins and low gene density in proximal bins were found. Two interstitial gene-rich islands flanked by relatively gene-poor regions on both the short and long arms and having good synteny with rice were discovered. The map locations of two ESTs indicated the possible presence of a small pericentric inversion on chromosome 2B. Wheat chromosome group 2 was shown to share syntenous blocks with rice chromosomes 4 and 7.

Efficient selection of unique and popular oligos for large EST databases

MOTIVATION

Expressed sequence tag (EST) databases have grown exponentially in recent years and now represent the largest collection of genetic sequences. An important application of these databases is that they contain information useful for the design of gene-specific oligonucleotides (or simply, oligos) that can be used in PCR primer design, microarray experiments and genomic library screening.

RESULTS

In this paper, we study two complementary problems concerning the selection of short oligos, e.g. 20-50 bases, from a large database of tens of thousands of ESTs: (i) selection of oligos each of which appears (exactly) in one unigene but does not appear (exactly or approximately) in any other unigene and (ii) selection of oligos that appear (exactly or approximately) in many unigenes. The first problem is called the unique oligo problem and has applications in PCR primer and microarray probe designs, and library screening for gene-rich clones. The second is called the popular oligo problem and is also useful in screening genomic libraries. We present an efficient algorithm to identify all unique oligos in the unigenes and an efficient heuristic algorithm to enumerate the most popular oligos. By taking into account the distribution of the frequencies of the words in the unigene database, the algorithms have been engineered carefully to achieve remarkable running times on regular PCs. Each of the algorithms takes only a couple of hours (on a 1.2 GHz CPU, 1 GB RAM machine) to run on a dataset 28 Mb of barley unigenes from the HarvEST database. We present simulation results on the synthetic data and a preliminary analysis of the barley unigene database.

AVAILABILITY

Available on request from the authors.

A new resource for cereal genomics: 22K barley GeneChip comes of age

In recent years, access to complete genomic sequences, coupled with rapidly accumulating data related to RNA and protein expression patterns, has made it possible to determine comprehensively how genes contribute to complex phenotypes. However, for major crop plants, publicly available, standard platforms for parallel expression analysis have been limited. We report the conception and design of the new publicly available, 22K Barley1 GeneChip probe array, a model for plants without a fully sequenced genome. Array content was derived from worldwide contribution of 350,000 high-quality ESTs from 84 cDNA libraries, in addition to 1,145 barley (Hordeum vulgare) gene sequences from the National Center for Biotechnology Information nonredundant database. Conserved sequences expressed in seedlings of wheat (Triticum aestivum), oat (Avena strigosa), rice (Oryza sativa), sorghum (Sorghum bicolor), and maize (Zea mays) were identified that will be valuable in the design of arrays across grasses. To enhance the usability of the data, BarleyBase, a MIAME-compliant, MySQL relational database, serves as a public repository for raw and normalized expression data from the Barley1 GeneChip probe array. Interconnecting links with PlantGDB and Gramene allow BarleyBase users to perform gene predictions using the 21,439 non-redundant Barley1 exemplar sequences or cross-species comparison at the genome level, respectively. We expect that this first generation array will accelerate hypothesis generation and gene discovery in disease defense pathways, responses to abiotic stresses, development, and evolutionary diversity in monocot plants.

Comparative DNA sequence analysis of wheat and rice genomes

The use of DNA sequence-based comparative genomics for evolutionary studies and for transferring information from model species to crop species has revolutionized molecular genetics and crop improvement strategies. This study compared 4485 expressed sequence tags (ESTs) that were physically mapped in wheat chromosome bins, to the public rice genome sequence data from 2251 ordered BAC/PAC clones using BLAST. A rice genome view of homologous wheat genome locations based on comparative sequence analysis revealed numerous chromosomal rearrangements that will significantly complicate the use of rice as a model for cross-species transfer of information in nonconserved regions.

The organization and rate of evolution of wheat genomes are correlated with recombination rates along chromosome arms

Genes detected by wheat expressed sequence tags (ESTs) were mapped into chromosome bins delineated by breakpoints of 159 overlapping deletions. These data were used to assess the organizational and evolutionary aspects of wheat genomes. Relative gene density and recombination rate increased with the relative distance of a bin from the centromere. Single-gene loci present once in the wheat genomes were found predominantly in the proximal, low-recombination regions, while multigene loci tended to be more frequent in distal, high-recombination regions. One-quarter of all gene motifs within wheat genomes were represented by two or more duplicated loci (paralogous sets). For 40 such sets, ancestral loci and loci derived from them by duplication were identified. Loci derived by duplication were most frequently located in distal, high-recombination chromosome regions whereas ancestral loci were most frequently located proximal to them. It is suggested that recombination has played a central role in the evolution of wheat genome structure and that gradients of recombination rates along chromosome arms promote more rapid rates of genome evolution in distal, high-recombination regions than in proximal, low-recombination regions.

The binding of maize DHN1 to lipid vesicles. Gain of structure and lipid specificity

Dehydrins (DHNs; late embryogenesis abundant D-11) are a family of plant proteins induced in response to abiotic stresses such as drought, low temperature, and salinity or during the late stages of embryogenesis. Spectral and thermal properties of these proteins in purified form suggest that they are "intrinsically unstructured." However, DHNs contain at least one copy of a consensus 15-amino acid sequence, the "K segment," which resembles a class A2 amphipathic alpha-helical, lipid-binding domain found in other proteins such as apolipoproteins and alpha-synuclein. The presence of the K segment raises the question of whether DHNs bind lipids, bilayers, or phospholipid vesicles. Here, we show that maize (Zea mays) DHN DHN1 can bind to lipid vesicles that contain acidic phospholipids. We also observe that DHN1 binds more favorably to vesicles of smaller diameter than to larger vesicles, and that the association of DHN1 with vesicles results in an apparent increase of alpha-helicity of the protein. Therefore, DHNs, and presumably somewhat similar plant stress proteins in the late embryogenesis abundant and cold-regulated classes may undergo function-related conformational changes at the water/membrane interface, perhaps related to the stabilization of vesicles or other endomembrane structures under stress conditions.

Identification and mapping of a putative stress response regulator gene in barley

Plants respond to environmental stress with a number of physiological and developmental changes. Water deficit is one of the major factors limiting plant growth and development and crop productivity. One response of plants to water deficit is accumulation of abscisic acid (ABA). An increase of ABA is responsible for the induction of many genes, presumably some of which contribute to drought tolerance. Analysis of gene expression in barley seedling shoots by differential display reverse transcriptase polymerase chain reaction (DDRT-PCR) led to the isolation of several drought-, cold- and ABA-induced partial cDNA fragments. Here we extensively characterize one of these cDNAs, designated DD6. First, a larger cDNA was extended from DD6 by 5'-RACE (rapid amplification of cDNA ends). Subsequently, the corresponding gene was isolated by screening a barley BAC library, and the sequences of the transcribed and flanking regions were determined. The deduced amino acid sequence has similarity to an Arabidopsis hypothetical protein and to a human and mouse DNA-binding protein. The corresponding gene, named Srg6 (stress-responsive gene), was mapped in a barley doubled haploid mapping population to chromosome 7H between markers ABC455 and salfp76, within a region that previously has been linked to osmotic adaptation in barley and other grass genomes.

Barley Cbf3 gene identification, expression pattern, and map location

Although cold and drought adaptation in cereals and other plants involve the induction of a large number of genes, inheritance studies in Triticeae (wheat [Triticum aestivum], barley [Hordeum vulgare], and rye [Secale cereale]) have revealed only a few major loci for frost or drought tolerance that are consistent across multiple genetic backgrounds and environments. One might imagine that these loci could encode highly conserved regulatory factors that have global effects on gene expression; therefore, genes encoding central regulators identified in other plants might be orthologs of these Triticeae stress tolerance genes. The CBF/DREB1 regulators, identified originally in Arabidopsis as key components of cold and drought regulation, merit this consideration. We constructed barley cDNA libraries, screened these libraries and a barley bacterial artificial chromosome library using rice (Oryza sativa) and barley Cbf probes, found orthologs of Arabidopsis CBF/DREB1 genes, and examined the expression and genetic map location of the barley Cbf3 gene, HvCbf3. HvCbf3 was induced by a chilling treatment. HvCbf3 is located on barley chromosome 5H between markers WG364b and saflp58 on the barley cv Dicktoo x barley cv Morex genetic linkage map. This position is some 40 to 50 cM proximal to the winter hardiness quantitative trait locus that includes the Vrn-1H gene, but may coincide with the wheat 5A Rcg1 locus, which governs the threshold temperature at which cor genes are induced. From this, it remains possible that HvCbf3 is the basis of a minor quantitative trait locus in some genetic backgrounds, though that possibility remains to be thoroughly explored.

An improved genetic linkage map for cowpea (Vigna unguiculata L.) combining AFLP, RFLP, RAPD, biochemical markers, and biological resistance traits

An improved genetic linkage map has been constructed for cowpea (Vigna unguiculata L. Walp.) based on the segregation of various molecular markers and biological resistance traits in a population of 94 recombinant inbred lines (RILs) derived from the cross between 'IT84S-2049' and '524B'. A set of 242 molecular markers, mostly amplified fragment length polymorphism (AFLP), linked to 17 biological resistance traits, resistance genes, and resistance gene analogs (RGAs) were scored for segregation within the parental and recombinant inbred lines. These data were used in conjunction with the 181 random amplified polymorphic DNA (RAPD), restriction fragment length polymorphism (RFLP), AFLP, and biochemical markers previously mapped to construct an integrated linkage map for cowpea. The new genetic map of cowpea consists of 11 linkage groups (LGs) spanning a total of 2670 cM, with an average distance of 6.43 cM between markers. Astonishingly, a large, contiguous portion of LG1 that had been undetected in previous mapping work was discovered. This region, spanning about 580 cM, is composed entirely of AFLP markers (54 in total). In addition to the construction of a new map, molecular markers associated with various biological resistance and (or) tolerance traits, resistance genes, and RGAs were also placed on the map, including markers for resistance to Striga gesnerioides races 1 and 3, CPMV, CPSMV, B1CMV, SBMV, Fusarium wilt, and root-knot nematodes. These markers will be useful for the development of tools for marker-assisted selection in cowpea breeding, as well as for subsequent map-based cloning of the various resistance genes.

Expression of the barley dehydrin multigene family and the development of freezing tolerance

Dehydrins (DHNs; LEA D11) are one of the typical families of plant proteins that accumulate in response to dehydration, low temperature, osmotic stress or treatment with abscisic acid (ABA), or during seed maturation. We previously found that three genes encoding low-molecular-weight DHNs (Dhn1, Dhn2 and Dhn9) map within a 15-cM region of barley chromosome 5H that overlaps a QTL for winterhardiness, while other Dhn genes encoding low- and high-molecular-weight DHNs are located on chromosomes 3H, 4H and 6H. Here we examine the expression of specific Dhn genes under conditions associated with expression of the winterhardiness phenotype. Plants grown at 4 degrees C or in the field in Riverside, California developed similar, modest levels of freezing tolerance, coinciding with little low-MW Dhn gene activity. Dicktoo (the more tolerant cultivar) and Morex (the less tolerant) grown in Saskatoon, Canada expressed higher levels of expression of genes for low-MW DHNs than did the same cultivars in Riverside, with expression being higher in Dicktoo than Morex. Dehydration or freeze-thaw also evoked expression of genes for low MW DHNs, suggesting that the dehydration component of freeze-thaw in the field induces low expression of genes encoding low-MW DHNs. These observations are consistent with the hypothesis that the major chilling-induced DHNs help to prime plant cells for acclimation to more intense cold, which then involves adaptation to dehydration during freeze-thaw cycling. A role for chromosome 5H-encoded DHNs in acclimation to more intense cold seems possible, even though it is not the basis of the major heritable variation in winterhardiness within the Dicktoo x Morex population.

Allelic variation of a dehydrin gene cosegregates with chilling tolerance during seedling emergence

Dehydrins (DHNs, LEA D-11) are plant proteins present during environmental stresses associated with dehydration or low temperatures and during seed maturation. Functions of DHNs have not yet been defined. Earlier, we hypothesized that a approximately 35-kDa DHN and membrane properties that reduce electrolyte leakage from seeds confer chilling tolerance during seedling emergence of cowpea (Vigna unguiculata L. Walp.) in an additive and independent manner. Evidence for this hypothesis was not rigorous because it was based on correlations of presence/absence of the DHN and slow electrolyte leakage with chilling tolerance in closely related cowpea lines that have some other genetic differences. Here, we provide more compelling genetic evidence for involvement of the DHN in chilling tolerance of cowpea. We developed near-isogenic lines by backcrossing. We isolated and determined the sequence of a cDNA corresponding to the approximately 35-kDa DHN and used gene-specific oligonucleotides derived from it to test the genetic linkage between the DHN presence/absence trait and the DHN structural gene. We tested for association between the DHN presence/absence trait and both low-temperature seed emergence and electrolyte leakage. We show that allelic differences in the Dhn structural gene map to the same position as the DHN protein presence/absence trait and that the presence of the approximately 35-kDa DHN is indeed associated with chilling tolerance during seedling emergence, independent of electrolyte leakage effects. Two types of allelic variation in the Dhn gene were identified in the protein-coding region, deletion of one Phi-segment from the DHN-negative lines and two single amino acid substitutions.

Purification and partial characterization of a dehydrin involved in chilling tolerance during seedling emergence of cowpea

Dehydrins are a family of proteins (LEA [late-embryogenesis abundant] D11) commonly induced by environmental stresses associated with low temperature or dehydration and during seed maturation drying. Our previous genetic studies suggested an association of an approximately 35-kD protein (by immunological evidence a dehydrin) with chilling tolerance during emergence of seedlings of cowpea (Vigna unguiculata) line 1393-2-11. In the present study we found that the accumulation of this protein in developing cowpea seeds is coordinated with the start of the dehydration phase of embryo development. We purified this protein from dry seeds of cowpea line 1393-2-11 by using the characteristic high-temperature solubility of dehydrins as an initial enrichment step, which was followed by three chromatography steps involving cation exchange, hydrophobic interaction, and anion exchange. Various characteristics of this protein confirmed that indeed it is a dehydrin, including total amino acid composition, partial amino acid sequencing, and the adoption of alpha-helical structure in the presence of sodium dodecyl sulfate. The propensity of dehydrins to adopt alpha-helical structure in the presence of sodium dodecyl sulfate, together with the apparent polypeptide adhesion property of this cowpea dehydrin, suggests a role in stabilizing other proteins or membranes. Taken together, the genetic, physiological, and physicochemical data are at this stage consistent with a cause-and-effect relationship between the presence in mature seeds of the approximately 35-kD dehydrin, which is the product of a single member of a multigene family, and an increment of chilling tolerance during emergence of cowpea seedlings.

A ca. 40 kDa maize (Zea mays L.) embryo dehydrin is encoded by the dhn2 locus on chromosome 9

Dehydrins (LEA D11 proteins) are the products of multigene families in a number of higher plants. To date, however, only one dehydrin locus, dhn1 (a major embryo and drought-induced protein of ca. 18 kDa) has been placed on chromosome 6L of the genetic linkage map of maize. The presence of a larger, ca. 40 kDa embryo protein that is also specifically detected by anti-dehydrin antibodies had been observed in some maize inbreds, including B73, suggesting that other dhn loci may exist. The ca. 22 kDa and ca. 40 kDa immunopositive proteins were purified from B73 and their amino acid compositions determined. The two proteins' amino acid compositions are typical of dehydrins, yet they differ from each other, indicating that they are distinct dhn gene products. Different size alleles for both proteins, or presence/absence in the case of the ca. 40 kDa protein, were evident from comparisons of embryo proteins of various maize inbreds. Analysis of segregating F2 progeny derived from self-pollination of F1 hybrids from four crosses (B73 x OH43, Mo17 x A632, AHO x A632, Latente x A632) revealed that alleles of the two genes assort independently. Map positions of the two dhn loci were then determined using two maize recombinant inbred line (RIL) mapping populations. The predicted map position of the gene controlling production of the ca. 22 kDa protein confirmed that this protein is the product of the dhn1 gene. The gene encoding the ca. 40 kDa dehydrin-like protein maps to a new locus on chromosome 9S near wx1, which we have named dhn2.

Purification of a maize dehydrin protein expressed in Escherichia coli

A maize dehydrin protein (Dhn1) containing 167 amino acids with a predicted molecular weight of 17.0 kDa was produced in the Escherichia coli overexpression strain BL21 (DE3)pLysS. Site-directed mutagenesis was used to construct a plasmid with a protein coding region corresponding exactly to the original cDNA. Protein production was induced by IPTG. Dhn1 was enriched from total soluble protein by heat-fractionation and centrifugation and then purified by sequential cation exchange and hydrophobic interaction chromatography. The purified protein was visualized by SDS-PAGE and immunoblot analysis using a polyclonal antibody to the dehydrin consensus region. Expression in E. coli resulted in approximately 1.2 mg of purified protein per liter of induced culture.

Purification of a maize dehydrin

A maize dehydrin with an apparent molecular weight of 20 kDa was purified from whole kernels of maize inbred line G50. Kernels were ground in a seed mill, stirred overnight in extraction buffer, and centrifuged to extract soluble proteins. The sample was heated to 89 degrees C and centrifuged to remove heat-insoluble proteins. The remaining soluble proteins were fractionated in a three-step chromatographic process. Following cation exchange, hydrophobic interaction, and gel filtration chromatography, pure dehydrin samples were obtained.

A view of plant dehydrins using antibodies specific to the carboxy terminal peptide

Dehydrins are characterized by the consensus KIKEKLPG amino acid sequence found near the carboxy terminus, and usually repeated from one to many times within the protein. A synthetic peptide containing this consensus sequence was used to produce specific antibodies that recognize dehydrins in a wide range of plants. This range covered two families of monocots, viz. Gramineae (Hordeum vulgare L., Triticum aestivum L., Zea mays L., Oryza sativa L.) and Liliaceae (Allium sativa L.), and five families of dicots, Malvaceae (Gossypium hirsutum L.), Solanaceae (Lycopersicon esculentum L.), Brassicaceae (Raphanus sativus L.), Fabaceae (Vigna unguiculata L.), and Cucurbitaceae (Cucumis sativus L.). Two families of gymnosperms, Pinaceae (Pinus edulis Engelm.) and Ginkgoaceae (Ginkgo biloba L.), were also included. For several plants in which dehydrin cDNA and genomic clones have previously been characterized, it now appears that the dehydrin family of proteins is larger, and the regulation of dehydrin expression much more complex, than earlier studies have shown.

An osmotic stress protein of cyanobacteria is immunologically related to plant dehydrins

Dehydrins are a family of desiccation proteins that were identified originally in plants (T.J. Close, A.A. Kortt, P.M. Chandler [1989] Plant Mol Biol 13: 95-108; G. Galau, T.J. Close [1992] Plant Physiol 98: 1523-1525). Dehydrins are characterized by the consensus amino acid sequence domain EKKGIMDKIKEKLPG found at or near the carboxy terminus; the core of this domain (KIKEKLPG) may be repeated from one to many times within the complete polypeptide. Dehydrins generally accumulate in plants in response to dehydration stress, regardless of whether the stimulus is evaporation, chilling, or a decrease in external osmotic potential. Polyclonal antibodies highly specific to the consensus carboxy terminus of plant dehydrins were used to search for dehydrins in cyanobacteria, many of which are known to survive desiccation. A 40-kD osmotic-stress-induced protein was identified in Anabaena sp. strain PCC 7120. The 40-kD protein was usually not detected in logarithmic cultures and was induced by shifting the growth medium to higher solute concentrations. Several solutes have inductive effects, including sucrose, sorbitol, and polyethylene glycol (PEG). Measurements of osmotic potential suggest that a shift of -0.5 MPa (sucrose and PEG) or -1.2 MPa (sorbitol) is sufficient to induce synthesis of the 40-kD protein. Glycerol, which is highly permeable, was not an inducer at -1.2 MPa (0.5 M), nor was the plant hormone abscisic acid. Induction appears to be evoked by a shift in osmotic potential approximately equal in absolute magnitude to the expected turgor pressure of bacterial cells in logarithmic phase growth. A dehydrin-like polypeptide was also identified among osmotically induced proteins from two other filamentous, heterocyst-forming cyano-bacteria. A 40-kD protein was observed in Calothrix sp. strain PCC 7601, and in Nostoc sp. strain Mac-R2, an osmotic-induced doublet at 39 and 40 kD was observed. From these data, it appears that cyanobacteria produce a dehydrin-like protein under osmotic stress.

Control of transient expression of chimaeric genes by gibberellic acid and abscisic acid in protoplasts prepared from mature barley aleurone layers

Gibberellic acid (GA3) and abscisic acid (ABA) control the transcription of alpha-amylase genes in barley aleurone cells. This control is likely to be exerted through cis-acting hormone-responsive elements in the promoter region of the gene. In order to further define these elements, we have developed procedures for obtaining transient expression of chimaeric genes in protoplasts prepared from mature barley aleurone layers. Constructs with heterologous constitutive promoters and with heterologous and homologous GA3- and ABA-regulated promoters were expressed specifically by these cells. This system would appear to offer great potential in gene regulation studies especially for hormonally regulated homologous genes. Functional analysis of a barley alpha-amylase gene has been performed using this system. A 2050 bp fragment from a high-pI alpha-amylase gene was fused to a reporter gene (GUS) and control of its expression was examined. Deletion analysis of this promoter fragment showed that major GA- and ABA-responsive elements occurred between 174 and 41 bp upstream from the transcription initiation site.

Nucleotide sequences of genes encoding the type II chloramphenicol acetyltransferases of Escherichia coli and Haemophilus influenzae, which are sensitive to inhibition by thiol-reactive reagents

Sensitivity of enzymes to inhibition by thiol-reactive reagents is often presented as evidence for the possible involvement of cysteine residues in substrate binding and catalysis or to highlight possible important differences in structure and mechanism between closely related enzymes. The primary phenotypic distinction between the enterobacterial type II chloramphenicol acetyltransferase (CATII; typified by the enzyme encoded by the incW transmissible plasmid pSa) and the CATI and CATIII variants is the greatly enhanced susceptibility of CATII to inactivation by thiol-specific modifying reagents. Determination of the nucleotide sequence of the gene, catII, present on pSa and that of a related determinant, catIIH, isolated from Haemophilus influenzae indicates that sensitivity to such reagents cannot be due to the presence of additional reactive cysteine residues in CATII. Comparative analysis of the inactivation of CATII and CATIII by 5,5'-dithiobis-(2-nitrobenzoic acid) (DTNB), 4,4'-dithiodipyridine (DTDP) and methyl methanethiosulphonate (MMTS) suggests that (i) inactivation occurs as a result of chemical modification of the same residue (Cys-31) in each enzyme, (ii) reagents that inactivate via a pseudo-first-order process (DTNB and DTDP) appear to bind with a greater affinity to CATII, and (iii) the intrinsic reactivity of Cys-31 in CATII greatly exceeds that of the corresponding residue in CATIII. The results lead to the conclusion that a striking difference in chemical reactivity of a unique and conserved thiol group between closely related enzyme variants may not be easily explained even when a high-resolution tertiary structure is available for one of them. Plausible explanations include more favourable access of reagents to Cys-31 in CATII or an enhanced reactivity of its thiol group imposed by the side chains of residues that are not in immediate contact with it.

A cDNA-based comparison of dehydration-induced proteins (dehydrins) in barley and corn

Several cDNAs related to an ABA-induced cDNA from barley aleurone were isolated from barley and corn seedlings that were undergoing dehydration. Four different barley polypeptides with sizes of 22.6, 16.2, 14.4 and 14.2 kDa and a single corn polypeptide with a size of 17.0 kDa were predicted from the nucleotide sequences of the cDNAs. These dehydration-induced proteins (dehydrins) are very similar to each other and to a previously identified rice protein induced by ABA and salt, and have at least some similarity to a previously identified cotton embryo protein. Each dehydrin is extremely hydrophilic, glycine-rich, cysteine- and tryptophan-free and contains repeated units in a conserved linear order. A lysine-rich repeating unit occurs twice in each protein, once at the carboxy terminus and once partway through the polypeptide, adjacent to a succession of serines. This repeating unit and the adjacent flanking run of serines are conserved with minimal variation among all dehydrins. Another repeating unit is flanked by the two copies of the lysine-rich unit, and varies in number from one to five copies. This latter repeating unit is less conserved than the former, varying even within a singly dehydrin. The messenger RNAs corresponding to each cDNA are abundant in dehydrating, but not in well-watered seedlings. The amino acid sequence of tryptic peptides from purified dehydration-induced proteins of corn established that the corn cDNAs correspond to a protein that is produced in abundance during the response of corn seedlings to dehydration.

High levels of double-stranded transferred DNA (T-DNA) processing from an intact nopaline Ti plasmid

To obtain bacterial-mediated oncogenic transformation of plants, the transferred DNA (T-DNA) of the tumor-inducing (Ti) plasmid of Agrobacterium tumefaciens is transferred to its plant host cells during infection. The initial phases of transformation involve the processing of the T-DNA in the bacterial cell after induction of the vir genes located on the Ti plasmid. The kinetics and conditions of this processing were examined and upon induction with acetosyringone up to 40% of the left and right borders of the T-DNA were cleaved. This cleavage was dependent upon virA, virG, and VirD and was rec-independent. Processed T-DNA was observed within 30 min after induction and was delayed by an increased concentration of phosphate in the induction medium. When DNA was isolated in the absence of protease treatment, the DNA fragment corresponding to the left side of the cut at both the left and right border region exhibited gel retardation, suggesting one or more "pilot" proteins may be involved in T-DNA transfer. Although the relative abundance of a processed product does not necessarily imply relative importance, the preponderance of double-stranded cleavage products suggests that double-stranded T-DNA should be considered as a possible intermediate in T-DNA transfer.

Regulation of the vir genes of Agrobacterium tumefaciens plasmid pTiC58

The virulence (vir) region of pTiC58 was screened for promoter activities by using gene fusions to a promoterless lux operon in the broad-host-range vector pUCD615. Active vir fragments contained the strongly acetosyringone-inducible promoters of virB, virC, virD, and virE and the weakly inducible promoters of virA and virG. Identical induction patterns were obtained with freshly sliced carrot disks, suggesting that an inducer is released after plant tissue is wounded. Optimal conditions for vir gene induction were pH 5.7 for 50 microM acetosyringone or sinapic acid. The induction of virB and virE by acetosyringone was strictly dependent on intact virA and virG loci. An increase in the copy number of virG resulted in a proportional, acetosyringone-independent increase in vir gene expression, and a further increase occurred only if an inducing compound and virA were present.

Dual control of Agrobacterium tumefaciens Ti plasmid virulence genes

The virulence genes of nopaline (pTiC58) and octopine (pTiA6NC) Ti plasmids are similarly affected by the Agrobacterium tumefaciens ros mutation. Of six vir region complementation groups (virA, virB, virG, virC, virD, and virE) examined by using fusions to reporter genes, the promoters of only two (virC and virD) responded to the ros mutation. For each promoter that was affected by ros, the level of expression of its associated genes was substantially elevated in the mutant. This increase was not influenced by Ti plasmid-encoded factors, and the mutation did not interfere with the induction of pTiC58 vir genes by phenolic compounds via the VirA/VirG regulatory control mechanism. The effects of the ros mutation and acetosyringone were cumulative for all vir promoters examined. The pleiotropic characteristics of the ros mutant include the complete absence of the major acidic capsular polysaccharide.

Molecular characterization of the virC genes of the Ti plasmid

The virC (formerly bak) complementation group of the nopaline-type Ti plasmid pTiC58 encodes two proteins, VirC1 and VirC2. According to the primary structure of the polypeptides predicted by the nucleotide sequence, VirC1 is composed of 231 amino acids with a total molecular mass of 25.5 kilodaltons, and VirC2 is composed of 202 amino acids with a molecular mass of 22.1 kilodaltons. The pTiC58 VirC1 and VirC2 polypeptides are equal in length to VirC1 and VirC2 of the octopine-type plasmid pTiA6NC. VirC1 proteins of pTiC58 and pTiA6NC are identical at 202 (87.4%) of the amino acid residues, and this homology is distributed fairly evenly throughout the protein. VirC2 identities occur at 142 residues (70.3%), but fall predominantly into two blocks of higher homology (84.6 and 78.5%) separated by a 41-residue segment of much lower homology (29.3%). Mutations in virC resulted in attenuated virulence on all hosts tested, the severity of attenuation varying markedly depending on the type of plant inoculated. For example, the attenuation was more pronounced on Kalanchoe than on sunflower or jimson weed. Virulence was restored to normal on all hosts by in-trans complementation with corresponding nonmutant DNA fragments of pTiC58 or of the octopine-type plasmid pTi15955. Two oligopeptides from within the predicted pTiC58 VirC1 polypeptide were synthesized and used to raise antibodies. These antibodies were used to detect the VirC1 product of both pTiC58 and pTi15955. In both cases, virC was expressed constitutively in the Agrobacterium tumefaciens ros mutant. The homology between virC genes of octopine- and nopaline-type Ti plasmids thus includes a conservation of genetic regulatory control mechanisms as well as considerable conservation of the primary structure of the protein products.

Regulation of Ti plasmid virulence genes by a chromosomal locus of Agrobacterium tumefaciens

We isolated a mutant strain of Agrobacterium tumefaciens, designated Ros, that has a pleiotropic phenotype which includes elevated levels of expression of certain genes in the virulence (Vir) region of tumor-inducing plasmid pTiC58. This mutant and others were isolated by fusing the promoter of the Vir bak gene to a promoterless gene (cat) that encodes chloramphenicol acetyltransferase and then selecting for increased expression of cat in A. tumefaciens. The ros mutation is chromosomal in nature and is characterized by a more-than-300-fold increase in the level of expression of bak and a 12-fold increase in the level of expression of an adjacent divergent operon containing the hdv genes, which are involved in some aspect of host specificity. The Ros mutant is fully virulent and is typified by its unusual colony morphology; the colonies have rough surfaces, uneven edges, and a pit in the center at 24 degrees C and vary markedly in appearance from one growth temperature to another. The Ros mutant is not able to form colonies at 12 degrees C, a temperature at which the wild-type strain forms colonies in 14 days. The ros mutation occurs spontaneously with a frequency of 5 X 10(-8). Genetic and biochemical evidence indicates that the product of the ros locus is a negative regulator of Ti plasmid genes and is related to undefined chromosomally encoded functions that are involved in the mutant phenotype.

Design and development of amplifiable broad-host-range cloning vectors: analysis of the vir region of Agrobacterium tumefaciens plasmid pTiC58

The construction of a set of new plasmids that are suitable as general cloning vectors in Escherichia coli and Agrobacterium tumefaciens is described. Plasmid pUCD2 is amplifiable in E. coli, replicates in a wide range of gram-negative hosts and contains a number of useful restriction endonuclear cleavage sites and antibiotic resistance genes. This includes unique sites for KpnI, SacI, SacII, PstI, ClaI, SalI, EcoRV, and PvuII and the genes for resistance to kanamycin, tetracycline, ampicillin, and spectinomycin/streptomycin. Derivatives of pUCD2 include pUCD4, which has a unique XbaI site and the cosmid pUCD5, which also contains a unique EcoRI site. Two smaller plasmids pUCD9P and pUCD9X, contain many of the same unique sites as pUCD2 and pUCD4, but carry only the pBR322 replication origin and therefore do not display the extensive host-range of pSa. These plasmids were used to isolate and manipulate fragments of the A. tumefaciens pTiC58 plasmid in both E. coli and A. tumefaciens. Fragments from the virulence (vir) region of pTiC58 inserted immediately upstream of the spectinomycin resistance gene of pUCD2 resulted in spectinomycin resistance levels that varied greatly depending on the particular fragment and its orientation of insertion. Using this property we find that a major portion of the vir region of pTiC58 is transcribed in A. tumefaciens and E. coli from left to right toward the T region.

Genetic complementation of Agrobacterium tumefaciens Ti plasmid mutants in the virulence region

Mutants with Tn5 insertions in the vir region of the Agrobacterium tumefaciens TiC58 plasmid are unable to form crown-gall tumors. Complementation tests of these vir region mutants were carried out by constructing merodiploids in a recombination-deficient strain. Each merodiploid possessed a mutant TiC58 plasmid and a recombinant plasmid containing either the homologous wild-type DNA region or the homologous region containing a second Tn5 insertion. The analysis identified six complementation groups. Mutations in one of these complementation groups were not complemented in trans and represent a cis-dominant locus. The mutation in one complementation group showed variation in host range.

M13 bacteriophage and pUC plasmids containing DNA inserts but still capable of beta-galactosidase alpha-complementation

A DNA fragment encoding the transposon Tn9 chloramphenicol acetyltransferase gene (cat) was inserted into M13 phage and pUC plasmid cloning vehicles. When the cat gene was inserted in the same orientation as the lacZ gene, two new polypeptides were produced. One polypeptide possessed chloramphenicol acetyltransferase activity, while the other expressed beta-galactosidase alpha-donor activity. Both new polypeptides were translated from a hybrid messenger RNA initiating from the lac promoter. These observations may help explain why not all inserts produce white plaques.

Isolation of the origin of replication of the IncW-group plasmid pSa

The origin of replication of the IncW plasmid pSa has been cloned and the function of this origin in Escherichia coli examined. A 1.9-kb region of DNA is required for efficient autonomous replication, and a 0.47-kb fragment within this region can initiate replication only in the presence of an autonomously replicating derivative of pSa. An Mr 35,000 protein (repA) is encoded adjacent to the origin and is required for efficient initiation of replication. The derivatives examined provide information suggesting a direct role of partition factors in plasmid replication and incompatibility.

Genetic analysis of a pleiotropic mutant of Klebsiella pneumoniae affected in nitrogen metabolism

Genetic and reversion analyses of a thermosensitive pleiotropic mutant strain of Klebsiella pneumoniae with defects in nitrogen fixation and nitrogen metabolism have shown that the pleiotropie behaviour of mutant is due to a single mutation in a gene designated nim. This gene is contransducible with trp at a frequency of about 30% (using bacteriophage P1) and with cys at a frequency of about 14%. The gene order is cys, trp, nim. The defect in the nim mutant is complemented by the E. coli F' element, F'148.