The miniaturized nuclear genome of eukaryotic endosymbiont contains genes that overlap, genes that are cotranscribed, and the smallest known spliceosomal introns.

Chlorarachniophyte algae contain a complex, multi-membraned chloroplast derived from the endosymbiosis of a eukaryotic alga. The vestigial nucleus of the endosymbiont, called the nucleomorph, contains only three small linear chromosomes with a haploid genome size of 380 kb and is the smallest known eukaryotic genome. Nucleotide sequence data from a subtelomeric fragment of chromosome III were analyzed as a preliminary investigation of the coding capacity of this vestigial genome. Several housekeeping genes including U6 small nuclear RNA (snRNA), ribosomal proteins S4 and S13, a core protein of the spliceosome [small nuclear ribonucleoprotein (snRNP) E], and a cip-like protease (clpP) were identified. Expression of these genes was confirmed by combinations of Northern blot analysis, in situ hybridization, immunocytochemistry, and cDNA analysis. The protein-encoding genes are typically eukaryotic in overall structure and their messenger RNAs are polyadenylylated. A novel feature is the abundance of 18-, 19-, or 20-nucleotide introns; the smallest spliceosomal introns known. Two of the genes, U6 and S13, overlap while another two genes, snRNP E and clpP, are cotranscribed in a single mRNA. The overall gene organization is extraordinarily compact, making the nucleomorph a unique model for eukaryotic genomics.

Fluorescence in situ hybridization analysis of keratinocyte growth factor gene amplification and dispersion in evolution of great apes and humans

Keratinocyte growth factor (KGF) is a member of the fibroblast growth factor family. Portions of the gene encoding KGF were amplified during primate evolution and are present in multiple nonprocessed copies in the human genome. Nucleotide analysis of a representative sampling of these KGF-like sequences indicated that they were at least 95% identical to corresponding regions of the KGF gene. To localize these sequences to specific chromosomal sites in human and higher primates, we used fluorescence in situ hybridization. In human, using a cosmid probe encoding KGF exon 1, we assigned the location of the KGF gene to chromosome 15q15–21.1. In addition, copies of KGF-like sequences hybridizing only with a cosmid probe encoding exons 2 and 3 were localized to dispersed sites on chromosome 2q21, 9p11, 9q12–13, 18p11, 18q11, 21q11, and 21q21.1. The distribution of KGF-like sequences suggests a role for alphoid DNA in their amplification and dispersion. In chimpanzee, KGF-like sequences were observed at five chromosomal sites, which were each homologous to sites in human, while in gorilla, a subset of four of these homologous sites was identified; in orangutan two sites were identified, while gibbon exhibited only a single site. The chromosomal localization of KGF sequences in human and great ape genomes indicates that amplification and dispersion occurred in multiple discrete steps, with initial KGF gene duplication and dispersion taking place in gibbon and involving loci corresponding to human chromosomes 15 and 21. These findings support the concept of a closer evolutionary relationship of human and chimpanzee and a possible selective pressure for such dispersion during the evolution of higher primates.

Specific cleavage of chromosomal and plasmid DNA strands in gram-positive and gram-negative bacteria can be detected with nucleotide resolution

A sensitive and precise in vitro technique for detecting DNA strand discontinuities produced in vivo has been developed. The procedure, a form of runoff DNA synthesis on molecules released from lysed bacterial cells, mapped precisely the position of cleavage of the plasmid pMV158 leading strand origin in Streptococcus pneumoniae and the site of strand scission, nic, at the transfer origins of F and the F-like plasmid R1 in Escherichia coli. When high frequency of recombination strains of E. coli were examined, DNA strand discontinuities at the nic positions of the chromosomally integrated fertility factors were also observed. Detection of DNA strand scission at the nic position of F DNA in the high frequency of recombination strains, as well as in the episomal factors, was dependent on sexual expression from the transmissable element, but was independent of mating. These results imply that not only the transfer origins of extrachromosomal F and F-like fertility factors, but also the origins of stably integrated copies of these plasmids, are subject to an equilibrium of cleavage and ligation in vivo in the absence of DNA transfer.

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.

Calibrating bacterial evolution

Attempts to calibrate bacterial evolution have relied on the assumption that rates of molecular sequence divergence in bacteria are similar to those of higher eukaryotes, or to those of the few bacterial taxa for which ancestors can be reliably dated from ecological or geological evidence. Despite similarities in the substitution rates estimated for some lineages, comparisons of the relative rates of evolution at different classes of nucleotide sites indicate no basis for their universal application to all bacteria. However, there is evidence that bacteria have a constant genome-wide mutation rate on an evolutionary time scale but that this rate differs dramatically from the rate estimated by experimental methods.

Identification and classification of 16 new kinesin superfamily (KIF) proteins in mouse genome

KIF (kinesin superfamily) proteins are microtubule-dependent molecular motors that play important roles in intracellular transport and cell division. The extent to which KIFs are involved in various transporting phenomena, as well as their regulation mechanism, are unknown. The identification of 16 new KIFs in this report doubles the existing number of KIFs known in the mouse. Conserved nucleotide sequences in the motor domain were amplified by PCR using cDNAs of mouse nervous tissue, kidney, and small intestine as templates. The new KIFs were studied with respect to their expression patterns in different tissues, chromosomal location, and molecular evolution. Our results suggest that (i) there is no apparent tendency among related subclasses of KIFs of cosegregation in chromosomal mapping, and (ii) according to their tissue distribution patterns, KIFs can be divided into two classes–i.e., ubiquitous and specific tissue-dominant. Further characterization of KIFs may elucidate unknown fundamental phenomena underlying intracellular transport. Finally, we propose a straightforward nomenclature system for the members of the mouse kinesin superfamily.

The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.

Horizontal gene transfer and mutation: Ngrol genes in the genome of Nicotiana glauca

Ngrol genes (NgrolB, NgrolC, NgORF13, and NgORF14) that are similar in sequence to genes in the left transferred DNA (TL-DNA) of Agrobacterium rhizogenes have been found in the genome of untransformed plants of Nicotiana glauca. It has been suggested that a bacterial infection resulted in transformation of Ngrol genes early in the evolution of the genus Nicotiana. Although the corresponding four rol genes in TL-DNA provoked hairy-root syndrome in plants, present-day N. glauca and plants transformed with Ngrol genes did not exhibit this phenotype. Sequenced complementation analysis revealed that the NgrolB gene did not induce adventitious roots because it contained two point mutations. Single-base site-directed mutagenesis at these two positions restored the capacity for root induction to the NgrolB gene. When the NgrolB, with these two base substitutions, was positioned under the control of the cauliflower mosaic virus 35S promoter (P35S), transgenic tobacco plants exhibited morphological abnormalities that were not observed in P35s-RirolB plants. In contrast, the activity of the NgrolC gene may have been conserved after an ancient infection by bacteria. Discussed is the effect of the horizontal gene transfer of the Ngrol genes and mutations in the NgrolB gene on the phenotype of ancient plants during the evolution of N. glauca.

Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

Human cells compromised for p53 function exhibit defective global and transcription-coupled nucleotide excision repair, whereas cells compromised for pRb function are defective only in global repair

After exposure to DNA-damaging agents, the p53 tumor suppressor protects against neoplastic transformation by inducing growth arrest and apoptosis. A series of investigations has also demonstrated that, in UV-exposed cells, p53 regulates the removal of DNA photoproducts from the genome overall (global nucleotide excision repair), but does not participate in an overlapping pathway that removes damage specifically from the transcribed strand of active genes (transcription-coupled nucleotide excision repair). Here, the highly sensitive ligation-mediated PCR was employed to quantify, at nucleotide resolution, the repair of UVB-induced cyclobutane pyrimidine dimers (CPDs) in genetically p53-deficient Li–Fraumeni skin fibroblasts, as well as in human lung fibroblasts expressing the human papillomavirus (HPV) E6 oncoprotein that functionally inactivates p53. Lung fibroblasts expressing the HPV E7 gene product, which similarly inactivates the retinoblastoma tumor-suppressor protein (pRb), were also investigated. pRb acts downstream of p53 to mediate G1 arrest, but has no demonstrated role in DNA repair. Relative to normal cells, HPV E6-expressing lung fibroblasts and Li–Fraumeni skin fibroblasts each manifested defective CPD repair along both the transcribed and nontranscribed strands of the p53 and/or c-jun loci. HPV E7-expressing lung fibroblasts also exhibited reduced CPD removal, but only along the nontranscribed strand. Our results provide striking evidence that transcription-coupled repair, in addition to global repair, are p53-dependent in UV-exposed human fibroblasts. Moreover, the observed DNA-repair defect in HPV E7-expressing cells reveals a function for this oncoprotein in HPV-mediated carcinogenesis, and may suggest a role for pRb in global nucleotide excision repair.

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

RNA editing in Arabidopsis mitochondria effects 441 C to U changes in ORFs

On the basis of the sequence of the mitochondrial genome in the flowering plant Arabidopsis thaliana, RNA editing events were systematically investigated in the respective RNA population. A total of 456 C to U, but no U to C, conversions were identified exclusively in mRNAs, 441 in ORFs, 8 in introns, and 7 in leader and trailer sequences. No RNA editing was seen in any of the rRNAs or in several tRNAs investigated for potential mismatch corrections. RNA editing affects individual coding regions with frequencies varying between 0 and 18.9% of the codons. The predominance of RNA editing events in the first two codon positions is not related to translational decoding, because it is not correlated with codon usage. As a general effect, RNA editing increases the hydrophobicity of the coded mitochondrial proteins. Concerning the selection of RNA editing sites, little significant nucleotide preference is observed in their vicinity in comparison to unedited C residues. This sequence bias is, per se, not sufficient to specify individual C nucleotides in the total RNA population in Arabidopsis mitochondria.

Interactive cloning with the SH3 domain of N-src identifies a new brain specific ion channel protein, with homology to Eag and cyclic nucleotide-gated channels

We have isolated a novel cDNA, that appears to represent a new class of ion channels, by using the yeast two-hybrid system and the SH3 domain of the neural form of Src (N-src) as a bait. The encoded polypeptide, BCNG-1, is distantly related to cyclic nucleotide-gated channels and the voltage-gated channels, Eag and H-erg. BCNG-1 is expressed exclusively in the brain, as a glycosylated protein of ≈132 kDa. Immunohistochemical analysis indicates that BCNG-1 is preferentially expressed in specific subsets of neurons in the neocortex, hippocampus, and cerebellum, in particular pyramidal neurons and basket cells. Within individual neurons, the BCNG-1 protein is localized to either the dendrites or the axon terminals depending on the cell type. Southern blot analysis shows that several other BCNG-related sequences are present in the mouse genome, indicating the emergence of an entire subfamily of ion channel coding genes. These findings suggest the existence of a new type of ion channel, which is potentially able to modulate membrane excitability in the brain and could respond to regulation by cyclic nucleotides.

Engineering the largest RNA virus genome as an infectious bacterial artificial chromosome

The construction of cDNA clones encoding large-size RNA molecules of biological interest, like coronavirus genomes, which are among the largest mature RNA molecules known to biology, has been hampered by the instability of those cDNAs in bacteria. Herein, we show that the application of two strategies, cloning of the cDNAs into a bacterial artificial chromosome and nuclear expression of RNAs that are typically produced within the cytoplasm, is useful for the engineering of large RNA molecules. A cDNA encoding an infectious coronavirus RNA genome has been cloned as a bacterial artificial chromosome. The rescued coronavirus conserved all of the genetic markers introduced throughout the sequence and showed a standard mRNA pattern and the antigenic characteristics expected for the synthetic virus. The cDNA was transcribed within the nucleus, and the RNA translocated to the cytoplasm. Interestingly, the recovered virus had essentially the same sequence as the original one, and no splicing was observed. The cDNA was derived from an attenuated isolate that replicates exclusively in the respiratory tract of swine. During the engineering of the infectious cDNA, the spike gene of the virus was replaced by the spike gene of an enteric isolate. The synthetic virus replicated abundantly in the enteric tract and was fully virulent, demonstrating that the tropism and virulence of the recovered coronavirus can be modified. This demonstration opens up the possibility of employing this infectious cDNA as a vector for vaccine development in human, porcine, canine, and feline species susceptible to group 1 coronaviruses.

Nucleotide sequence of the Kaposi sarcoma-associated herpesvirus (HHV8)

The genome of the Kaposi sarcoma-associated herpesvirus (KSHV or HHV8) was mapped with cosmid and phage genomic libraries from the BC-1 cell line. Its nucleotide sequence was determined except for a 3-kb region at the right end of the genome that was refractory to cloning. The BC-1 KSHV genome consists of a 140.5-kb-long unique coding region flanked by multiple G+C-rich 801-bp terminal repeat sequences. A genomic duplication that apparently arose in the parental tumor is present in this cell culture-derived strain. At least 81 ORFs, including 66 with homology to herpesvirus saimiri ORFs, and 5 internal repeat regions are present in the long unique region. The virus encodes homologs to complement-binding proteins, three cytokines (two macrophage inflammatory proteins and interleukin 6), dihydrofolate reductase, bcl-2, interferon regulatory factors, interleukin 8 receptor, neural cell adhesion molecule-like adhesin, and a D-type cyclin, as well as viral structural and metabolic proteins. Terminal repeat analysis of virus DNA from a KS lesion suggests a monoclonal expansion of KSHV in the KS tumor.