Ancient repeated sequences in the pea and mung bean genomes and implications for genome evolution

Innovations in functional genomics and molecular breeding of pea: exploring advances and opportunities

The garden pea (Pisum sativum L.) is a significant cool-season legume, serving as crucial food sources, animal feed, and industrial raw materials. The advancement of functional genomics over the past two decades has provided substantial theoretical foundations and progress to pea breeding. Notably, the release of the pea reference genome has enhanced our understanding of plant architecture, symbiotic nitrogen fixation (SNF), flowering time, floral organ development, seed development, and stress resistance. However, a considerable gap remains between pea functional genomics and molecular breeding. This review summarizes the current advancements in pea functional genomics and breeding while highlighting the future challenges in pea molecular breeding.

A hypervariable intron of the STAYGREEN locus provides excellent discrimination among Pisum fulvum accessions and reveals evidence for a relatively recent hybridization event with Pisum sativum

An analysis of 82 non-synonymous Pisum fulvum accessions for sequence variation in a fragment of the STAYGREEN (SGR) locus revealed 57 alleles, most of which differed in indel structure. Eight additional P. fulvum accessions, each supposedly synonymous with a different accession of the initial group, were also analyzed. In every case the paired synonymous accessions possessed the same SGR sequence but varied slightly for a 6-trait morphological phenotype, indicating that SGR sequence is a much more reliable indicator of accession identity than is a morphological characterization. SGR sequence analysis confirmed our previous finding that P. fulvum accessions separate into two allele groups. This division was not supported by results of previous studies that were based on sequences distributed across the entire genome, suggesting that the division may have been produced by selection at a nearby locus and that the SGR phylogeny may not be good indicator of overall relationships within the species. One P. fulvum accession, PI 595941 (=JI1796), displayed an SGR sequence outside the variation typical of the species. Instead, its allele resembled alleles limited to a set of Pisum sativum landraces from the Middle East, suggesting hybridization between ancestors of PI 595941 and some primitive form of domesticated P. sativum. With one exception from the extreme northwest corner of Israel, P. fulvum accessions collected north of latitude 35.5° N were fixed for alleles from group A. These northern accessions also displayed greatly reduced SGR sequence diversity compared to group A accessions collected from other regions, suggesting that the northern populations may represent recent extensions of the range of the species. Group B accessions were distributed from Lake Tiberias south and were generally sympatric with the southern group A accessions. Although group B accessions occupied a smaller area than group A, the SGR sequence diversity in this group (28 alleles in 33 accessions) exceeded that for group A.

Genomic region associated with pod color variation in pea (Pisum sativum)

Pea (Pisum sativum) was chosen as the research material by Gregor Mendel to discover the laws of inheritance. Out of seven traits studied by Mendel, genes controlling three traits including pod shape, pod color, and flower position have not been identified to date. With the aim of identifying the genomic region controlling pod color, we determined the genome sequence of a pea line with yellow pods. Genome sequence reads obtained using a Nanopore sequencing technology were assembled into 117,981 contigs (3.3 Gb), with an N50 value of 51.2 kb. A total of 531,242 potential protein-coding genes were predicted, of which 519,349 (2.8 Gb) were located within repetitive sequences (2.8 Gb). The assembled sequences were ordered using a reference as a guide to build pseudomolecules. Subsequent genetic and association analyses led to the identification of a genomic region that controls pea pod color. DNA sequences at this genomic location and transcriptome profiles of green and yellow pod lines were analyzed, and genes encoding 3' exoribonucleases were selected as potential candidates controlling pod color. The results presented in this study are expected to accelerate pan-genome studies in pea and facilitate the identification of the gene controlling one of the traits studied by Mendel.

A reference genome for pea provides insight into legume genome evolution

We report the first annotated chromosome-level reference genome assembly for pea, Gregor Mendel's original genetic model. Phylogenetics and paleogenomics show genomic rearrangements across legumes and suggest a major role for repetitive elements in pea genome evolution. Compared to other sequenced Leguminosae genomes, the pea genome shows intense gene dynamics, most likely associated with genome size expansion when the Fabeae diverged from its sister tribes. During Pisum evolution, translocation and transposition differentially occurred across lineages. This reference sequence will accelerate our understanding of the molecular basis of agronomically important traits and support crop improvement.

Genomic repeat abundances contain phylogenetic signal

A large proportion of genomic information, particularly repetitive elements, is usually ignored when researchers are using next-generation sequencing. Here we demonstrate the usefulness of this repetitive fraction in phylogenetic analyses, utilizing comparative graph-based clustering of next-generation sequence reads, which results in abundance estimates of different classes of genomic repeats. Phylogenetic trees are then inferred based on the genome-wide abundance of different repeat types treated as continuously varying characters; such repeats are scattered across chromosomes and in angiosperms can constitute a majority of nuclear genomic DNA. In six diverse examples, five angiosperms and one insect, this method provides generally well-supported relationships at interspecific and intergeneric levels that agree with results from more standard phylogenetic analyses of commonly used markers. We propose that this methodology may prove especially useful in groups where there is little genetic differentiation in standard phylogenetic markers. At the same time as providing data for phylogenetic inference, this method additionally yields a wealth of data for comparative studies of genome evolution.

Transcriptome sequencing for high throughput SNP development and genetic mapping in Pea

BACKGROUND

Pea has a complex genome of 4.3 Gb for which only limited genomic resources are available to date. Although SNP markers are now highly valuable for research and modern breeding, only a few are described and used in pea for genetic diversity and linkage analysis.

RESULTS

We developed a large resource by cDNA sequencing of 8 genotypes representative of modern breeding material using the Roche 454 technology, combining both long reads (400 bp) and high coverage (3.8 million reads, reaching a total of 1,369 megabases). Sequencing data were assembled and generated a 68 K unigene set, from which 41 K were annotated from their best blast hit against the model species Medicago truncatula. Annotated contigs showed an even distribution along M. truncatula pseudochromosomes, suggesting a good representation of the pea genome. 10 K pea contigs were found to be polymorphic among the genetic material surveyed, corresponding to 35 K SNPs.We validated a subset of 1538 SNPs through the GoldenGate assay, proving their ability to structure a diversity panel of breeding germplasm. Among them, 1340 were genetically mapped and used to build a new consensus map comprising a total of 2070 markers. Based on blast analysis, we could establish 1252 bridges between our pea consensus map and the pseudochromosomes of M. truncatula, which provides new insight on synteny between the two species.

CONCLUSIONS

Our approach created significant new resources in pea, i.e. the most comprehensive genetic map to date tightly linked to the model species M. truncatula and a large SNP resource for both academic research and breeding.

Next-generation sequencing reveals the impact of repetitive DNA across phylogenetically closely related genomes of Orobanchaceae

We used next-generation sequencing to characterize the genomes of nine species of Orobanchaceae of known phylogenetic relationships, different life forms, and including a polyploid species. The study species are the autotrophic, nonparasitic Lindenbergia philippensis, the hemiparasitic Schwalbea americana, and seven nonphotosynthetic parasitic species of Orobanche (Orobanche crenata, Orobanche cumana, Orobanche gracilis (tetraploid), and Orobanche pancicii) and Phelipanche (Phelipanche lavandulacea, Phelipanche purpurea, and Phelipanche ramosa). Ty3/Gypsy elements comprise 1.93%-28.34% of the nine genomes and Ty1/Copia elements comprise 8.09%-22.83%. When compared with L. philippensis and S. americana, the nonphotosynthetic species contain higher proportions of repetitive DNA sequences, perhaps reflecting relaxed selection on genome size in parasitic organisms. Among the parasitic species, those in the genus Orobanche have smaller genomes but higher proportions of repetitive DNA than those in Phelipanche, mostly due to a diversification of repeats and an accumulation of Ty3/Gypsy elements. Genome downsizing in the tetraploid O. gracilis probably led to sequence loss across most repeat types.

Graph-based clustering and characterization of repetitive sequences in next-generation sequencing data

BACKGROUND

The investigation of plant genome structure and evolution requires comprehensive characterization of repetitive sequences that make up the majority of higher plant nuclear DNA. Since genome-wide characterization of repetitive elements is complicated by their high abundance and diversity, novel approaches based on massively-parallel sequencing are being adapted to facilitate the analysis. It has recently been demonstrated that the low-pass genome sequencing provided by a single 454 sequencing reaction is sufficient to capture information about all major repeat families, thus providing the opportunity for efficient repeat investigation in a wide range of species. However, the development of appropriate data mining tools is required in order to fully utilize this sequencing data for repeat characterization.

RESULTS

We adapted a graph-based approach for similarity-based partitioning of whole genome 454 sequence reads in order to build clusters made of the reads derived from individual repeat families. The information about cluster sizes was utilized for assessing the proportion and composition of repeats in the genomes of two model species, Pisum sativum and Glycine max, differing in genome size and 454 sequencing coverage. Moreover, statistical analysis and visual inspection of the topology of the cluster graphs using a newly developed program tool, SeqGrapheR, were shown to be helpful in distinguishing basic types of repeats and investigating sequence variability within repeat families.

CONCLUSIONS

Repetitive regions of plant genomes can be efficiently characterized by the presented graph-based analysis and the graph representation of repeats can be further used to assess the variability and evolutionary divergence of repeat families, discover and characterize novel elements, and aid in subsequent assembly of their consensus sequences.

Patterns of tandem repetition in plant whole genome assemblies

Tandem repeats often confound large genome assemblies. A survey of tandemly arrayed repetitive sequences was carried out in whole genome sequences of the green alga Chlamydomonas reinhardtii, the moss Physcomitrella patens, the monocots rice and sorghum, and the dicots Arabidopsis thaliana, poplar, grapevine, and papaya, in order to test how these assemblies deal with this fraction of DNA. Our results suggest that plant genome assemblies preferentially include tandem repeats composed of shorter monomeric units (especially dinucleotide and 9-30-bp repeats), while higher repetitive units pose more difficulties to assemble. Nevertheless, notwithstanding that currently available sequencing technologies struggle with higher arrays of repeated DNA, major well-known repetitive elements including centromeric and telomeric repeats as well as high copy-number genes, were found to be reasonably well represented. A database including all tandem repeat sequences characterized here was created to benefit future comparative genomic analyses.

Evolution of genome size and complexity in Pinus

BACKGROUND

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA.

CONCLUSIONS/SIGNIFICANCE

Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.

Expressed sequence tags from the Yukon ecotype of Thellungiella reveal that gene expression in response to cold, drought and salinity shows little overlap

Thellungiella salsuginea (also known as T. halophila) is a close relative of Arabidopsis that is very tolerant of drought, freezing, and salinity and may be an appropriate model to identify the molecular mechanisms underlying abiotic stress tolerance in plants. We produced 6578 ESTs, which represented 3628 unique genes (unigenes), from cDNA libraries of cold-, drought-, and salinity-stressed plants from the Yukon ecotype of Thellungiella. Among the unigenes, 94.1% encoded products that were most similar in amino acid sequence to Arabidopsis and 1.5% had no match with a member of the family Brassicaceae. Unigenes from the cold library were more similar to Arabidopsis sequences than either drought- or salinity-induced sequences, indicating that latter responses may be more divergent between Thellungiella and Arabidopsis. Analysis of gene ontology using the best matched Arabidopsis locus showed that the Thellungiella unigenes represented all biological processes and all cellular components, with the highest number of sequences attributed to the chloroplast and mitochondria. Only 140 of the unigenes were found in all three abiotic stress cDNA libraries. Of these common unigenes, 70% have no known function, which demonstrates that Thellungiella can be a rich resource of genetic information about environmental responses. Some of the ESTs in this collection have low sequence similarity with those in Genbank suggesting that they may encode functions that may contribute to Thellungiella's high degree of stress tolerance when compared with Arabidopsis. Moreover, Thellungiella is a closer relative of agriculturally important Brassica spp. than Arabidopsis, which may prove valuable in transferring information to crop improvement programs.

Chromosome rearrangements and transposable elements

There has been limited corroboration to date for McClintock's vision of gene regulation by transposable elements (TEs), although her proposition on the origin of species by TE-induced complex chromosome reorganizations in combination with gene mutations, i.e., the involvement of both factors in relatively sudden formations of species in many plant and animal genera, has been more promising. Moreover, resolution is in sight for several seemingly contradictory phenomena such as the endless reshuffling of chromosome structures and gene sequences versus synteny and the constancy of living fossils (or stasis in general). Recent wide-ranging investigations have confirmed and enlarged the number of earlier cases of TE target site selection (hot spots for TE integration), implying preestablished rather than accidental chromosome rearrangements for nonhomologous recombination of host DNA. The possibility of a partly predetermined generation of biodiversity and new species is discussed. The views of several leading transposon experts on the rather abrupt origin of new species have not been synthesized into the macroevolutionary theory of the punctuated equilibrium school of paleontology inferred from thoroughly consistent features of the fossil record.

Molecular and cytogenetic analysis of repetitive DNA in pea (Pisum sativum L.)

A set of pea DNA sequences representing the most abundant genomic repeats was obtained by combining several approaches. Dispersed repeats were isolated by screening a short-insert genomic library using genomic DNA as a probe. Thirty-two clones ranging from 149 to 2961 bp in size and from 1000 to 39 000/1C in their copy number were sequenced and further characterized. Fourteen clones were identified as retrotransposon-like sequences, based on their homologies to known elements. Fluorescence in situ hybridization using clones of reverse transcriptase and integrase coding sequences as probes revealed that corresponding retroelements were scattered along all pea chromosomes. Two novel families of tandem repeats, named PisTR-A and PisTR-B, were isolated by screening a genomic DNA library with Cot-1 DNA and by employing genomic self-priming PCR, respectively. PisTR-A repeats are 211–212 bp long, their abundance is 2 × 104 copies/1C, and they are partially clustered in a secondary constriction of one chromosome pair with the rest of their copies dispersed on all chromosomes. PisTR-B sequences are of similar abundance (104 copies/1C) but differ from the "A" family in their monomer length (50 bp), high A/T content, and chromosomal localization in a limited number of discrete bands. These bands are located mainly in (sub)telomeric and pericentromeric regions, and their patterns, together with chromosome morphology, allow discrimination of all chromosome types within the pea karyotype. Whereas both tandem repeat families are mostly specific to the genus Pisum, many of the dispersed repeats were detected in other legume species, mainly those in the genus Vicia.Key words: repetitive DNA, plant genome, retroelements, satellite DNA, Pisum sativum.

DNA content for Asian pines parallels New World relatives

This hypothesis is based on an observed correlation between DNA content and thermal regime for fish, zooplankton, salamanders, and some flowering plants. The Asian pine results provided no support for the latitudinal hypothesis; Asian tropical pine species did not have smaller genomes than their temperate or boreal relatives. DNA content of haploid megagametophyte tissue varied from 21.85 pg/C for hard pine Pinus densiflora Sieb. & Zucc. to 29.59 pg/C for soft pine Pinus bungeana Zucc. Pinus merkusii Jung. & De Vriese (29.63 pg/C) was the exceptional hard pine, with a genome size larger than many soft pines. The mean DNA content of Asian soft pines exceeded Asian hard pines (Δ 3.22 pg/C), a parallel to the previously reported trends for New World pines. No continental effect was detected. Based on 46 pines species sampled in centers of species diversity in Asian and the New World, soft pines had mean DNA content which exceeded hard pines by 4.97 pg/C.Key words: gymnosperms, conifers, laser flow cytometry, megagametophytes, C values, phylogeny.

The gene distribution in the genomes of pea, tomato and date palm

The vast majority of genes of maize, rice, barley and wheat are contained in long gene-rich regions (collectively called the 'gene space') separated by long gene-empty regions. The gene space covers a narrow, 0.8-1.6%, GC range, possibly because of the presence of abundant transposons. Here we report that the gene space is not an exclusive property of Gramineae, because it also exists in the large genome of pea (5000 Mb). Moreover, the gene space is not just dependent upon genome size, since a gene space is found in rice (415 Mb), but not in Arabidopsis (120 Mb), nor in two other plants investigated in the present work, date palm (250 Mb) and tomato (1000 Mb).

Characterization of the tomato (Lycopersicon esculentum) genome using in vitro and in situ DNA reassociation

A detailed in vitro study of the kinetics of DNA renaturation, i.e., a C0t analysis, can be used to determine the size of a genome, the relative proportions of single-copy and repetitive sequences, and the complexity of genome components. Despite the dual importance of tomato (Lycopersicon esculentum) as a model for basic plant research and as a crop plant, to the best of our knowledge a C0t analysis has never been published for this species. This is probably due to difficulties associated with isolating sufficient quantities of polyphenol-free nuclear DNA from tomato. Recently we developed a technique for isolating milligram quantities of purified DNA from tomato nuclei, and we used DNA isolated in this manner to prepare a C0t curve for the tomato genome. Analysis of the C0t data indicates that the tomato genome (1C) consists of approximately 0.86 pg of DNA. In agreement with earlier molecular studies, the C0t analysis suggests that most (~73%) of the tomato genome is composed of single-copy sequences. Since 77% of the DNA in tomato chromosomes is found in constitutive heterochromatin, many of the single-copy sequences must reside in heterochromatin, an unexpected arrangement, considering that the constitutive heterochromatin of most species is predominantly repetitive DNA. To determine the distribution of repetitive and single-copy DNA along tomato pachytene chromosomes, we used hydroxyapatite-purified C0t fractions as probes for fluorescence in situ hybridization (FISH). Our FISH results indicate that highly repetitive DNA hybridizes almost exclusively with heterochromatin. While single-copy DNA comprises most of the DNA in euchromatin, heterochromatin contains the majority of single-copy DNA sequences, an observation consistent with our C0t data and previous cytological studies.Key words: tomato, Lycopersicon esculentum, genome size, heterochromatin, euchromatin, DNA reassociation, fluorescence in situ hybridization, FISH, C0t.

Molecular evolution: codes, clocks, genes and genomes

The discoveries, advancements and continuing controversies in the field of molecular evolution are reviewed. Topics summarized are (1) the evolution of the genetic code, (2) gene evolution including the demonstration of homology, estimation of sequence divergence, phylogenetic trees, the molecular clock and the origin of genes and gene families by various genetic mechanisms, and (3) eukaryotic genome evolution, including the highly repeated satellite sequences, the interspersed and potentially mobile repeated sequences and the unique sequence fraction of the genome.

Two classes of differentially regulated glutamine synthetase genes are expressed in the soybean nodule: a nodule-specific class and a constitutively expressed class

We have characterized two sets of cDNA clones representing the glutamine synthetase (GS) mRNA in soybean nodules. Using the 3'-untranslated regions of a representative member of each set, as gene member(s) specific probes, we have shown that one set of the GS genes are expressed in a nodule-specific manner, while the other set is expressed in other tissues, besides the nodules. The nodule-specific GS genes are expressed in a developmentally regulated manner in the nodules, independent of the onset of nitrogen fixation. The other class of GS genes is expressed constitutively in all tissues tested, but its expression level is dramatically enhanced in nodules following onset of N2 fixation. The latter set of genes is also expressed in cotyledons of germinating seedlings in a developmentally regulated manner. Analysis of hybrid select translation products and genomic Southern blots suggests that multiple gene members in each class are expressed in the nodules.

Modulation of glutamine synthetase gene expression in tobacco by the introduction of an alfalfa glutamine synthetase gene in sense and antisense orientation: molecular and biochemical analysis

A glutamine synthetase (GS) cDNA isolated from an alfalfa cell culture cDNA library was found to represent a cytoplasmic GS. The full-length alfalfa GS1 coding sequence, in both sense and antisense orientation and under the transcriptional control of the cauliflower mosaic virus 35S promoter, was introduced into tobacco. Leaves of tobacco plants transformed with the sense construct contained greatly elevated levels of GS transcript and GS polypeptide which assembled into active enzyme. Leaves of the plants transformed with the antisense GS1 construct showed a significant decrease in the level of both GS1 and GS2 polypeptides and GS activity, but did not show any significant decrease in the level of endogenous GS mRNA. We have proposed that antisense inhibition using a heterologous antisense GS RNA occurs at the level of translation. Our results also suggest that the post-translational assembly of GS subunits into a holoenzyme requires an additional factor(s) and is under regulatory control.

Synthetic cryIIIA gene from Bacillus thuringiensis improved for high expression in plants

A 1974 bp synthetic gene was constructed from chemically synthesized oligonucleotides in order to improve transgenic protein expression of the cryIIIA gene from Bacillus thuringiensis var. tenebrionis in transgenic tobacco. The crystal toxin genes (cry) from B. thuringiensis are difficult to express in plants even when under the control of efficient plant regulatory sequences. We identified and eliminated five classes of sequence found throughout the cryIIIA gene that mimic eukaryotic processing signals and which may be responsible for the low levels of transcription and translation. Furthermore, the GC content of the gene was raised from 36% to 49% and the codon usage was changed to be more plant-like. When the synthetic gene was placed behind the cauliflower mosaic virus 35S promoter and the alfalfa mosaic virus translational enhancer, up to 0.6% of the total protein in transgenic tobacco plants was cryIIIA as measured from immunoblot analysis. Bioassay data using potato beetle larvae confirmed this estimate.

Isolation and characterization of cDNA clones for RNA species induced by substituted benzenesulfonamides in corn

A search of compounds capable of inducing specific gene expression in plants without affecting growth and development led to the examination of changes in the pattern of gene expression in corn after treatment with substituted benzenesulfonamide herbicide safeners. Following hydroponic treatment of corn with the safener N-(aminocarbonyl)-2-chlorobenzenesulfonamide (2-CBSU), the specific induction of new translatable mRNA species was observed. Replicate copies of a cDNA library made using RNA from 2-CBSU-treated corn roots were differentially screened with cDNA probes made from either the same mRNA fraction used for library construction or mRNA isolated from roots treated with 2-chlorobenzenesulfonamide (2-CBSA), an inactive analog of the safener. Colonies showing hybridization only with the probe made using mRNA from 2-CBSU-treated roots were further characterized to assess the specificity of the induction and decay of the corresponding induced RNA species. RNA blot analyses showed two clones, designated In2-1 and In2-2, contained plasmids that hybridized to RNAs that were induced from an undetectable background in corn roots within 30 minutes after treatment with 2-CBSU. Leaf and meristem tissues showed similar inductions of the In2-1 and In2-2 RNA species after a delay of several hours. In addition, both RNA species were induced in corn by foliar application of 2-CBSU. In contrast, neither RNA species was induced following stress treatments of plants. These results indicate a substituted benzenesulfonamide safener might be used with the promoters from the In2-1 and In2-2 genes to develop a new inducible gene expression system for plants.

Down-regulation of phytochrome mRNA abundance by red light and benzyladenine in etiolated cucumber cotyledons

Northern blot analysis revealed that a single 4.2 kb phytochrome mRNA species was detectable in cotyledons excised from five-day-old etiolated cucumber seedlings. Intact etiolated five-day-old cucumber seedlings were given a red light or benzyladenine treatment, and cotyledons were harvested at various times following treatment. The abundance of phytochrome mRNA in the cotyledons was quantitated using 32P-labeled RNA probes and slot blot analysis. By 2 h after irradiation the phytochrome mRNA level was reduced to 40% of the initial abundance and reaccumulation began by 3 h after irradiation. Reaccumulation of phytochrome mRNA to the time-zero dark control level was achieved by 10 h after treatment. A decrease in phytochrome mRNA abundance was evident by 2 h after benzyladenine treatment, and a maximal reduction to 45% of the time-zero dark control was attained by 4 h after treatment. No recovery of the phytochrome mRNA level was evident by 8 h after benzyladenine treatment. The abundance of actin mRNA was unaffected by benzyladenine treatment.

Characterization of the kafirin gene family from sorghum reveals extensive homology with zein from maize

Electrophoretic analysis of translation products of polyadenylated RNA isolated from mid-maturation sorghum seed in the presence of [(35)S]met, [(3)H]leu, or [(3)H]val revealed two major proteins of kDa and 21 kDa. These products were not detected when [(3)H]lys was supplied as the radioactive substrate. Under similar electrophoretic conditions, kafirin (a major seed storage prolamin of sorghum), migrated as two bands of 22 kDa and 19 kDa. Sequence analysis of two cDNA clones (pSK8 and pSKR2) from sorghum seed mRNA revealed them to be highly homologous with each other and to the 22 kDa zeins from maize, suggesting that they represented kafirin cDNAs. Compared with pSKR2, pSK8 had an insertion of 24 nucleotides and a deletion of 24 nucleotides, so that the coding regions were nearly identical in length. The deduced amino acid sequence for these cDNA clones reveals that kafirin, like zein, is rich in glutamine and nonpolar amino acids, but contains no lysine. Both kafirin and zein have a 21 amino acid signal peptide exhibiting 80% homology and eight copies of a repetitive amino acid block in the C-terminal domain with the consensus: infI (supP) LL finP (supA) LN infQ (supP) LALANPAAYLQQQQ.The kafirin cDNAs were used as probes to screen a sorghum genomic library; one genomic clone (λGK.1) was sequenced and found to be very similar (97.8%) to the pSK8 cDNA clone. Clone λGK.1 contains features typical for a functional gene in that the intronless open reading frame encoding 268 amino acids is flanked at the 5' end by sequences corresponding to the CAAT and TATA promoter boxes (positioned at about -60 and -30 bp, respectively, from the transcriptional initiation site), and at the 3' end by a consensus polyadenylation signal. In common with zein genomic clones, kafirin clones contain a 15 basepair consensus sequence centered at postion -320 relative to the transcriptional initiation site. Under similar hybridization conditions, genomic reconstruction analysis using an oligonucleotide probe indicated the presence of less than 20 copies of kafirin per haploid sorghum genome compared with approximatley 140 copies of zein per haploid maize genome.

Isolation, characterization and sequencing of a novel repetitive DNA from the mung bean Vigna radiata

A family of highly reiterated, small (approx. 300 bp) sequences has been identified in DNA of the mung bean Vigna radiata. The members are extensively interspersed throughout the chromosomes with some clustering. They also occur extrachromosomally. There is no tissue-specificity to the repeat family but it is highly species-specific. The repetitive DNA hybridizes to total RNA as well as to polyadenylated RNA isolated from germinated mung beans. It has analogy with the human AluI family in the mode of isolation, size, genomic distribution, copy number and transcribability though they do not share any sequence homology. A repetitive DNA clone was selected from a shotgun genomic library of mung bean DNA in pBR322. The average copy number of the cloned repeat is estimated to be 8 x 10(4) per haploid genome, and thus constitutes approx. 5% of the total mung bean genome. The genomic organization and transcription of the cloned repeat is reported. Sequencing of the cloned repetitive DNA reveals the presence of the number of direct and inverted repeats and some short palindromic sequences.

Cloning of cDNA for phytochrome from etiolated Cucurbita and coordinate photoregulation of the abundance of two distinct phytochrome transcripts

We have isolated several cDNA clones for phytochrome from a dicot, Cucurbita pepo L. cv. Black Beauty (zucchini), and have used them to study the regulation of Cucurbita phytochrome mRNA levels. A cDNA library was constructed from poly(A)(+) RNA isolated from etiolated Cucurbita hypocotyl hooks and enriched for phytochrome mRNA by size fractionation. This library was screened with a (32)P-labeled fragment isolated from an Avena phytochrome cDNA clone. Several putative phytochrome clones were isolated and mapped by restriction endonuclease analysis. On the basis of this analysis there is no evidence for the expression of multiple phytochrome genes in Cucurbita. Recent sequence analysis has confirmed that the largest of these clones, pFMD1 (∼3.6 kb), does indeed encode phytochrome and that it contains the entire amino acid coding sequence for Cucurbita phytochrome (33). RNA blot analysis has revealed that two polyadenylated phytochrome transcripts (∼5.6 kb and ∼4.2 kb) are present in both cotyledons and hypocotyl hooks of Cucurbita. In etiolated Cucurbita seedlings given a saturating pulse of red light, the abundance of both transcripts coordinately declines to 50-60% of the dark levels within 3 h and reaccumulates to dark levels within 24 h. Reversal of induction of this response by a far-red light pulse immediately following red light treatment is not observed, which is in contrast to the far-red reversibility of the red light promoted decrease in phytochrome mRNA abundance observed in Avena (6). Etiolated seedlings transferred to continuous white light also show a coordinate decrease in the levels of the two RNAs to ∼40% of the dark levels within 3 h. The magnitude of the light-induced decline in phytochrome mRNA abundance in Cucurbita is substantially less than the decrease previously reported for Avena (6).

Light-stimulated accumulation of the peroxisomal enzymes hydroxypyruvate reductase and serine:glyoxylate aminotransferase and their translatable mRNAs in cotyledons of cucumber seedlings

The development of peroxisomal enzymes in cotyledons of cucumber seedlings is strongly dependent on light. In light-grown seedlings, activities of two peroxisomal enzymes, hydroxypyruvate reductase (HPR) and serine: glyoxylate aminotransferase (SGAT), were barely detectable until three days postimbibition, after which time both activities increased rapidly and linearly for at least three days. In the dark, the activities of these enzymes increased slightly over the same time period, but only to about 5% to 10% of 7-day light-induced levels. When 51/2-day dark-grown seedlings were transferred into white light, activities of HPR and SGAT began to increase after approximately 8 h. HPR protein was shown by an immunoprecipitation assay to increase concurrently with enzymatic activity in both light- and dark-grown cotyledons. Immunoblotting results suggested that the amounts of SGAT-A and SGAT-B, the two subunits of SGAT, also developed along with SGAT activity. The relative levels of translatable mRNAs encoding HPR, SGAT-A, and SGAT-B were also light-dependent, and increased with a developmental pattern similar to enzyme activity and protein levels in light- and dark-grown cotyledons. In 51/2-day dark-grown cotyledons that were transferred to the light, translatable mRNAs for SGAT-A and SGAT-B began to increase within 1 h of illumination and continued of increase rapidly and linearly for the next 24 h in the light to a new steady-state level that was 45 times that of dark controls. Translatable HPR mRNA exhibited a biphasic pattern of accumulation, with a three-fold increase during the first 6 h of illumination, followed by an additional six-fold increase between 8 and 24 h. The accumulation of translationally active mRNA for both enzymes preceded the accumulation of the corresponding protein and enzyme activity by about 8 h. Our data suggest that the rise in enzyme activity depends on an increase in translatable mRNA for these enzymes and is regulated at a pretranslational level, most likely involving transcription of new mRNA.

Use of sodium trichloroacetate and mung bean nuclease to increase sensitivity and precision during transcript mapping

An improved method for mapping RNA transcript boundaries by the nuclease protection technique is presented. This method exploits the large (greater than 20 degrees C) difference in the thermal stability of RNA:DNA and DNA:DNA duplexes in concentrated chaotropic salt solutions. At 45 degrees C in 3.0 M sodium trichloroacetate RNA:DNA hybridization is very efficient but DNA:DNA duplexes remain completely denatured. For many applications, this solvent system can eliminate the need to prepare probes that are free of competing or irrelevant DNA molecules. Fifty- to 100-fold more RNA:DNA hybridization is observed when reassociation is performed in 3.0 M sodium trichloroacetate than in solutions containing high concentrations of formamide. A comparison of the use of S1 nuclease or mung bean nuclease suggests that mung bean nuclease can produce more precise and less ambiguous nuclease protection patterns.

Expression of nodule-specific glutamine synthetase genes during nodule development in soybeans

A cDNA clone (pcPvNGS-01) to glutamine synthetase (GS) mRNA from root nodules of Phaseolus vulgaris showed cross-hybridization to GS and mRNA from soybean root nodules, thus allowing its use as a probe to study the expression of GS genes during root nodule development in soybeans. Hybrid-select translation of root and nodule RNA of soybean with DNA from pcPvNGS-01, followed by 2D gel electrophoresis, showed six peptides in the root and an additional four peptides in the nodule which represent nodule-specific glutamine synthetase (GSn) gene products. The GSn gene products appeared for the first time between day 11 and 12 after infection, either concomitant with the onset of nitrogenase activity or immediately following it. The levels of expression of the GSn and leghemoglobin genes were not affected in young Fix(-) nodules formed by Bradyrhizobium japonicum strains that are defective in nitrogenase activity, suggesting that the induction of these two sets of host genes take place independent of nitrogenase activity. However, in Fix(-) nodules that are incapable of maintaining the peribacteroid membrane, GSn gene products were not detected while 1ba, 1bc2 and 1bc3 appeared. In both the timing of appearance during root nodule development and the effect of different bacterial mutations on the expression, GSn genes differ from most other nodulin genes examined (30), suggesting different regulatory mechanisms.

Genome size variation among north american minnows (Cyprinidae). I. Distribution of the variation in five species

Genome sizes (nuclear DNA contents) were examined spectrophotometrically from ten individuals of each of five species of North American cyprinid fishes (minnows). The distributions of DNA values both within and between the five species were essentially continuous and normal. Differences between individuals within populations were significant and contributed to approximately 16 per cent of the total variation. Variation between individuals within species ranged from 4.7-13.5 per cent and averaged ca. 7.4 per cent. Variation between species ranged from 0-9.5 per cent and the average difference between any species pair was ca. 4.6 per cent. Statistical analyses showed that the methodology used was sufficient to detect significant differences in genome size as small as 2-3 per cent. Consideration of these data lead to the following tentative conclusions: (i) changes in genome size in cyprinids appear small in amount, frequent in occurrence, to involve both gains and losses of DNA, and to be cumulative and independent in effect; (ii) differences within and between cyprinid taxa are likely the result of accumulations of small changes in DNA quantity; and (iii) the primary focus of quantitative DNA variation in cyprinids is between individuals within populations. The extent of DNA quantity variation which occurs within species would appear to preclude any direct relationship between genome size variation and many of the organismal parameters (including speciation) which differentiate the five species. In short, the data suggest that a significant fraction of the cyprinid genome, perhaps more than 10 per cent, is free to vary quantitatively without phenotypic constraint or biological consequence. This fraction is considerably larger than that theoretically needed for the structural gene component.

Phytochrome regulation of phytochrome mRNA abundance

Pure phytochrome RNA sequence synthesized in an SP6-derived in vitro transcription system has been used as a standard to quantitate phytochrome mRNA abundance in Avena seedlings using a filter hybridization assay. In 4-day-old etiolated Avena seedlings phytochrome mRNA represents ∼0.1% of the total poly(A)(+) RNA. Irradiation of such seedlings with a saturating red-light pulse or continuous white light induces a decline in this mRNA that is detectable within 30 min and results in a 50% reduction by ∼60 min and >90% reduction within 5 h. The effect of the red-light pulse is reversed, approximately to the level of the far-red control, by an immediately subsequent far-red pulse. In seedlings maintained in extended darkness after the red-light pulse, the initial rapid decline in phytochrome mRNA level is followed by a slower reaccumulation such that 50-60% of the initial abundance is reached by 48 h. White-light grown seedlings transferred to darkness exhibit a similar accumulation of phytochrome mRNA that is accelerated by removal of residual Pfr with a far-red light pulse at the start of the dark period. The data establish that previously reported phytochrome-regulated changes in translatable phytochrome mRNA levels result from changes in the physical abundance of this mRNA rather than from altered translatability.

Use of a recombination-deficient phage lambda system to construct wheat genomic libraries

The poor cloning efficiency of wheat (Triticum aestivum cv. Yamhill) DNA in conventional cloning vectors has previously prevented the preparation of complete genomic libraries. We show here that while wheat DNA does not clone efficiently using the vector Ch4A, it can be cloned efficiently using Ch32. Ch32 clones are red- gam+ and therefore can be propagated on recombination-deficient hosts. These results suggest that instability of wheat sequences in conventional lambda vector systems has frustrated previous attempts to prepare libraries.

Structure of soybean Kunitz trypsin inhibitor mRNA determined from cDNA by using oligodeoxynucleotide primers

Oligodeoxynucleotides complementary to the deduced mRNA sequence of soybean Kunitz trypsin inhibitor (KTI) were used to prime the synthesis of cDNA from soybean cotyledon total poly(A) RNA. The primed cDNA was used to select clones from a Glycine max cotyledon cDNA library. Two out of twelve hybridizing clones were shown to contain KTI cDNA. The nucleotide sequence of one clone, pSTI 9-2, was determined and it was found to encompass the complete protein coding region of KTI excet for three C-terminal residues. Trypsin inhibitor is synthesized with a 25 amino acid hydrophobic N-terminal sequence presumed to be a signal peptide. The mature polypeptide encoded by pSTI 9-2 agrees with the published amino acid composition of KTI, but contains two discrepancies at the peptide sequence level.

Does the concentration of DNA (Co) and the time of incubation (t) as parameters of Cot influence the thermal stability of the DNA duplexes?

It has been shown in a previous paper (8) that the prime product of reassociation of related DNA sequences under open experimental conditions are mismatched duplexes which undergo 'maturation' upon further incubation. Due to this feature, the Tm value of the duplexes of a large number of DNAs is strongly dependent on the Cot value. Here we present data showing that the Tm of the duplexes of such type of DNAs depends also on the concentration of DNA in the range of one and the same Cot value. The significance of this finding in studying the taxonomic relationship by DNA-DNA hybridisation is discussed.

Phytochrome control of RNA levels in developing pea and mung-bean leaves

We have examined phytochrome effects on the abundance of transcripts from several nuclear and chloroplast genes in buds of dark-grown pea seedlings and primary leaves of dark-grown mung-bean seedlings. Probes for nuclear-coded RNAs were selected from a library of cDNA clones and included those corresponding to the small subunit (SS) of ribulosebisphosphate carboxylase and a chlorophyll a/b binding protein (AB). Transcripts from chloroplast genes for RuBP carboxylase large subunit (LS) and a 32,000-dalton photosystem II polypeptide (PII) were assayed with cloned fragments of the chloroplast genome. In addition, we present data on transcripts from a number of other nuclear genes of unknown function, several of which change in abundance during light-induced development. Transcript levels were measured as a proportion of total RNA by a dot blot assay in which RNA from different tissues or stages is fixed to nitrocellulose and hybridized with (32)P-labeled probes prepared from cloned DNAs. Several patterns of induction can be seen. For example, although both SS and AB RNAs show positive, red/far-red reversible responses in both pea and mung bean, in pea buds the induction ratio for SS RNA is much higher than that for AB RNA, while just the reverse is true for mung-bean leaves. In addition, treatment with lowfluence red light produces full induction of the pea AB RNA, while SS RNA in the same tissue does not reach a maximum steady-state level until after about 24 h of supplementary high-intensity white light. In pea buds, chloroplast genes (LS, PII) also show clear responses to phytochrome, as measured by the steady-state levels of their RNA products. Chloroplast DNA levels (as a fraction of the total cellular DNA) show the same response pattern, which may indicate that in peas many of the light effects we see are related to a general stimulation of chloroplast development. In mung beans, the levels of plastid DNA and RNA are already quite high in the leaves of 7-d dark-grown seedlings, and light effects are much less pronounced. The results are consistent with the notion that chloroplast development is arrested at a later stage in dark-grown mung-bean leaves than in etiolated pea buds.

Improved yield of full-length phaseolin cDNA clones by controlling premature anticomplementary DNA synthesis

The manner in which the known enzymatic properties of reverse transcriptase may limit the length of double stranded cDNAs to be used in cloning was studied. Results here suggest that the well-documented ability of reverse transcriptase to synthesize anticomplementary DNA can, if unrecognized, seriously limit the final yield of full-length cDNA clones. Under conditions which permitted anticomplementary DNA synthesis during the synthesis of the first cDNA strand, no full-length cDNA clones for phaseolin, the principal storage proteins ofPhaseolus vulgaris, were detected among 19 phaseolin-positive cDNA clones. When anticomplementary DNA synthesis was inhibited with 4 mM sodium pyrophosphate, 5 full-length cDNA clones were identified among 45 phaseolin-positive clones. In both cases, the products of the first strand synthesis were C-tailed and the second strand synthesized by reverse transcriptase using oligo(dG) as a primer. The implications of anticomplementary synthesis in cloning methods involving the use of S1 nuclease are discussed. In addition, a rapid, one-step procedure for obtaining partial clones which equally represent the 5' and 3' ends of the RNA is presented.

Repeat sequence interspersion in coding DNA of peas does not reflect that in total pea DNA

The pattern of sequence organization in the regions of the pea genome near sequences coding for mRNA differs significantly from that in total DNA. Interspersion of repeated and single copy sequences is so extensive that 85% of 1300 nucleotide-long fragments contain highly repetitive sequences (about 5000 copies per haploid genome). However, data presented here demonstrate that sequences which code for mRNA are enriched in the small fraction of fragments which do not contain these highly repetitive sequences. Thus, in contrast to the great majority of other sequences in the genome, most mRNA coding sequences are not located within 1300 nucleotides of highly repetitive elements. Moreover, our data indicate that those repeats (if any) which are closely associated with mRNA coding sequences belong to low copy number families characterized by an unusually low degree of sequence divergence.

Chloroplast DNA rearrangements are more frequent when a large inverted repeat sequence is lost

We examined the arrangement of sequences common to seven angiosperm chloroplast genomes. The chloroplast DNAs of spinach, petunia and cucumber are essentially colinear. They share with the corn chloroplast genome a large inversion of approximately 50 kb relative to the genomes of three legumes--mung bean, pea and broad bean. There is one additional rearrangement, a second, smaller inversion within the 50 kb inversion, which is specific to the corn genome. These two changes are the only detectable rearrangements that have occurred during the evolution of the species examined (corn, spinach, petunia, cucumber and mung bean) whose chloroplast genomes contain a large inverted repeat sequence of 22-25 kb. In contrast, we find extensive sequence rearrangements in comparing the pea and broad bean genomes, both of which have deleted one entire segment of the inverted repeat, and also in comparing each of these to the mung bean genome. Thus there is a relatively stable arrangement of sequences in those genomes with the inverted repeat and a much more dynamic arrangement in those that have lost it. We discuss several explanations for this correlation, including the possibility that the inverted repeat may play a direct role in maintaining a conserved arrangement of chloroplast DNA sequences.