Comparative genomic analysis of the mouse and rat amylase multigene family

The rat and mouse amylase gene families were characterized using sequence data from the UCSC genome assembly. We found that the rat genome contains one amylase-1 and two amylase-2 genes, lying close to one another on the same chromosome. Detailed analysis revealed at least six additional amylase pseudogenes in the rat genome in the region adjacent to the amylase-2 genes. In contrast, the mouse has one amylase-1 gene and five amylase-2 genes; the latter are tandemly and systematically arranged on the same chromosome and were generated by segmental duplication. Detailed analysis revealed that the mouse has two amylase pseudogenes, located 5' to the five amylase-2 segments. Thus, the amylase genes of mouse and rat tend to be amplified; the sequences of some of them are fixed while others have become pseudogenes during evolution. This is the second report of amylase genomic organization in mammals and the first in the rodents.     

High-resolution mapping of repetitive DNA by in situ hybridization: molecular and chromosomal features of prominent dispersed and discretely localized DNA families from the wild beet species Beta procumbens

Members of three prominent DNA families of Beta procumbens have been isolated as Sau3A repeats. Two families consisting of repeats of about 158 bp and 312 bp are organized as satellite DNAs (Sau3A satellites I and II), whereas the third family with a repeat length of 202 bp is interspersed throughout the genome. Multi-colour fluorescence in situ hybridization was used for physical mapping of the DNA families, and has shown that these tandemly organized families occur in large heterochromatic and DAPI positive blocks. The Sau3A satellite I hybridized exclusively around or near the centromeres of 10, 11 or 12 chromosomes. The Sau3A satellite family I showed high intraspecific variability and high-resolution physical mapping was performed on pachytene chromosomes using differentially labelled repeats. The physical order of satellite subfamily arrays along a chromosome was visualized and provided evidence that large arrays of plant satellite repeats are not contiguous and consist of distinct subfamily domains. Re-hybridization of a heterologous rRNA probe to mitotic metaphase chromosomes revealed that the 18S-5.8S-25S rRNA genes are located at subterminal position on one chromosome pair missing repeat clusters of the Sau3A satellite family I. It is known that arrays of Sau3A satellite I repeats are tightly linked to a nematode (Heterodera schachtii) resistance gene and our results show that the gene might be located close to the centromere. Large arrays of the Sau3A satellite II were found in centromeric regions of 16 chromosomes and, in addition, a considerable interspersion of repeats over all chromosomes was observed. The family of interspersed 202 bp repeats is uniformly distributed over all chromosomes and largely excluded from the rRNA gene cluster but shows local amplification in some regions. Southern hybridization has shown that all three families are specific for genomes of the section Procumbentes of the genus Beta.     

Localization of three DNA segments encompassing tRNA genes to human chromosomes 1, 5, and 16: proposed mechanism and significance of tRNA gene dispersion

The chromosomal locations of three cloned human DNA fragments encompassing tRNA genes have been determined by Southern analysis of human-rodent somatic cell hybrid DNAs with subfragments from these cloned genes and flanking sequences used as hybridization probes. These three DNA segments have been assigned to human chromosomes 1, 5, and 16, and homologous sequences are probably located on chromosome 14 and a separate locus on chromosome 1. These studies, combined with previous results, indicate that tRNA genes and pseudogenes are dispersed on at least seven different human chromosomes and suggest that these sequences will probably be found on most, if not all, human chromosomes. Short (8-12 nucleotide) direct terminal repeats flank many of the dispersed tRNA genes. The presence of these flanking repeats, combined with the dispersion of tRNA genes throughout the human genome, suggests that many of these genes may have arisen by an RNA-mediated retroposition mechanism. The possible functional significance of this gene dispersion is considered.     

Gene decay in archaea

The gene-dense chromosomes of archaea and bacteria were long thought to be devoid of pseudogenes, but with the massive increase in available genome sequences, whole genome comparisons between closely related species have identified mutations that have rendered numerous genes inactive. Comparative analyses of sequenced archaeal genomes revealed numerous pseudogenes, which can constitute up to 8.6% of the annotated coding sequences in some genomes. The largest proportion of pseudogenes is created by gene truncations, followed by frameshift mutations. Within archaeal genomes, large numbers of pseudogenes contain more than one inactivating mutation, suggesting that pseudogenes are deleted from the genome more slowly in archaea than in bacteria. Although archaea seem to retain pseudogenes longer than do bacteria, most archaeal genomes have unique repertoires of pseudogenes.     

The organization of cytoplasmic ribosomal protein genes in the Arabidopsis genome

Eukaryotic ribosomes are made of two components, four ribosomal RNAs, and approximately 80 ribosomal proteins (r-proteins). The exact number of r-proteins and r-protein genes in higher plants is not known. The strong conservation in eukaryotic r-protein primary sequence allowed us to use the well-characterized rat (Rattus norvegicus) r-protein set to identify orthologues on the five haploid chromosomes of Arabidopsis. By use of the numerous expressed sequence tag (EST) accessions and the complete genomic sequence of this species, we identified 249 genes (including some pseudogenes) corresponding to 80 (32 small subunit and 48 large subunit) cytoplasmic r-protein types. None of the r-protein genes are single copy and most are encoded by three or four expressed genes, indicative of the internal duplication of the Arabidopsis genome. The r-proteins are distributed throughout the genome. Inspection of genes in the vicinity of r-protein gene family members confirms extensive duplications of large chromosome fragments and sheds light on the evolutionary history of the Arabidopsis genome. Examination of large duplicated regions indicated that a significant fraction of the r-protein genes have been either lost from one of the duplicated fragments or inserted after the initial duplication event. Only 52 r-protein genes lack a matching EST accession, and 19 of these contain incomplete open reading frames, confirming that most genes are expressed. Assessment of cognate EST numbers suggests that r-protein gene family members are differentially expressed.     

Genomic rearrangements resulting in PLP1 deletion occur by nonhomologous end joining and cause different dysmyelinating phenotypes in males and females

In the majority of patients with Pelizaeus-Merzbacher disease, duplication of the proteolipid protein gene PLP1 is responsible, whereas deletion of PLP1 is infrequent. Genomic mechanisms for these submicroscopic chromosomal rearrangements remain unknown. We identified three families with PLP1 deletions (including one family described elsewhere) that arose by three distinct processes. In one family, PLP1 deletion resulted from a maternal balanced submicroscopic insertional translocation of the entire PLP1 gene to the telomere of chromosome 19. PLP1 on the 19qtel is probably inactive by virtue of a position effect, because a healthy male sibling carries the same der(19) chromosome along with a normal X chromosome. Genomic mapping of the deleted segments revealed that the deletions are smaller than most of the PLP1 duplications and involve only two other genes. We hypothesize that the deletion is infrequent, because only the smaller deletions can avoid causing either infertility or lethality. Analyses of the DNA sequence flanking the deletion breakpoints revealed Alu-Alu recombination in the family with translocation. In the other two families, no homologous sequence flanking the breakpoints was found, but the distal breakpoints were embedded in novel low-copy repeats, suggesting the potential involvement of genome architecture in stimulating these rearrangements. In one family, junction sequences revealed a complex recombination event. Our data suggest that PLP1 deletions are likely caused by nonhomologous end joining.     

Evolution of immunoglobulin kappa chain variable region genes in vertebrates

The major source of immunoglobulin diversity is variation in DNA sequence among multiple copies of variable region (V) genes of the heavy- and light-chain multigene families. In order to clarify the evolutionary pattern of the multigene family of immunoglobulin light kappa chain V region (V kappa) genes, phylogenetic analyses of V kappa genes from humans and other vertebrate species were conducted. The results obtained indicate that the V kappa genes so far sequenced can be grouped into three major monophyletic clusters, the cartilaginous fish, bony fish and amphibian, and mammalian clusters, and that the cartilaginous fish cluster first separated from the rest of the V kappa genes and then the remaining two clusters diverged. The mammalian V kappa genes can further be divided into 10 V kappa groups, 7 of which are present in the human genome. Human and mouse V kappa genes from different V kappa groups are intermingled rather than clustered on the chromosome, and there are a large number of pseudogenes scattered on the chromosome. This indicates that the chromosomal locations of V kappa genes have been shuffled many times by gene duplication, deletion, and transposition in the evolutionary process and that many genes have become nonfunctional during this process. This mode of evolution is consistent with the model of birth-and-death evolution rather than with the model of concerted evolution. An analysis of duplicate V kappa functional genes and pseudogenes in the human genome has indicated that pseudogenes evolve faster than functional genes but that the rate of nonsynonymous nucleotide substitution in the complementarity-determining regions of V kappa genes has been enhanced by positive Darwinian selection.      

Structure, sequence, chromosomal location, and evolution of the human ferredoxin gene family

Ferredoxin is an iron-sulfur protein that serves as an electron transport intermediate for mitochondrial cytochromes P450 involved in steroid, vitamin D, and bile acid metabolism. We cloned and characterized the human ferredoxin gene family, which includes two expressed genes and two pseudogenes. Sequence analysis of this gene family revealed that it encodes only one protein product. The expressed genes were assigned to chromosome 11 and pseudogenes to chromosomes 20 and 21 by identifying single-copy probes from each gene segment and hybridizing them to DNA from rodent-human hybrid cells. The pseudogenes lacked introns and contained numerous mutations, including insertion, deletion, and substitution which rendered them inactive. They were 96% and 85% homologous to the expressed gene, yet they were only 78% homologous with each other. The intronless nature, higher diversity among themselves, and distinct chromosomal location of the pseudogenes suggests that they arose by independent, retroposon-mediated events.     

Chromosomal localization of three human poly(A)-binding protein genes and four related pseudogenes

In humans, the poly(A)-binding proteins (PABPs) comprise a small nuclear isoform and a conserved gene family that displays at least three functional proteins: PABP1, inducible PABP (iPABP), and PABP3, plus four pseudogenes (1, 2, 3, and PABP4). In situ hybridization of PABP3 cDNA as the probe on metaphasic chromosomes have revealed five possible loci for this gene family at 2q21-q22, 13q11-q12, 12q13.3-q15, 8q22, and 3q24-q25. Amplifications of specific DNA fragments from a human-rodent somatic cell hybrid panel have allowed us to associate PABP1 and PABP3 with 8q22 and 13q11-q12, respectively. The iPABP gene has been assigned to chromosome 1. This result, compared with radiation hybrid database information, strengthens the location of this gene to 1p32-p36. The pseudogenes PABP4, 1, and 2 have been assigned to chromosomes 15, 4, and 14, respectively. Three loci detected on chromosome spreads are not associated with any amplified fragment. They might represent other related PABP genes not yet identified.     

Ten families of variant genes encoded in subtelomeric regions of multiple chromosomes of Plasmodium chabaudi,a malaria species that undergoes antigenic variation in the laboratory mouse

The chromosome ends of human malaria parasites harbour many genes encoding proteins that are exported to the surface of infected red cells, often being involved in host-parasite interactions and immune evasion. Unlike other murine malaria parasites Plasmodium chabaudi undergoes antigenic variation during passage in the laboratory mouse and hence is a model suitable for investigation of switching mechanisms. However, little is known about the subtelomeric regions of P. chabaudi chromosomes and its variable antigens. Here we report 80 kb of sequence from an end of one P. chabaudi chromosome. Hybridization of probes spanning this region to two dimensional pulsed field gels of the genome revealed 10 multicopy gene families located exclusively in subtelomeric regions of multiple P. chabaudi chromosomes, interspersed amongst multicopy intergenic regions. Hence all chromosomes share a common subtelomeric structure, presumably playing a similar role in spatial positioning as the P. falciparum Rep20 sequence. Expression in blood stages, domains characteristic of surface antigens and copy numbers between four and several hundred per genome, indicate a functional role in antigenic variation for some of these families. We identify members of the cir family, as well as novel genes, that although clearly homologous to cir have large low complexity regions in the predicted extracellular domains. Although all families have homologues in other rodent Plasmodium species, four were previously not known to be subtelomeric. Six have homologues in human and simian malarias.     

Molecular evolution of receptor-like kinase genes in hexaploid wheat. Independent evolution of orthologs after polyploidization and mechanisms of local rearrangements at paralogous loci

Hexaploid wheat is a young polyploid species and represents a good model to study mechanisms of gene evolution after polyploidization. Recent studies at the scale of the whole genome have suggested rapid genomic changes after polyploidization but so far the rearrangements that have occurred in terms of gene content and organization have not been analyzed at the microlevel in wheat. Here, we have isolated members of a receptor kinase (Lrk) gene family in hexaploid and diploid wheat, Aegilops tauschii, and barley (Hordeum vulgare). Phylogenetic analysis has allowed us to establish evolutionary relationships (orthology versus paralogy) between the different members of this gene family in wheat as well as with Lrk genes from barley. It also demonstrated that the sequences of the homoeologous Lrk genes evolved independently after polyploidization. In addition, we found evidence for gene loss during the evolution of wheat and barley. Analysis of large genomic fragments isolated from nonorthologous Lrk loci showed a high conservation of the gene content and gene organization at these loci on the homoeologous group 1 chromosomes of wheat and barley. Finally, sequence comparison of two paralogous fragments of chromosome 1B showed a large number of local events (sequence duplications, deletions, and insertions), which reveal rearrangements and mechanisms for genome enlargement at the microlevel.     

A bacterial genome in flux: the twelve linear and nine circular extrachromosomal DNAs in an infectious isolate of the Lyme disease spirochete Borrelia burgdorferi

We have determined that Borrelia burgdorferi strain B31 MI carries 21 extrachromosomal DNA elements, the largest number known for any bacterium. Among these are 12 linear and nine circular plasmids, whose sequences total 610 694 bp. We report here the nucleotide sequence of three linear and seven circular plasmids (comprising 290 546 bp) in this infectious isolate. This completes the genome sequencing project for this organism; its genome size is 1 521 419 bp (plus about 2000 bp of undetermined telomeric sequences). Analysis of the sequence implies that there has been extensive and sometimes rather recent DNA rearrangement among a number of the linear plasmids. Many of these events appear to have been mediated by recombinational processes that formed duplications. These many regions of similarity are reflected in the fact that most plasmid genes are members of one of the genome's 161 paralogous gene families; 107 of these gene families, which vary in size from two to 41 members, contain at least one plasmid gene. These rearrangements appear to have contributed to a surprisingly large number of apparently non-functional pseudogenes, a very unusual feature for a prokaryotic genome. The presence of these damaged genes suggests that some of the plasmids may be in a period of rapid evolution. The sequence predicts 535 plasmid genes >/=300 bp in length that may be intact and 167 apparently mutationally damaged and/or unexpressed genes (pseudogenes). The large majority, over 90%, of genes on these plasmids have no convincing similarity to genes outside Borrelia, suggesting that they perform specialized functions.     

Genome of papaya, a fast growing tropical fruit tree

Papaya is a major fruit crop in tropical and subtropical regions worldwide. It has long been recognized as a nutritious and healthy fruit rich in vitamins A and C. Its small genome, unique aspects of nascent sex chromosomes, and agricultural importance are justifications for sequencing the genome. A female plant of the transgenic variety SunUp was selected for sequencing its genome to avoid the complication of assembling the XY chromosomes in a male or hermaphrodite plant. The draft genome revealed fewer genes than sequenced genomes of flowering plants, partly due to its lack of genome wide duplication since the ancient triplication event shared by eudicots. Most gene families have fewer members in papaya, including significantly fewer disease resistance genes. However, striking amplifications in gene number were found in some functional groups, including MADS-box genes, starch synthases, and volatiles that might affect the speciation and adaptation of papaya. The draft genome was used to clone a gene controlling fruit flesh color and to accelerate the construction of physical maps of sex chromosomes in papaya. Finishing the papaya genome and re-sequencing selected genomes in the family will further facilitate papaya improvement and the study of genome and sex chromosome evolution in angiosperms, particularly in Caricaceae.     

The XRCC2 and XRCC3 repair genes are required for chromosome stability in mammalian cells.

The irs1 and irs1SF hamster cell lines are mutated for the XRCC2 and XRCC3 genes, respectively. Both show heightened sensitivity to ionizing radiation and particularly to the DNA cross-linking chemical mitomycin C (MMC). Frequencies of spontaneous chromosomal aberration have previously been reported to be higher in these two cell lines than in parental, wild-type cell lines. Microcell-mediated chromosome transfer was used to introduce complementing or non-complementing human chromosomes into each cell line. irs1 cells received human chromosome 7 (which contains the human XRCC2 gene) or, as a control, human chromosome 4. irs1SF cells received human chromosome 14 (which contains the XRCC3 gene) or human chromosome 7. For each set of hybrid cell lines, clones carrying the complementing human chromosome recovered MMC resistance to near-wild-type levels, while control clones carrying noncomplementing chromosomes remained sensitive to MMC. Fluorescence in situ hybridization with a human-specific probe revealed that the human chromosome in complemented clones remained intact in almost all cells even after extended passage. However, the human chromosome in noncomplemented clones frequently underwent chromosome rearrangements including breaks, deletions, and translocations. Chromosome aberrations accumulated slowly in the noncomplemented clones over subsequent passages, with some particular deletions and unbalanced translocations persistently transmitted throughout individual subclones. Our results indicate that the XRCC2 and XRCC3 genes, which are now considered members of the RAD51 gene family, play essential roles in maintaining chromosome stability during cell division. This may reflect roles in DNA repair, possibly via homologous recombination.     

Comparative Analysis of Evolution in a Rodent Histone H2a Pseudogene

Sequences were obtained from five species of rodents that are orthologous to an H2a histone pseudogene from Mus musculus. The pseudogene is part of the cluster of replication-dependent histone genes found on Mus musculus chromosome 13. Comparative analysis of these five sequences together with the previously published sequence from M. musculus shows that this gene has likely been a pseudogene throughout the evolution of the genus Mus, while the gene from Rattus norvegicus is likely functional. Three large (>20 bp) deletions were found among the Mus pseudogenes, a feature that is very unusual compared to surveys of processed pseudogenes. In addition, there are two single-base deletions and one 4-bp insertion among the Mus pseudogenes. The species distributions of one of the large deletions and the 4-bp insertion require either independent insertions of an identical sequence, independent deletions with identical boundaries, or a deletion followed by precise reintegration of the original sequence. The evidence favors the hypothesis of multiple deletions with identical boundaries. The ``coding' regions of the Mus pseudogenes show a much reduced level of among-species variability in the 3′ half of the pseudogene, compared both to the 5′ half and to flanking sequences. This supports a hypothesis that the 3′ end of the pseudogene is the target of frequent gene conversion by functional H2a genes. Received: 1 April 1997 / Accepted: 12 June 1997