Structural assignments to the Mycoplasma genitalium proteins show extensive gene duplications and domain rearrangements

The parasitic bacterium Mycoplasma genitalium has a small, reduced genome with close to a basic set of genes. As a first step toward determining the families of protein domains that form the products of these genes, we have used the multiple sequence programs psi-blast and geanfammer to match the sequences of the 467 gene products of M. genitalium to the sequences of the domains that form proteins of known structure [Protein Data Bank (PDB) sequences]. PDB sequences (274) match all of 106 M. genitalium sequences and some parts of another 85; thus, 41% of its total sequences are matched in all or part. The evolutionary relationships of the PDB domains that match M. genitalium are described in the structural classification of proteins (SCOP) database. Using this information, we show that the domains in the matched M. genitalium sequences come from 114 superfamilies and that 58% of them have arisen by gene duplication. This level of duplication is more than twice that found by using pairwise sequence comparisons. The PDB domain matches also describe the domain structure of the matched sequences: just over a quarter contain one domain and the rest have combinations of two or more domains.

Root of the universal tree of life based on ancient aminoacyl-tRNA synthetase gene duplications.

Universal trees based on sequences of single gene homologs cannot be rooted. Iwabe et al. [Iwabe, N., Kuma, K.-I., Hasegawa, M., Osawa, S. & Miyata, T. (1989) Proc. Natl. Acad. Sci. USA 86, 9355-9359] circumvented this problem by using ancient gene duplications that predated the last common ancestor of all living things. Their separate, reciprocally rooted gene trees for elongation factors and ATPase subunits showed Bacteria (eubacteria) as branching first from the universal tree with Archaea (archaebacteria) and Eucarya (eukaryotes) as sister groups. Given its topical importance to evolutionary biology and concerns about the appropriateness of the ATPase data set, an evaluation of the universal tree root using other ancient gene duplications is essential. In this study, we derive a rooting for the universal tree using aminoacyl-tRNA synthetase genes, an extensive multigene family whose divergence likely preceded that of prokaryotes and eukaryotes. An approximately 1600-bp conserved region was sequenced from the isoleucyl-tRNA synthetases of several species representing deep evolutionary branches of eukaryotes (Nosema locustae), Bacteria (Aquifex pyrophilus and Thermotoga maritima) and Archaea (Pyrococcus furiosus and Sulfolobus acidocaldarius). In addition, a new valyl-tRNA synthetase was characterized from the protist Trichomonas vaginalis. Different phylogenetic methods were used to generate trees of isoleucyl-tRNA synthetases rooted by valyl- and leucyl-tRNA synthetases. All isoleucyl-tRNA synthetase trees showed Archaea and Eucarya as sister groups, providing strong confirmation for the universal tree rooting reported by Iwabe et al. As well, there was strong support for the monophyly (sensu Hennig) of Archaea. The valyl-tRNA synthetase gene from Tr. vaginalis clustered with other eukaryotic ValRS genes, which may have been transferred from the mitochondrial genome to the nuclear genome, suggesting that this amitochondrial trichomonad once harbored an endosymbiotic bacterium.

Extent of genomic rearrangement after genome duplication in yeast

Whole-genome duplication approximately 108 years ago was proposed as an explanation for the many duplicated chromosomal regions in Saccharomyces cerevisiae. Here we have used computer simulations and analytic methods to estimate some parameters describing the evolution of the yeast genome after this duplication event. Computer simulation of a model in which 8% of the original genes were retained in duplicate after genome duplication, and 70–100 reciprocal translocations occurred between chromosomes, produced arrangements of duplicated chromosomal regions very similar to the map of real duplications in yeast. An analytical method produced an independent estimate of 84 map disruptions. These results imply that many smaller duplicated chromosomal regions exist in the yeast genome in addition to the 55 originally reported. We also examined the possibility of determining the original order of chromosomal blocks in the ancestral unduplicated genome, but this cannot be done without information from one or more additional species. If the genome sequence of one other species (such as Kluyveromyces lactis) were known it should be possible to identify 150–200 paired regions covering the whole yeast genome and to reconstruct approximately two-thirds of the original order of blocks of genes in yeast. Rates of interchromosome translocation in yeast and mammals appear similar despite their very different rates of homologous recombination per kilobase.

Loss of heterozygosity induced by a chromosomal double-strand break

The repair of chromosomal double-strand breaks (DSBs) is necessary for genomic integrity in all organisms. Genetic consequences of misrepair include chromosomal loss, deletion, and duplication resulting in loss of heterozygosity (LOH), a common finding in human solid tumors. Although work with radiation-sensitive cell lines suggests that mammalian cells primarily rejoin DSBs by nonhomologous mechanisms, alternative mechanisms that are implicated in chromosomal LOH, such as allelic recombination, may also occur. We have examined chromosomal DSB repair between homologs in a gene targeted mammalian cell line at the retinoblastoma (Rb) locus. We have found that allelic recombinational repair occurs in mammalian cells and is increased at least two orders of magnitude by the induction of a chromosomal DSB. One consequence of allelic recombination is LOH at the Rb locus. Some of the repair events also resulted in other types of genetic instability, including deletions and duplications. We speculate that mammalian cells may have developed efficient nonhomologous DSB repair processes to bypass allelic recombination and the potential for reduction to homozygosity.

Nuclear punctate distribution of ALL-1 is conferred by distinct elements at the N terminus of the protein

The ALL-1 gene positioned at 11q23 is directly involved in human acute leukemia either through a variety of chromosome translocations or by partial tandem duplications. ALL-1 is the human homologue of Drosophila trithorax which plays a critical role in maintaining proper spatial and temporal expression of the Antennapedia-bithorax homeotic genes determining the fruit fly’s body pattern. Utilizing specific antibodies, we found that the ALL-1 protein distributes in cultured cells in a nuclear punctate pattern. Several chimeric ALL-1 proteins encoded by products of the chromosome translocations and expressed in transfected cells showed similar speckles. Dissection of the ALL-1 protein identified within its ≈1,100 N-terminal residues three polypeptides directing nuclear localization and at least two main domains conferring distribution in dots. The latter spanned two short sequences conserved with TRITHORAX. Enforced nuclear expression of other domains of ALL-1, such as the PHD (zinc) fingers and the SET motif, resulted in uniform nonpunctate patterns. This indicates that positioning of the ALL-1 protein in subnuclear structures is mediated via interactions of ALL-1 N-terminal elements. We suggest that the speckles represent protein complexes which contain multiple copies of the ALL-1 protein and are positioned at ALL-1 target sites on the chromatin. Therefore, the role of the N-terminal portion of ALL-1 is to direct the protein to its target genes.

Pinus banksiana has at least seven expressed alcohol dehydrogenase genes in two linked groups

The alcohol dehydrogenase (Adh) gene family is much more complex in Pinus banksiana than in angiosperms, with at least seven expressed genes organized as two tightly linked clusters. Intron number and position are highly conserved between P. banksiana and angiosperms. Unlike angiosperm Adh genes, numerous duplications, as large as 217 bp, were observed within the noncoding regions of P. banksiana Adh genes and may be a common feature of conifer genes. A high frequency of duplication over a wide range of scales may contribute to the large genome size of conifers.

Recombination-dependent deletion formation in mammalian cells deficient in the nucleotide excision repair gene ERCC1

Nucleotide excision repair proteins have been implicated in genetic recombination by experiments in Saccharomyces cerevisiae and Drosophila melanogaster, but their role, if any, in mammalian cells is undefined. To investigate the role of the nucleotide excision repair gene ERCC1, the hamster homologue to the S. cerevisiae RAD10 gene, we disabled the gene by targeted knockout. Partial tandem duplications of the adenine phosphoribosyltransferase (APRT) gene then were constructed at the endogenous APRT locus in ERCC1− and ERCC1+ cells. To detect the full spectrum of gene-altering events, we used a loss-of-function assay in which the parental APRT+ tandem duplication could give rise to APRT− cells by homologous recombination, gene rearrangement, or point mutation. Measurement of rates and analysis of individual APRT− products indicated that gene rearrangements (principally deletions) were increased at least 50-fold, whereas homologous recombination was affected little. The formation of deletions is not caused by a general effect of the ERCC1 deficiency on gene stability, because ERCC1− cell lines with a single wild-type copy of the APRT gene yielded no increase in deletions. Thus, deletion formation is dependent on the tandem duplication, and presumably the process of homologous recombination. Recombination-dependent deletion formation in ERCC1− cells is supported by a significant decrease in a particular class of crossover products that are thought to arise by repair of a heteroduplex intermediate in recombination. We suggest that the ERCC1 gene product in mammalian cells is involved in the processing of heteroduplex intermediates in recombination and that the misprocessed intermediates in ERCC1− cells are repaired by illegitimate recombination.

Duplicated genes evolve independently after polyploid formation in cotton

Of the many processes that generate gene duplications, polyploidy is unique in that entire genomes are duplicated. This process has been important in the evolution of many eukaryotic groups, and it occurs with high frequency in plants. Recent evidence suggests that polyploidization may be accompanied by rapid genomic changes, but the evolutionary fate of discrete loci recently doubled by polyploidy (homoeologues) has not been studied. Here we use locus-specific isolation techniques with comparative mapping to characterize the evolution of homoeologous loci in allopolyploid cotton (Gossypium hirsutum) and in species representing its diploid progenitors. We isolated and sequenced 16 loci from both genomes of the allopolyploid, from both progenitor diploid genomes and appropriate outgroups. Phylogenetic analysis of the resulting 73.5 kb of sequence data demonstrated that for all 16 loci (14.7 kb/genome), the topology expected from organismal history was recovered. In contrast to observations involving repetitive DNAs in cotton, there was no evidence of interaction among duplicated genes in the allopolyploid. Polyploidy was not accompanied by an obvious increase in mutations indicative of pseudogene formation. Additionally, differences in rates of divergence among homoeologues in polyploids and orthologues in diploids were indistinguishable across loci, with significant rate deviation restricted to two putative pseudogenes. Our results indicate that most duplicated genes in allopolyploid cotton evolve independently of each other and at the same rate as those of their diploid progenitors. These indications of genic stasis accompanying polyploidization provide a sharp contrast to recent examples of rapid genomic evolution in allopolyploids.

AP-2-null cells disrupt morphogenesis of the eye, face, and limbs in chimeric mice

The homozygous disruption of the mouse AP-2 gene yields a complex and lethal phenotype that results from defective development of the neural tube, head, and body wall. The severe and pleiotropic developmental abnormalities observed in the knockout mouse suggested that AP-2 may regulate several morphogenic pathways. To uncouple the individual developmental mechanisms that are dependent on AP-2, we have now analyzed chimeric mice composed of both wild-type and AP-2-null cells. The phenotypes obtained from these chimeras indicate that there is an independent requirement for AP-2 in the formation of the neural tube, body wall, and craniofacial skeleton. In addition, these studies reveal that AP-2 exerts a major influence on eye formation, which is a critical new role for AP-2 that was masked previously in the knockout mice. Furthermore, we also have uncovered an unexpected influence of AP-2 on limb pattern formation; this influence is typified by major limb duplications. The range of phenotypes observed in the chimeras displays a significant overlap with those caused by teratogenic levels of retinoic acid, strongly suggesting that AP-2 is an important component of the mechanism of action of this morphogen.

Comparative mapping of the human 22q11 chromosomal region and the orthologous region in mice reveals complex changes in gene organization

The region of human chromosome 22q11 is prone to rearrangements. The resulting chromosomal abnormalities are involved in Velo-cardio-facial and DiGeorge syndromes (VCFS and DGS) (deletions), “cat eye” syndrome (duplications), and certain types of tumors (translocations). As a prelude to the development of mouse models for VCFS/DGS by generating targeted deletions in the mouse genome, we examined the organization of genes from human chromosome 22q11 in the mouse. Using genetic linkage analysis and detailed physical mapping, we show that genes from a relatively small region of human 22q11 are distributed on three mouse chromosomes (MMU6, MMU10, and MMU16). Furthermore, although the region corresponding to about 2.5 megabases of the VCFS/DGS critical region is located on mouse chromosome 16, the relative organization of the region is quite different from that in humans. Our results show that the instability of the 22q11 region is not restricted to humans but may have been present throughout evolution. The results also underscore the importance of detailed comparative mapping of genes in mice and humans as a prerequisite for the development of mouse models of human diseases involving chromosomal rearrangements.

Eight novel families of miniature inverted repeat transposable elements in the African malaria mosquito, Anopheles gambiae

Eight novel families of miniature inverted repeat transposable elements (MITEs) were discovered in the African malaria mosquito, Anopheles gambiae, by using new software designed to rapidly identify MITE-like sequences based on their structural characteristics. Divergent subfamilies have been found in two families. Past mobility was demonstrated by evidence of MITE insertions that resulted in the duplication of specific TA, TAA, or 8-bp targets. Some of these MITEs share the same target duplications and similar terminal sequences with MITEs and other DNA transposons in human and other organisms. MITEs in A. gambiae range from 40 to 1340 copies per genome, much less abundant than MITEs in the yellow fever mosquito, Aedes aegypti. Statistical analyses suggest that most A. gambiae MITEs are in highly AT-rich regions, many of which are closely associated with each other. The analyses of these novel MITEs underscored interesting questions regarding their diversity, origin, evolution, and relationships to the host genomes. The discovery of diverse families of MITEs in A. gambiae has important practical implications in light of current efforts to control malaria by replacing vector mosquitoes with genetically modified refractory mosquitoes. Finally, the systematic approach to rapidly identify novel MITEs should have broad applications for the analysis of the ever-growing sequence databases of a wide range of organisms.

Evolution of vertebrate steroid receptors from an ancestral estrogen receptor by ligand exploitation and serial genome expansions

The evolution of novelty in tightly integrated biological systems, such as hormones and their receptors, seems to challenge the theory of natural selection: it has not been clear how a new function for any one part (such as a ligand) can be selected for unless the other members of the system (e.g., a receptor) are already present. Here I show—based on identification and phylogenetic analysis of steroid receptors in basal vertebrates and reconstruction of the sequences and functional attributes of ancestral proteins—that the first steroid receptor was an estrogen receptor, followed by a progesterone receptor. Genome mapping and phylogenetic analyses indicate that the full complement of mammalian steroid receptors evolved from these ancient receptors by two large-scale genome expansions, one before the advent of jawed vertebrates and one after. Specific regulation of physiological processes by androgens and corticoids are relatively recent innovations that emerged after these duplications. These findings support a model of ligand exploitation in which the terminal ligand in a biosynthetic pathway is the first for which a receptor evolves; selection for this hormone also selects for the synthesis of intermediates despite the absence of receptors, and duplicated receptors then evolve affinity for these substances. In this way, novel hormone-receptor pairs are created, and an integrated system of increasing complexity elaborated. This model suggests that ligands for some “orphan” receptors may be found among intermediates in the synthesis of ligands for phylogenetically related receptors.

The elusive function of metallothioneins

Biochemistry and genetics are both required to elucidate the function of macromolecules. There is no question that metallothioneins (MTs) have unique biochemical properties, but genetic experiments have not substantiated the importance of MTs under physiological conditions. Even after thousands of studies describing the structure, biochemical characteristics, tissue distribution, induction, and consequences of genetic disruption and deliberate overexpression, the evolutionary forces that led to the initial appearance, gene duplications, and nearly ubiquitous expression of MTs remain enigmatic.

Maize as a model for the evolution of plant nuclear genomes

The maize genome is replete with chromosomal duplications and repetitive DNA. The duplications resulted from an ancient polyploid event that occurred over 11 million years ago. Based on DNA sequence data, the polyploid event occurred after the divergence between sorghum and maize, and hence the polyploid event explains some of the difference in DNA content between these two species. Genomic rearrangement and diploidization followed the polyploid event. Most of the repetitive DNA in the maize genome is retrotransposable elements, and they comprise 50% of the genome. Retrotransposon multiplication has been relatively recent—within the last 5–6 million years—suggesting that the proliferation of retrotransposons has also contributed to differences in DNA content between sorghum and maize. There are still unanswered questions about repetitive DNA, including the distribution of repetitive DNA throughout the genome, the relative impacts of retrotransposons and chromosomal duplication in plant genome evolution, and the hypothesized correlation of duplication events with transposition. Population genetic processes also affect the evolution of genomes. We discuss how centromeric genes should, in theory, contain less genetic diversity than noncentromeric genes. In addition, studies of diversity in the wild relatives of maize indicate that different genes have different histories and also show that domestication and intensive breeding have had heterogeneous effects on genetic diversity across genes.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.