Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

The chicken beta 2-microglobulin gene is located on a non-major histocompatibility complex microchromosome: a small, G+C-rich gene with X and Y boxes in the promoter.

beta 2-Microglobulin is an essential subunit of major histocompatibility complex (Mhc) class I molecules, which present antigenic peptides to T lymphocytes. We sequenced a number of cDNAs and two genomic clones corresponding to chicken beta 2-microglobulin. The chicken beta 2-microglobulin gene has a similar genomic organization but smaller introns and higher G+C content than mammalian beta 2-microglobulin genes. The promoter region is particularly G+C-rich and contains, in addition to interferon regulatory elements, potential S/W, X, and Y boxes that were originally described for mammalian class II but not class I alpha or beta 2-microglobulin genes. There is a single chicken beta 2-microglobulin gene that has little polymorphism in the coding region. Restriction fragment length polymorphisms from Mhc homozygous lines, Mhc congenic lines, and backcross families, as well as in situ hybridization, show that the beta 2-microglobulin gene is located on a microchromosome different from the one that contains the chicken Mhc. We propose that the structural similarities between the beta 2-microglobulin and Mhc genes in the chicken are due to their presence on microchromosomes and suggest that these features and the microchromosomes appeared by deletion of DNA in the lineage leading to the birds.

Genomic DNA fingerprinting by restriction fragment end labeling.

A typing method for bacteria was developed and applied to several species, including Escherichia coli and Actinobacillus actinomycetemcomitans. Total genomic DNA was digested with a restriction endonuclease, and fragments were enabled with [alpha-32P]dATP by using the Klenow fragment of DNA polymerase and separated by electrophoresis in 6% polyacrylamide/8 M urea (sequencing gel). Depending on the restriction endonuclease and the bacterium, the method produced approximately 30-50 well-separated fragments in the size range of 100-400 nucleotides. For A. actinomycetemcomitans, all strains had bands in common. Nevertheless, many polymorphisms could be observed, and the 31 strains tested could be classified into 29 distinct types. Furthermore, serotype-specific fragments could be assigned for the three serotypes investigated. The method described is very sensitive, allowing more distinct types to be distinguished than other commonly used typing methods. When the method was applied to 10 other clinically relevant bacterial species, both species-specific bands and strain-specific bands were found. Isolates from different locations of one patient showed indistinguishable patterns. Computer-assisted analysis of the DNA fingerprints allowed the determination of similarity coefficients. It is concluded that genomic fingerprinting by restriction fragment end labeling (RFEL) is a powerful and generally applicable technique to type bacterial species.

Genomic evolution during a 10,000-generation experiment with bacteria

Molecular methods are used widely to measure genetic diversity within populations and determine relationships among species. However, it is difficult to observe genomic evolution in action because these dynamics are too slow in most organisms. To overcome this limitation, we sampled genomes from populations of Escherichia coli evolving in the laboratory for 10,000 generations. We analyzed the genomes for restriction fragment length polymorphisms (RFLP) using seven insertion sequences (IS) as probes; most polymorphisms detected by this approach reflect rearrangements (including transpositions) rather than point mutations. The evolving genomes became increasingly different from their ancestor over time. Moreover, tremendous diversity accumulated within each population, such that almost every individual had a different genetic fingerprint after 10,000 generations. As has been often suggested, but not previously shown by experiment, the rates of phenotypic and genomic change were discordant, both across replicate populations and over time within a population. Certain pivotal mutations were shared by all descendants in a population, and these are candidates for beneficial mutations, which are rare and difficult to find. More generally, these data show that the genome is highly dynamic even over a time scale that is, from an evolutionary perspective, very brief.

Detecting immigration by using multilocus genotypes

Immigration is an important force shaping the social structure, evolution, and genetics of populations. A statistical method is presented that uses multilocus genotypes to identify individuals who are immigrants, or have recent immigrant ancestry. The method is appropriate for use with allozymes, microsatellites, or restriction fragment length polymorphisms (RFLPs) and assumes linkage equilibrium among loci. Potential applications include studies of dispersal among natural populations of animals and plants, human evolutionary studies, and typing zoo animals of unknown origin (for use in captive breeding programs). The method is illustrated by analyzing RFLP genotypes in samples of humans from Australian, Japanese, New Guinean, and Senegalese populations. The test has power to detect immigrant ancestors, for these data, up to two generations in the past even though the overall differentiation of allele frequencies among populations is low.

Chromosome-specific molecular organization of maize (Zea mays L.) centromeric regions

A set of oat–maize chromosome addition lines with individual maize (Zea mays L.) chromosomes present in plants with a complete oat (Avena sativa L.) chromosome complement provides a unique opportunity to analyze the organization of centromeric regions of each maize chromosome. A DNA sequence, MCS1a, described previously as a maize centromere-associated sequence, was used as a probe to isolate cosmid clones from a genomic library made of DNA purified from a maize chromosome 9 addition line. Analysis of six cosmid clones containing centromeric DNA segments revealed a complex organization. The MCS1a sequence was found to comprise a portion of the long terminal repeats of a retrotransposon-like repeated element, termed CentA. Two of the six cosmid clones contained regions composed of a newly identified family of tandem repeats, termed CentC. Copies of CentA and tandem arrays of CentC are interspersed with other repetitive elements, including the previously identified maize retroelements Huck and Prem2. Fluorescence in situ hybridization revealed that CentC and CentA elements are limited to the centromeric region of each maize chromosome. The retroelements Huck and Prem2 are dispersed along all maize chromosomes, although Huck elements are present in an increased concentration around centromeric regions. Significant variation in the size of the blocks of CentC and in the copy number of CentA elements, as well as restriction fragment length variations were detected within the centromeric region of each maize chromosome studied. The different proportions and arrangements of these elements and likely others provide each centromeric region with a unique overall structure.

Yersinia pestis, the cause of plague, is a recently emerged clone of Yersinia pseudotuberculosis

Plague, one of the most devastating diseases of human history, is caused by Yersinia pestis. In this study, we analyzed the population genetic structure of Y. pestis and the two other pathogenic Yersinia species, Y. pseudotuberculosis and Y. enterocolitica. Fragments of five housekeeping genes and a gene involved in the synthesis of lipopolysaccharide were sequenced from 36 strains representing the global diversity of Y. pestis and from 12–13 strains from each of the other species. No sequence diversity was found in any Y. pestis gene, and these alleles were identical or nearly identical to alleles from Y. pseudotuberculosis. Thus, Y. pestis is a clone that evolved from Y. pseudotuberculosis 1,500–20,000 years ago, shortly before the first known pandemics of human plague. Three biovars (Antiqua, Medievalis, and Orientalis) have been distinguished by microbiologists within the Y. pestis clone. These biovars form distinct branches of a phylogenetic tree based on restriction fragment length polymorphisms of the locations of the IS100 insertion element. These data are consistent with previous inferences that Antiqua caused a plague pandemic in the sixth century, Medievalis caused the Black Death and subsequent epidemics during the second pandemic wave, and Orientalis caused the current plague pandemic.

Major recent and independent changes in levels and patterns of expression have occurred at the b gene, a regulatory locus in maize

The b locus encodes a transcription factor that regulates the expression of genes that produce purple anthocyanin pigment. Different b alleles are expressed in distinct tissues, causing tissue-specific anthocyanin production. Understanding how phenotypic diversity is produced and maintained at the b locus should provide models for how other regulatory genes, including those that influence morphological traits and development, evolve. We have investigated how different levels and patterns of pigmentation have evolved by determining the phenotypic and evolutionary relationships between 18 alleles that represent the diversity of b alleles in Zea mays. Although most of these alleles have few phenotypic differences, five alleles have very distinct tissue-specific patterns of pigmentation. Superimposing the phenotypes on the molecular phylogeny reveals that the alleles with strong and distinctive patterns of expression are closely related to alleles with weak expression, implying that the distinctive patterns have arisen recently. We have identified apparent insertions in three of the five phenotypically distinct alleles, and the fourth has unique upstream restriction fragment length polymorphisms relative to closely related alleles. The insertion in B-Peru has been shown to be responsible for its unique expression and, in the other two alleles, the presence of the insertion correlates with the phenotype. These results suggest that major changes in gene expression are probably the result of large-scale changes in DNA sequence and/or structure most likely mediated by transposable elements.

BodyMap incorporated PCR-based expression profiling data and a gene ranking system

BodyMap is a human and mouse gene expression database that is based on site-directed 3′-expressed sequence tags generated at Osaka University. To date, it contains more than 300 000 tag sequences from 64 human and 39 mouse tissues. For the recent release, the precise anatomical expression patterns for more than half of the human gene entries were generated by introduced amplified fragment length polymorphism (iAFLP), which is a PCR-based high-throughput expression profiling method. The iAFLP data incorporated into BodyMap describe the relative contents of more than 12 000 transcripts across 30 tissue RNAs. In addition, a newly developed gene ranking system helps users obtain lists of genes that have desired expression patterns according to their significance. BodyMap supports complete transfer of unique data sets and provides analysis that is accessible through the WWW at http://bodymap.ims.u-tokyo.ac.jp.

Protein-coding genes are epigenetically regulated in Arabidopsis polyploids

The fate of redundant genes resulting from genome duplication is poorly understood. Previous studies indicated that ribosomal RNA genes from one parental origin are epigenetically silenced during interspecific hybridization or polyploidization. Regulatory mechanisms for protein-coding genes in polyploid genomes are unknown, partly because of difficulty in studying expression patterns of homologous genes. Here we apply amplified fragment length polymorphism (AFLP)–cDNA display to perform a genome-wide screen for orthologous genes silenced in Arabidopsis suecica, an allotetraploid derived from Arabidopsis thaliana and Cardaminopsis arenosa. We identified ten genes that are silenced from either A. thaliana or C. arenosa origin in A. suecica and located in four of the five A. thaliana chromosomes. These genes represent a variety of RNA and predicted proteins including four transcription factors such as TCP3. The silenced genes in the vicinity of TCP3 are hypermethylated and reactivated by blocking DNA methylation, suggesting epigenetic regulation is involved in the expression of orthologous genes in polyploid genomes. Compared with classic genetic mutations, epigenetic regulation may be advantageous for selection and adaptation of polyploid species during evolution and development.

Isolation of a superfamily of candidate disease-resistance genes in soybean based on a conserved nucleotide-binding site.

The tobacco N and Arabidopsis RPS2 genes, among several recently cloned disease-resistance genes, share highly conserved structure, a nucleotide-binding site (NBS). Using degenerate oligonucleotide primers for the NBS region of N and RPS2, we have amplified and cloned the NBS sequences from soybean. Each of these PCR-derived NBS clones detected low-or moderate-copy soybean DNA sequences and belongs to 1 of 11 different classes. Sequence analysis showed that all PCR clones encode three motifs (P-loop, kinase-2, and kinase-3a) of NBS nearly identical to those in N and RPS2. The intervening region between P-loop and kinase-3a of the 11 classes has high (26% average) amino acid sequence similarity to the N gene although not as high (19% average) to RPS2. These 11 classes represent a superfamily of NBS-containing soybean genes that are homologous to N and RPS2. Each class or subfamily was assessed for its positional association with known soybean disease-resistance genes through near-isogenic line assays, followed by linkage analysis in F2 populations using restriction fragment length polymorphisms. Five of the 11 subfamilies have thus far been mapped to the vicinity of known soybean genes for resistance to potyviruses (Rsv1 and Rpv), Phytophthora root rot (Rps1, Rps2, and Rps3), and powdery mildew (rmd). The conserved N- or RPS2-homologous NBS sequences and their positional associations with mapped soybean-resistance genes suggest that a number of the soybean disease-resistance genes may belong to this superfamily. The candidate subfamilies of NBS-containing genes identified by genetic mapping should greatly facilitate the molecular cloning of disease-resistance genes.

GSP-1 genes are linked to the grain hardness locus (Ha) on wheat chromosome 5D.

An important determinant of wheat grain quality is the hardness of the grain. The trait is controlled by a major locus, Ha, on the short arm of chromosome 5D. Purified starch granules from soft-grained wheats have associated with them 15-kDa polypeptides called grain softness proteins (GSPs) or "friabilins." Genes that encode one family of closely related GSP polypeptides - GSP-1 genes - were mapped using chromosome substitution lines to the group 5 chromosomes. An F2 population segregating for hard and soft alleles at the Ha locus on a near-isogenic background was used in a single-seed study of the inheritance of grain softness and of GSP-1 alleles. Grain softness versus grain hardness was inherited in a 3:1 ratio. The presence versus absence of GSPs in single seed starch preparations was coinherited with grain softness versus hardness. This showed that grain softness is primarily determined by seed, and not by maternal, genotype. In addition, no recombination was detected in 44 F2 plants between GSP-1 restriction fragment length polymorphisms and Ha alleles. Differences between hard and soft wheat grains in membrane structure and lipid extractability have been described and, of the three characterized proteins that are part of the mixture of 15-kDa polypeptides called GSPs, at least two, and probably all three, are proteins that bind polar lipids. The data are interpreted to suggest that the Ha locus may encode one or more members of a large family of lipid-binding proteins.

Barley telomeres shorten during differentiation but grow in callus culture.

Eukaryotic chromosomes terminate with long stretches of short, guanine-rich repeats. These repeats are added de novo by a specialized enzyme, telomerase. In humans telomeres shorten during differentiation, presumably due to the absence of telomerase activity in somatic cells. This phenomenon forms the basis for several models of telomere role in cellular senescence. Barley (Hordeum vulgare L.) telomeres consist of thousands of TTTAGGG repeats, closely resembling other higher eukaryotes. In vivo differentiation and aging resulted in reduction of terminal restriction fragment length paralleled by a decrease of telomere repeat number. Dedifferentiation in callus culture resulted in an increase of the terminal restriction fragment length and in the number of telomere repeats. Long-term callus cultures had very long telomeres. Absolute telomere lengths were genotype dependent, but the relative changes due to differentiation, dedifferentiation, and long-term callus culture were consistent among genotypes. A model is presented to describe the potential role of the telomere length in regulation of a cell's mitotic activity and senescence.

A knob-associated tandem repeat in maize capable of forming fold-back DNA segments: Are chromosome knobs megatransposons?

A class of tandemly repeated DNA sequences (TR-1) of 350-bp unit length was isolated from the knob DNA of chromosome 9 of Zea mays L. Comparative fluorescence in situ hybridization revealed that TR-1 elements are also present in cytologically detectable knobs on other maize chromosomes in different proportions relative to the previously described 180-bp repeats. At least one knob on chromosome 4 is composed predominantly of the TR-1 repeat. In addition, several small clusters of the TR-1 and 180-bp repeats have been found in different chromosomes, some not located in obvious knob heterochromatin. Variation in restriction fragment fingerprints and copy number of the TR-1 elements was found among maize lines and among maize chromosomes. TR-1 tandem arrays up to 70 kilobases in length can be interspersed with stretches of 180-bp tandem repeat arrays. DNA sequence analysis and restriction mapping of one particular stretch of tandemly arranged TR-1 units indicate that these elements may be organized in the form of fold-back DNA segments. The TR-1 repeat shares two short segments of homology with the 180-bp repeat. The longest of these segments (31 bp; 64% identity) corresponds to the conserved region among 180-bp repeats. The polymorphism and complex structure of knob DNA suggest that, similar to the fold-back DNA-containing giant transposons in Drosophila, maize knob DNA may have some properties of transposable elements.

Perturbation of the T cell repertoire in rheumatoid arthritis

Aberrations in the T cell repertoire with the emergence of oligoclonal populations have been described in patients with rheumatoid arthritis (RA). However, the extent of the repertoire perturbations as well as the underlying mechanisms are not known. We now have examined the diversity of the peripheral CD4 T cell repertoire by determining the frequencies of arbitrarily selected T cell receptor (TCR) β-chain sequences. Healthy individuals displayed a highly diverse repertoire, with a median frequency of individual TCR β-chain sequences of 1 in 2.4 × 107 CD4 T cells. In RA patients, the median TCR β-chain frequency was increased 10-fold, indicating marked contraction of the repertoire (P < 0.001). The loss in TCR diversity was not limited to CD4 memory T cells but also involved the compartment of naive T cells, suggesting that it reflected an abnormality in T cell repertoire formation and not a consequence of antigen recognition in the synovium. Also, control patients with chronic inflammatory disease such as hepatitis C expressed a diverse repertoire indistinguishable from that of normals. Telomere length studies indicated an increased replicative history of peripheral CD4 T cells in RA patients, suggesting an enhanced turnover within the CD4 compartment. Compared with age-matched controls, terminal restriction fragment sizes were 1.7 kilobases shorter (P < 0.001). These data demonstrate an altered CD4 T cell homeostasis in RA that may contribute to the autoimmune response as well as to the immunodeficiency in these patients.