Direct genetic analysis by matrix-assisted laser desorption/ionization mass spectrometry

An approach to analyzing single-nucleotide polymorphisms (SNPs) found in the human genome has been developed that couples a recently developed invasive cleavage assay for nucleic acids with detection by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The invasive cleavage assay is a signal amplification method that enables the analysis of SNPs by MALDI-TOF MS directly from human genomic DNA without the need for initial target amplification by PCR. The results presented here show the successful genotyping by this approach of twelve SNPs located randomly throughout the human genome. Conventional Sanger sequencing of these SNP positions confirmed the accuracy of the MALDI-TOF MS analysis results. The ability to unambiguously detect both homozygous and heterozygous genotypes is clearly demonstrated. The elimination of the need for target amplification by PCR, combined with the inherently rapid and accurate nature of detection by MALDI-TOF MS, gives this approach unique and significant advantages in the high-throughput genotyping of large numbers of SNPs, useful for locating, identifying, and characterizing the function of specific genes.

Toward functional genomics in bacteria: Analysis of gene expression in Escherichia coli from a bacterial artificial chromosome library of Bacillus cereus

As the study of microbes moves into the era of functional genomics, there is an increasing need for molecular tools for analysis of a wide diversity of microorganisms. Currently, biological study of many prokaryotes of agricultural, medical, and fundamental scientific interest is limited by the lack of adequate genetic tools. We report the application of the bacterial artificial chromosome (BAC) vector to prokaryotic biology as a powerful approach to address this need. We constructed a BAC library in Escherichia coli from genomic DNA of the Gram-positive bacterium Bacillus cereus. This library provides 5.75-fold coverage of the B. cereus genome, with an average insert size of 98 kb. To determine the extent of heterologous expression of B. cereus genes in the library, we screened it for expression of several B. cereus activities in the E. coli host. Clones expressing 6 of 10 activities tested were identified in the library, namely, ampicillin resistance, zwittermicin A resistance, esculin hydrolysis, hemolysis, orange pigment production, and lecithinase activity. We analyzed selected BAC clones genetically to identify rapidly specific B. cereus loci. These results suggest that BAC libraries will provide a powerful approach for studying gene expression from diverse prokaryotes.

V(D)J recombination is not activated by demethylation of the kappa locus

V(D)J recombination is thought to be regulated by changes in the accessibility of target sites, such as modulation of methylation. To test whether demethylation of the kappa locus can activate recombination, we generated two recombinationally active B cell lines in which the gene for maintenance of genomic DNA methylation, Dnmt1, was flanked with loxP sites. Transduction with a retrovirus expressing both the cre recombinase and green fluorescent protein allowed us to purify recombinationally active cells devoid of methylation. Loss of methylation of the kappa locus was not sufficient to activate recombination, although transcription was activated in one line. It appears that demethylation of the kappa locus is not the rate-limiting step for altering accessibility and thus regulated demethylation does not generate specificity of recombination.

An ethylene-induced cDNA encoding a lipase expressed at the onset of senescence

A cDNA clone encoding a lipase (lipolytic acyl hydrolase) expressed at the onset of petal senescence has been isolated by screening a cDNA expression library prepared from carnation flowers (Dianthus caryophyllus). The cDNA contains the lipase consensus sequence, ITFAGHSLGA, and encodes a 447-amino acid polypeptide with a calculated molecular mass of 50.2 kDa that appears to be a cytosolic protein. Over-expression of the clone in Escherichia coli yielded a protein of the expected molecular weight that proved capable of deesterifying fatty acids from p-nitrophenylpalmitate, tri-linolein, soybean phospholipid, and Tween in both in vitro and in situ assays of enzyme activity. The abundance of the lipase mRNA increases just as carnation flowers begin to senesce, and expression of the gene is also induced by treatment with ethylene. Southern blot analyses of carnation genomic DNA have indicated that the lipase is a single copy gene. The lipase gene is also expressed in carnation leaves and is up-regulated when the leaves are treated with ethylene. Deesterification of membrane lipids and ensuing loss of membrane structural integrity are well established early events of plant senescence, and the expression pattern of this lipase gene together with the lipolytic activity of its cognate protein indicate that it plays a fundamentally central role in mediating the onset of senescence.

Virus-encoded suppressor of posttranscriptional gene silencing targets a maintenance step in the silencing pathway

Certain plant viruses encode suppressors of posttranscriptional gene silencing (PTGS), an adaptive antiviral defense response that limits virus replication and spread. The tobacco etch potyvirus protein, helper component-proteinase (HC-Pro), suppresses PTGS of silenced transgenes. The effect of HC-Pro on different steps of the silencing pathway was analyzed by using both transient Agrobacterium tumefaciens-based delivery and transgenic systems. HC-Pro inactivated PTGS in plants containing a preexisting silenced β-glucuronidase (GUS) transgene. PTGS in this system was associated with both small RNA molecules (21–26 nt) corresponding to the 3′ proximal region of the transcribed GUS sequence and cytosine methylation of specific sites near the 3′ end of the GUS transgene. Introduction of HC-Pro into these plants resulted in loss of PTGS, loss of small RNAs, and partial loss of methylation. These results suggest that HC-Pro targets a PTGS maintenance (as opposed to an initiation or signaling) component at a point that affects accumulation of small RNAs and methylation of genomic DNA.

The MITE family Heartbreaker (Hbr): Molecular markers in maize

Transposable elements are ubiquitous in plant genomes, where they frequently comprise the majority of genomic DNA. The maize genome, which is believed to be structurally representative of large plant genomes, contains single genes or small gene islands interspersed with much longer blocks of retrotransposons. Given this organization, it would be desirable to identify molecular markers preferentially located in genic regions. In this report, the features of a newly described family of miniature inverted repeat transposable elements (MITEs) (called Heartbreaker), including high copy number and polymorphism, stability, and preference for genic regions, have been exploited in the development of a class of molecular markers for maize. To this end, a modification of the AFLP procedure called transposon display was used to generate and display hundreds of genomic fragments anchored in Hbr elements. An average of 52 markers were amplified for each primer combination tested. In all, 213 polymorphic fragments were reliably scored and mapped in 100 recombinant inbred lines derived from a cross between the maize inbreds B73 × Mo17. In this mapping population, Hbr markers are distributed evenly across the 10 maize chromosomes. This procedure should be of general use in the development of markers for other MITE families in maize and in other plant and animal species where MITEs have been identified.

Methylation of the minimal promoter of an embryonic globin gene silences transcription in primary erythroid cells

Methylation of cytosines in the dinucleotide CpG has been shown to suppress transcription of a number of tissue-specific genes, yet the precise mechanism is not fully understood. The vertebrate globin genes were among the first examples in which an inverse correlation was shown between CpG methylation and transcription. We studied the methylation pattern of the 235-bp ρ-globin gene promoter in genomic DNA from primary chicken erythroid cells using the sodium bisulfite conversion technique and found all CpGs in the promoter to be methylated in erythroid cells from adult chickens in which the ρ-globin gene is silent but unmethylated in 5-day (primitive) embryonic red cells in which the gene is transcribed. To elucidate further the mechanism of methylation-induced silencing, an expression construct consisting of 235 bp of 5′ promoter sequence of the ρ-globin gene along with a strong 5′ erythroid enhancer driving a chloramphenicol acetyltransferase reporter gene, ρ-CAT, was transfected into primary avian erythroid cells derived from 5-day embryos. Methylation of just the 235-bp ρ-globin gene promoter fragment at every CpG resulted in a 20- to 30-fold inhibition of transcription, and this effect was not overridden by the presence of potent erythroid-specific enhancers. The ability of the 235-bp ρ-globin gene promoter to bind to a DNA Methyl Cytosine binding Protein Complex (MeCPC) was tested in electrophoretic mobility shift assays utilizing primary avian erythroid cell nuclear extract. The results were that fully methylated but not unmethylated 235-bp ρ-globin gene promoter fragment could compete efficiently for MeCPC binding. These results are a direct demonstration that site-specific methylation of a globin gene promoter at the exact CpGs that are methylated in vivo can silence transcription in homologous primary erythroid cells. Further, these data implicate binding of MeCPC to the promoter in the mechanism of silencing.

Bromination of deoxycytidine by eosinophil peroxidase: A mechanism for mutagenesis by oxidative damage of nucleotide precursors

Oxidants generated by eosinophils during chronic inflammation may lead to mutagenesis in adjacent epithelial cells. Eosinophil peroxidase, a heme enzyme released by eosinophils, generates hypobromous acid that damages tissue in inflammatory conditions. We show that human eosinophils use eosinophil peroxidase to produce 5-bromodeoxycytidine. Flow cytometric, immunohistochemical, and mass spectrometric analyses all demonstrated that 5-bromodeoxycytidine generated by eosinophil peroxidase was taken up by cultured cells and incorporated into genomic DNA as 5-bromodeoxyuridine. Although previous studies have focused on oxidation of chromosomal DNA, our observations suggest another mechanism for oxidative damage of DNA. In this scenario, peroxidase-catalyzed halogenation of nucleotide precursors yields products that subsequently can be incorporated into DNA. Because the thymine analog 5-BrUra mispairs with guanine in DNA, generation of brominated pyrimidines by eosinophils might constitute a mechanism for cytotoxicity and mutagenesis at sites of inflammation.

Diversity Arrays: a solid state technology for sequence information independent genotyping

Here we present the successful application of the microarray technology platform to the analysis of DNA polymorphisms. Using the rice genome as a model, we demonstrate the potential of a high-throughput genome analysis method called Diversity Array Technology, DArT‘. In the format presented here the technology is assaying for the presence (or amount) of a specific DNA fragment in a representation derived from the total genomic DNA of an organism or a population of organisms. Two different approaches are presented: the first involves contrasting two representations on a single array while the second involves contrasting a representation with a reference DNA fragment common to all elements of the array. The Diversity Panels created using this method allow genetic fingerprinting of any organism or group of organisms belonging to the gene pool from which the panel was developed. Diversity Arrays enable rapid and economical application of a highly parallel, solid-state genotyping technology to any genome or complex genomic mixtures.

GeneSplicer: a new computational method for splice site prediction

GeneSplicer is a new, flexible system for detecting splice sites in the genomic DNA of various eukaryotes. The system has been tested successfully using DNA from two reference organisms: the model plant Arabidopsis thaliana and human. It was compared to six programs representing the leading splice site detectors for each of these species: NetPlantGene, NetGene2, HSPL, NNSplice, GENIO and SpliceView. In each case GeneSplicer performed comparably to the best alternative, in terms of both accuracy and computational efficiency.

The Dictyostelium discoideum family of Rho-related proteins

Taking advantage of the ongoing Dictyostelium genome sequencing project, we have assembled >73 kb of genomic DNA in 15 contigs harbouring 15 genes and one pseudogene of Rho-related proteins. Comparison with EST sequences revealed that every gene is interrupted by at least one and up to four introns. For racC extensive alternative splicing was identified. Northern blot analysis showed that mRNAs for racA, racE, racG, racH and racI were present at all stages of development, whereas racJ and racL were expressed only at late stages. Amino acid sequences have been analysed in the context of Rho-related proteins of other organisms. Rac1a/1b/1c, RacF1/F2 and to a lesser extent RacB and the GTPase domain of RacA can be grouped in the Rac subfamily. None of the additional Dictyostelium Rho-related proteins belongs to any of the well-defined subfamilies, like Rac, Cdc42 or Rho. RacD and RacA are unique in that they lack the prenylation motif characteristic of Rho proteins. RacD possesses a 50 residue C-terminal extension and RacA a 400 residue C-terminal extension that contains a proline-rich region, two BTB domains and a novel C-terminal domain. We have also identified homologues for RacA in Drosophila and mammals, thus defining a new subfamily of Rho proteins, RhoBTB.

PCR candidate region mismatch scanning: adaptation to quantitative, high-throughput genotyping

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

Comparative analysis of the gene-dense ACHE/TFR2 region on human chromosome 7q22 with the orthologous region on mouse chromosome 5

Chromosome 7q22 has been the focus of many cytogenetic and molecular studies aimed at delineating regions commonly deleted in myeloid leukemias and myelodysplastic syndromes. We have compared a gene-dense, GC-rich sub-region of 7q22 with the orthologous region on mouse chromosome 5. A physical map of 640 kb of genomic DNA from mouse chromosome 5 was derived from a series of overlapping bacterial artificial chromosomes. A 296 kb segment from the physical map, spanning Ache to Tfr2, was compared with 267 kb of human sequence. We identified a conserved linkage of 12 genes including an open reading frame flanked by Ache and Asr2, a novel cation-chloride cotransporter interacting protein Cip1, Ephb4, Zan and Perq1. While some of these genes have been previously described, in each case we present new data derived from our comparative sequence analysis. Adjacent unfinished sequence data from the mouse contains an orthologous block of 10 additional genes including three novel cDNA sequences that we subsequently mapped to human 7q22. Methods for displaying comparative genomic information, including unfinished sequence data, are becoming increasingly important. We supplement our printed comparative analysis with a new, Web-based program called Laj (local alignments with java). Laj provides interactive access to archived pairwise sequence alignments via the WWW. It displays synchronized views of a dot-plot, a percent identity plot, a nucleotide-level local alignment and a variety of relevant annotations. Our mouse–human comparison can be viewed at http://web.uvic.ca/~bioweb/laj.html. Laj is available at http://bio.cse.psu.edu/, along with online documentation and additional examples of annotated genomic regions.

Defining the minimal length of sequence homology required for selective gene isolation by TAR cloning

The transformation-associated recombination (TAR) cloning technique allows selective and accurate isolation of chromosomal regions and genes from complex genomes. The technique is based on in vivo recombination between genomic DNA and a linearized vector containing homologous sequences, or hooks, to the gene of interest. The recombination occurs during transformation of yeast spheroplasts that results in the generation of a yeast artificial chromosome (YAC) containing the gene of interest. To further enhance and refine the TAR cloning technology, we determined the minimal size of a specific hook required for gene isolation utilizing the Tg.AC mouse transgene as a targeted region. For this purpose a set of vectors containing a B1 repeat hook and a Tg.AC-specific hook of variable sizes (from 20 to 800 bp) was constructed and checked for efficiency of transgene isolation by a radial TAR cloning. When vectors with a specific hook that was ≥60 bp were utilized, ∼2% of transformants contained circular YACs with the Tg.AC transgene sequences. Efficiency of cloning dramatically decreased when the TAR vector contained a hook of 40 bp or less. Thus, the minimal length of a unique sequence required for gene isolation by TAR is ∼60 bp. No transgene-positive YAC clones were detected when an ARS element was incorporated into a vector, demonstrating that the absence of a yeast origin of replication in a vector is a prerequisite for efficient gene isolation by TAR cloning.

Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.