Ancient mtDNA sequences in the human nuclear genome: A potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer’s disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent ρ° (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the ρ° cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five “AD” missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from ρ° cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.

Identification of a Second Myosin-II in Schizosaccharomyces pombe:: Myp2p Is Conditionally Required for Cytokinesis

As in many eukaryotic cells, fission yeast cytokinesis depends on the assembly of an actin ring. We cloned myp2+, a myosin-II in Schizosaccharomyces pombe, conditionally required for cytokinesis. myp2+, the second myosin-II identified in S. pombe, does not completely overlap in function with myo2+. The catalytic domain of Myp2p is highly homologous to known myosin-IIs, and phylogenetic analysis places Myp2p in the myosin-II family. The Myp2p sequence contains well-conserved ATP- and actin-binding motifs, as well as two IQ motifs. However, the tail sequence is unusual, since it is predicted to form two long coiled-coils separated by a stretch of sequence containing 19 prolines. Disruption of myp2+ is not lethal but under nutrient limiting conditions cells lacking myp2+ function are multiseptated, elongated, and branched, indicative of a defect in cytokinesis. The presence of salt enhances these morphological defects. Additionally, Δmyp2 cells are cold sensitive in high salt, failing to form colonies at 17°C. Thus, myp2+ is required under conditions of stress, possibly linking extracellular growth conditions to efficient cytokinesis and cell growth. GFP-Myp2p localizes to a ring in the middle of late mitotic cells, consistent with a role in cytokinesis. Additionally, we constructed double mutants of Δmyp2 with temperature-sensitive mutant strains defective in cytokinesis. We observed synthetic lethal interactions between Δmyp2 and three alleles of cdc11ts, as well as more modest synthetic interactions with cdc14ts and cdc16ts, implicating myp2+ function for efficient cytokinesis under normal conditions.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last ≈1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.

Cryptocyanin, a crustacean molting protein: Evolutionary link with arthropod hemocyanins and insect hexamerins

Cryptocyanin, a copper-free hexameric protein in crab (Cancer magister) hemolymph, has been characterized and the amino acid sequence has been deduced from its cDNA. It is markedly similar in sequence, size, and structure to hemocyanin, the copper-containing oxygen-transport protein found in many arthropods. Cryptocyanin does not bind oxygen, however, and lacks three of the six highly conserved copper-binding histidine residues of hemocyanin. Cryptocyanin has no phenoloxidase activity, although a phenoloxidase is present in the hemolymph. The concentration of cryptocyanin in the hemolymph is closely coordinated with the molt cycle and reaches levels higher than hemocyanin during premolt. Cryptocyanin resembles insect hexamerins in the lack of copper, molt cycle patterns of biosynthesis, and potential contributions to the new exoskeleton. Phylogenetic analysis of sequence similarities between cryptocyanin and other members of the hemocyanin gene family shows that cryptocyanin is closely associated with crustacean hemocyanins and suggests that cryptocyanin arose as a result of a hemocyanin gene duplication. The presence of both hemocyanin and cryptocyanin in one animal provides an example of how insect hexamerins might have evolved from hemocyanin. Our results suggest that multiple members of the hemocyanin gene family—hemocyanin, cryptocyanin, phenoloxidase, and hexamerins—may participate in two vital functions of molting animals, oxygen binding and molting. Cryptocyanin may provide important molecular data to further investigate evolutionary relationships among all molting animals.

Molecular studies suggest that cartilaginous fishes have a terminal position in the piscine tree

The Chondrichthyes (cartilaginous fishes) are commonly accepted as being sister group to the other extant Gnathostomata (jawed vertebrates). To clarify gnathostome relationships and to aid in resolving and dating the major piscine divergences, we have sequenced the complete mtDNA of the starry skate and have included it in phylogenetic analysis along with three squalomorph chondrichthyans—the common dogfish, the spiny dogfish, and the star spotted dogfish—and a number of bony fishes and amniotes. The direction of evolution within the gnathostome tree was established by rooting it with the most closely related non-gnathostome outgroup, the sea lamprey, as well as with some more distantly related taxa. The analyses placed the chondrichthyans in a terminal position in the piscine tree. These findings, which also suggest that the origin of the amniote lineage is older than the age of the oldest extant bony fishes (the lungfishes), challenge the evolutionary direction of several morphological characters that have been used in reconstructing gnathostome relationships. Applying as a calibration point the age of the oldest lungfish fossils, 400 million years, the molecular estimate placed the squalomorph/batomorph divergence at ≈190 million years before present. This dating is consistent with the occurrence of the earliest batomorph (skates and rays) fossils in the paleontological record. The split between gnathostome fishes and the amniote lineage was dated at ≈420 million years before present.

Identification of a third distinct estrogen receptor and reclassification of estrogen receptors in teleosts

This paper describes three distinct estrogen receptor (ER) subtypes: ERα, ERβ, and a unique type, ERγ, cloned from a teleost fish, the Atlantic croaker Micropogonias undulatus; the first identification of a third type of classical ER in vertebrate species. Phylogenetic analysis shows that ERγ arose through gene duplication from ERβ early in the teleost lineage and indicates that ERγ is present in other teleosts, although it has not been recognized as such. The Atlantic croaker ERγ shows amino acid differences in regions important for ligand binding and receptor activation that are conserved in all other ERγs. The three ER subtypes are genetically distinct and have different distribution patterns in Atlantic croaker tissues. In addition, ERβ and ERγ fusion proteins can each bind estradiol-17β with high affinity. The presence of three functional ERs in one species expands the role of ER multiplicity in estrogen signaling systems and provides a unique opportunity to investigate the dynamics and mechanisms of ER evolution.

Evolutionary specialization of the nuclear targeting apparatus

The α- and β-karyopherins (Kaps), also called importins, mediate the nuclear transport of proteins. All α-Kaps contain a central domain composed of eight approximately 40 amino acid, tandemly arranged, armadillo-like (Arm) repeats. The number and order of these repeats have not changed since the common origin of fungi, plants, and mammals. Phylogenetic analysis suggests that the various α-Kaps fall into two groups, α1 and α2. Whereas animals encode both types, the yeast genome encodes only an α1-Kap. The β-Kaps are characterized by 14–15 tandemly arranged HEAT motifs. We show that the Arm repeats of α-Kaps and the HEAT motifs of β-Kaps are similar, suggesting that the α-Kaps and β-Kaps (and for that matter, all Arm and HEAT repeat-containing proteins) are members of the same protein superfamily. Phylogenetic analysis indicates that there are at least three major groups of β-Kaps, consistent with their proposed cargo specificities. We present a model in which an α-independent β-Kap progenitor gave rise to the α-dependent β-Kaps and the α-Kaps.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Did homeodomain proteins duplicate before the origin of angiosperms, fungi, and metazoa?

Homeodomain proteins are transcription factors that play a critical role in early development in eukaryotes. These proteins previously have been classified into numerous subgroups whose phylogenetic relationships are unclear. Our phylogenetic analysis of representative eukaryotic sequences suggests that there are two major groups of homeodomain proteins, each containing sequences from angiosperms, metazoa, and fungi. This result, based on parsimony and neighbor-joining analyses of primary amino acid sequences, was supported by two additional features of the proteins. The two protein groups are distinguished by an insertion/deletion in the homeodomain, between helices I and II. In addition, an amphipathic alpha-helical secondary structure in the region N terminal of the homeodomain is shared by angiosperm and metazoan sequences in one group. These results support the hypothesis that there was at least one duplication of homeobox genes before the origin of angiosperms, fungi, and metazoa. This duplication, in turn, suggests that these proteins had diverse functions early in the evolution of eukaryotes. The shared secondary structure in angiosperm and metazoan sequences points to an ancient conserved functional domain.

Conserved interaction between distinct Krüppel-associated box domains and the transcriptional intermediary factor 1 β

The Krüppel-associated box (KRAB) domain, originally identified as a 75-aa sequence present in numerous Krüppel-type zinc-finger proteins, is a potent DNA-binding-dependent transcriptional repression domain that is believed to function through interaction with the transcriptional intermediary factor 1 (TIF1) β. On the basis of sequence comparison and phylogenetic analysis, we have recently defined three distinct subfamilies of KRAB domains. In the present study, individual members of each subfamily were tested for transcriptional repression and interaction with TIF1β and two other closely related family members (TIF1α and TIF1γ). All KRAB variants were shown, (i) to repress transcription when targeted to DNA through fusion to a heterologous DNA-binding domain in mammalian cells, and (ii) to interact specifically with TIF1β, but not with TIF1α or TIF1γ. Taken together, these results implicate TIF1β as a common transcriptional corepressor for the three distinct subfamilies of KRAB zinc-finger proteins and suggest a high degree of conservation in the molecular mechanism underlying their transcriptional repression activity.

A thermostable endonuclease III homolog from the archaeon Pyrobaculum aerophilum

Pyrimidine adducts in cellular DNA arise from modification of the pyrimidine 5,6-double bond by oxidation, reduction or hydration. The biological outcome includes increased mutation rate and potential lethality. A major DNA N-glycosylase responsible for the excision of modified pyrimidine bases is the base excision repair (BER) glycosylase endonuclease III, for which functional homologs have been identified and characterized in Escherichia coli, yeast and humans. So far, little is known about how hyperthermophilic Archaea cope with such pyrimidine damage. Here we report characterization of an endonuclease III homolog, PaNth, from the hyperthermophilic archaeon Pyrobaculum aerophilum, whose optimal growth temperature is 100°C. The predicted product of 223 amino acids shares significant sequence homology with several [4Fe-4S]-containing DNA N-glycosylases including E.coli endonuclease III (EcNth). The histidine-tagged recombinant protein was expressed in E.coli and purified. Under optimal conditions of 80–160 mM NaCl and 70°C, PaNth displays DNA glycosylase/β-lyase activity with the modified pyrimidine base 5,6-dihydrothymine (DHT). This activity is enhanced when DHT is paired with G. Our data, showing the structural and functional similarity between PaNth and EcNth, suggests that BER of modified pyrimidines may be a conserved repair mechanism in Archaea. Conserved amino acid residues are identified for five subfamilies of endonuclease III/UV endonuclease homologs clustered by phylogenetic analysis.

Vipp1 deletion mutant of Synechocystis: A connection between bacterial phage shock and thylakoid biogenesis?

Plant chloroplasts originated from an endosymbiotic event by which an ancestor of contemporary cyanobacteria was engulfed by an early eukaryotic cell and then transformed into an organelle. Oxygenic photosynthesis is the specific feature of cyanobacteria and chloroplasts, and the photosynthetic machinery resides in an internal membrane system, the thylakoids. The origin and genesis of thylakoid membranes, which are essential for oxygenic photosynthesis, are still an enigma. Vipp1 (vesicle-inducing protein in plastids 1) is a protein located in both the inner envelope and the thylakoids of Pisum sativum and Arabidopsis thaliana. In Arabidopsis disruption of the VIPP1 gene severely affects the plant's ability to form properly structured thylakoids and as a consequence to carry out photosynthesis. In contrast, Vipp1 in Synechocystis appears to be located exclusively in the plasma membrane. Yet, as in higher plants, disruption of the VIPP1 gene locus leads to the complete loss of thylakoid formation. So far VIPP1 genes are found only in organisms carrying out oxygenic photosynthesis. They share sequence homology with a subunit encoded by the bacterial phage shock operon (PspA) but differ from PspA by a C-terminal extension of about 30 amino acids. In two cyanobacteria, Synechocystis and Anabaena, both a VIPP1 and a pspA gene are present, and phylogenetic analysis indicates that VIPP1 originated from a gene duplication of the latter and thereafter acquired its new function. It also appears that the C-terminal extension that discriminates VIPP1 proteins from PspA is important for its function in thylakoid formation.

X-ray structure of the human hyperplastic discs protein: An ortholog of the C-terminal domain of poly(A)-binding protein

The poly(A)-binding protein (PABP) recognizes the 3′ mRNA poly(A) tail and plays an essential role in eukaryotic translation initiation and mRNA stabilization/degradation. PABP is a modular protein, with four N-terminal RNA-binding domains and an extensive C terminus. The C-terminal region of PABP is essential for normal growth in yeast and has been implicated in mediating PABP homo-oligomerization and protein–protein interactions. A small, proteolytically stable, highly conserved domain has been identified within this C-terminal segment. Remarkably, this domain is also present in the hyperplastic discs protein (HYD) family of ubiquitin ligases. To better understand the function of this conserved region, an x-ray structure of the PABP-like segment of the human HYD protein has been determined at 1.04-Å resolution. The conserved domain adopts a novel fold resembling a right-handed supercoil of four α-helices. Sequence profile searches and comparative protein structure modeling identified a small ORF from the Arabidopsis thaliana genome that encodes a structurally similar but distantly related PABP/HYD domain. Phylogenetic analysis of the experimentally determined (HYD) and homology modeled (PABP) protein surfaces revealed a conserved feature that may be responsible for binding to a PABP interacting protein, Paip1, and other shared interaction partners.

Characterization of Two Divergent Endo-β-1,4-Glucanase cDNA Clones Highly Expressed in the Nonclimacteric Strawberry Fruit

Two cDNAs clones (Cel1 and Cel2) encoding divergent endo-β-1,4-glucanases (EGases) have been isolated from a cDNA library obtained from ripe strawberry (Fragaria x ananassa Duch) fruit. The analysis of the amino acid sequence suggests that Cel1 and Cel2 EGases have different secondary and tertiary structures and that they differ in the presence of potential N-glycosylation sites. By in vitro translation we show that Cel1 and Cel2 bear a functional signal peptide, the cleavage of which yields mature proteins of 52 and 60 kD, respectively. Phylogenetic analysis revealed that the Cel2 EGase diverged early in evolution from other plant EGases. Northern analysis showed that both EGases are highly expressed in fruit and that they have different temporal patterns of accumulation. The Cel2 EGase was expressed in green fruit, accumulating as the fruit turned from green to white and remaining at an elevated, constant level throughout fruit ripening. In contrast, the Cel1 transcript was not detected in green fruit and only a low level of expression was observed in white fruit. The level of Cel1 mRNA increased gradually during ripening, reaching a maximum in fully ripe fruit. The high levels of Cel1 and Cel2 mRNA in ripe fruit and their overlapping patterns of expression suggest that these EGases play an important role in softening during ripening. In addition, the early expression of Cel2 in green fruit, well before significant softening begins, suggests that the product of this gene may also be involved in processes other than fruit softening, e.g. cell wall expansion.

A Millennial Myosin Census

The past decade has seen a remarkable explosion in our knowledge of the size and diversity of the myosin superfamily. Since these actin-based motors are candidates to provide the molecular basis for many cellular movements, it is essential that motility researchers be aware of the complete set of myosins in a given organism. The availability of cDNA and/or draft genomic sequences from humans, Drosophila melanogaster, Caenorhabditis elegans, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Dictyostelium discoideum has allowed us to tentatively define and compare the sets of myosin genes in these organisms. This analysis has also led to the identification of several putative myosin genes that may be of general interest. In humans, for example, we find a total of 40 known or predicted myosin genes including two new myosins-I, three new class II (conventional) myosins, a second member of the class III/ninaC myosins, a gene similar to the class XV deafness myosin, and a novel myosin sharing at most 33% identity with other members of the superfamily. These myosins are in addition to the recently discovered class XVI myosin with N-terminal ankyrin repeats and two human genes with similarity to the class XVIII PDZ-myosin from mouse. We briefly describe these newly recognized myosins and extend our previous phylogenetic analysis of the myosin superfamily to include a comparison of the complete or nearly complete inventories of myosin genes from several experimentally important organisms.