Origin and evolutionary relationships of giant Galápagos tortoises

Perhaps the most enduring debate in reptile systematics has involved the giant Galápagos tortoises (Geochelone nigra), whose origins and systematic relationships captivated Charles Darwin and remain unresolved to this day. Here we report a phylogenetic reconstruction based on mitochondrial DNA sequences from Galápagos tortoises and Geochelone from mainland South America and Africa. The closest living relative to the Galápagos tortoise is not among the larger-bodied tortoises of South America but is the relatively small-bodied Geochelone chilensis, or Chaco tortoise. The split between G. chilensis and the Galápagos lineage probably occurred 6 to 12 million years ago, before the origin of the oldest extant Galápagos island. Our data suggest that the four named southern subspecies on the largest island, Isabela, are not distinct genetic units, whereas a genetically distinct northernmost Isabela subspecies is probably the result of a separate colonization. Most unexpectedly, the lone survivor of the abingdoni subspecies from Pinta Island (“Lonesome George”) is very closely related to tortoises from San Cristóbal and Española, the islands farthest from the island of Pinta. To rule out a possible recent transplant of Lonesome George, we sequenced DNA from three tortoises collected on Pinta in 1906. They have sequences identical to Lonesome George, consistent with his being the last survivor of his subspecies. This finding may provide guidance in finding a mate for Lonesome George, who so far has failed to reproduce.

The diversity and evolutionary relationships of the pregnancy-associated glycoproteins, an aspartic proteinase subfamily consisting of many trophoblast-expressed genes

The pregnancy-associated glycoproteins (PAGs) are structurally related to the pepsins, thought to be restricted to the hooved (ungulate) mammals and characterized by being expressed specifically in the outer epithelial cell layer (chorion/trophectoderm) of the placenta. At least some PAGs are catalytically inactive as proteinases, although each appears to possess a cleft capable of binding peptides. By cloning expressed genes from ovine and bovine placental cDNA libraries, by Southern genomic blotting, by screening genomic libraries, and by using PCR to amplify portions of PAG genes from genomic DNA, we estimate that cattle, sheep, and most probably all ruminant Artiodactyla possess many, possibly 100 or more, PAG genes, many of which are placentally expressed. The PAGs are highly diverse in sequence, with regions of hypervariability confined largely to surface-exposed loops. Nonsynonymous (replacement) mutations in the regions of the genes coding for these hypervariable loop segments have accumulated at a higher rate than synonymous (silent) mutations. Construction of distance phylograms, based on comparisons of PAG and related aspartic proteinase amino acid sequences, suggests that much diversification of the PAG genes occurred after the divergence of the Artiodactyla and Perissodactyla, but that at least one gene is represented outside the hooved species. The results also suggest that positive selection of duplicated genes has acted to provide considerable functional diversity among the PAGs, whose presence at the interface between the placenta and endometrium and in the maternal circulation indicates involvement in fetal–maternal interactions.

Conserved structural features of the synaptic fusion complex: SNARE proteins reclassified as Q- and R-SNAREs

SNARE [soluble NSF (N-ethylmaleimide-sensitive fusion protein) attachment protein receptor] proteins are essential for membrane fusion and are conserved from yeast to humans. Sequence alignments of the most conserved regions were mapped onto the recently solved crystal structure of the heterotrimeric synaptic fusion complex. The association of the four α-helices in the synaptic fusion complex structure produces highly conserved layers of interacting amino acid side chains in the center of the four-helix bundle. Mutations in these layers reduce complex stability and cause defects in membrane traffic even in distantly related SNAREs. When syntaxin-4 is modeled into the synaptic fusion complex as a replacement of syntaxin-1A, no major steric clashes arise and the most variable amino acids localize to the outer surface of the complex. We conclude that the main structural features of the neuronal complex are highly conserved during evolution. On the basis of these features we have reclassified SNARE proteins into Q-SNAREs and R-SNAREs, and we propose that fusion-competent SNARE complexes generally consist of four-helix bundles composed of three Q-SNAREs and one R-SNARE.

On the origin of and phylogenetic relationships among living amphibians

The phylogenetic relationships among the three orders of modern amphibians (Caudata, Gymnophiona, and Anura) have been estimated based on both morphological and molecular evidence. Most morphological and paleontological studies of living and fossil amphibians support the hypothesis that salamanders and frogs are sister lineages (the Batrachia hypothesis) and that caecilians are more distantly related. Previous interpretations of molecular data based on nuclear and mitochondrial rRNA sequences suggested that salamanders and caecilians are sister groups to the exclusion of frogs. In an attempt to resolve this apparent conflict, the complete mitochondrial genomes of a salamander (Mertensiella luschani) and a caecilian (Typhlonectes natans) were determined (16,656 and 17,005 bp, respectively) and compared with previously published sequences from a frog (Xenopus laevis) and several other groups of vertebrates. Phylogenetic analyses of the mitochondrial data supported with high bootstrap values the monophyly of living amphibians with respect to other living groups of tetrapods, and a sister group relationship of salamanders and frogs. The lack of phylogenetically informative sites in the previous rRNA data sets (because of its shorter size and higher among-site rate variation) likely explains the discrepancy between our results and those based on previous molecular data. Strong support of the Batrachia hypothesis from both molecule- and morphology-based studies provides a robust phylogenetic framework that will be helpful to comparative studies among the three living orders of amphibians and will permit better understanding of the considerably divergent vertebral, brain, and digit developmental patterns found in frogs and salamanders.

Crossover and noncrossover recombination during meiosis: timing and pathway relationships.

During meiosis, crossovers occur at a high level, but the level of noncrossover recombinants is even higher. The biological rationale for the existence of the latter events is not known. It has been suggested that a noncrossover-specific pathway exists specifically to mediate chromosome pairing. Using a physical assay that monitors both crossovers and noncrossovers in cultures of yeast undergoing synchronous meiosis, we find that both types of products appear at essentially the same time, after chromosomes are fully synapsed at pachytene. We have also analyzed a situation in which commitment to meiotic recombination and formation of the synaptonemal complex are coordinately suppressed (mer1 versus mer1 MER2++). We find that suppression is due primarily to restoration of meiosis-specific double-strand breaks, a characteristic of the major meiotic recombination pathway. Taken together, the observations presented suggest that there probably is no noncrossover-specific pathway and that restoration of intermediate events in a single pairing/recombination pathway promotes synaptonemal complex formation. The biological significant of noncrossover recombination remains to be determined, however.

The small GTP-binding protein Rho links G protein-coupled receptors and Gα12 to the serum response element and to cellular transformation

Receptors coupled to heterotrimeric G proteins can effectively stimulate growth promoting pathways in a large variety of cell types, and if persistently activated, these receptors can also behave as dominant-acting oncoproteins. Consistently, activating mutations for G proteins of the Gαs and Gαi2 families were found in human tumors; and members of the Gαq and Gα12 families are fully transforming when expressed in murine fibroblasts. In an effort aimed to elucidate the molecular events involved in proliferative signaling through heterotrimeric G proteins we have focused recently on gene expression regulation. Using NIH 3T3 fibroblasts expressing m1 muscarinic acetylcholine receptors as a model system, we have observed that activation of this transforming G protein-coupled receptors induces the rapid expression of a variety of early responsive genes, including the c-fos protooncogene. One of the c-fos promoter elements, the serum response element (SRE), plays a central regulatory role, and activation of SRE-dependent transcription has been found to be regulated by several proteins, including the serum response factor and the ternary complex factor. With the aid of reporter plasmids for gene expression, we observed here that stimulation of m1 muscarinic acetylcholine receptors potently induced SRE-driven reporter gene activity in NIH 3T3 cells. In these cells, only the Gα12 family of heterotrimeric G protein α subunits strongly induced the SRE, while Gβ1γ2 dimers activated SRE to a more limited extent. Furthermore, our study provides strong evidence that m1, Gα12 and the small GTP-binding protein RhoA are components of a novel signal transduction pathway that leads to the ternary complex factor-independent transcriptional activation of the SRE and to cellular transformation.

Rhodopsin kinase: Expression in baculovirus-infected insect cells, and characterization of post-translational modifications*

Structure–function studies of rhodopsin kinase (RK; EC 2.7.1.125) require a variety of mutants. Therefore, there is need for a suitable system for the expression of RK mutant genes. Here we report on a study of expression of the RK gene in baculovirus-infected Sf21 cells and characterization of the enzyme produced as purified to near homogeneity. Particular attention has been paid to the post-translational modifications, autophosphorylation and isoprenylation, found in the native bovine RK. The protein produced has been purified using, successively, heparin-Sepharose, Mono Q, and Mono S FPLC (fast protein liquid chromatography) and was obtained in amounts of about 2 mg from 1 liter of cell culture. The enzyme from the last step of purification was obtained in two main fractions that differ in the level of phosphorylation. The protein peak eluted first carries two phosphate groups per protein, whereas the second protein peak is monophosphorylated. Further, while both peaks are isoprenylated, the isoprenyl groups consist of mixtures of C5, C10, C15, and C20 isoprenyl moieties. From these results, we conclude that the above expression system is suitable for some but not all aspects of structure–function studies.

Isolation of a brefeldin A-inhibited guanine nucleotide-exchange protein for ADP ribosylation factor (ARF) 1 and ARF3 that contains a Sec7-like domain

Brefeldin A (BFA) inhibited the exchange of ADP ribosylation factor (ARF)-bound GDP for GTP by a Golgi-associated guanine nucleotide-exchange protein (GEP) [Helms, J. B. & Rothman, J. E. (1992) Nature (London) 360, 352–354; Donaldson, J. G., Finazzi, D. & Klausner, R. D. (1992) Nature (London) 360, 350–352]. Cytosolic ARF GEP was also inhibited by BFA, but after purification from bovine brain and rat spleen, it was no longer BFA-sensitive [Tsai, S.-C., Adamik, R., Moss, J. & Vaughan, M. (1996) Proc. Natl. Acad. Sci. USA 93, 305–309]. We describe here purification from bovine brain cytosol of a BFA-inhibited GEP. After chromatography on DEAE–Sephacel, hydroxylapatite, and Mono Q and precipitation at pH 5.8, GEP was eluted from Superose 6 as a large molecular weight complex at the position of thyroglobulin (≈670 kDa). After SDS/PAGE of samples from column fractions, silver-stained protein bands of ≈190 and 200 kDa correlated with activity. BFA-inhibited GEP activity of the 200-kDa protein was demonstrated following electroelution from the gel and renaturation by dialysis. Four tryptic peptides from the 200-kDa protein had amino acid sequences that were 47% identical to sequences in Sec7 from Saccharomyces cerevisiae (total of 51 amino acids), consistent with the view that the BFA-sensitive 200-kDa protein may be a mammalian counterpart of Sec7 that plays a similar role in cellular vesicular transport and Sec7 may be a GEP for one or more yeast ARFs.

Major recent and independent changes in levels and patterns of expression have occurred at the b gene, a regulatory locus in maize

The b locus encodes a transcription factor that regulates the expression of genes that produce purple anthocyanin pigment. Different b alleles are expressed in distinct tissues, causing tissue-specific anthocyanin production. Understanding how phenotypic diversity is produced and maintained at the b locus should provide models for how other regulatory genes, including those that influence morphological traits and development, evolve. We have investigated how different levels and patterns of pigmentation have evolved by determining the phenotypic and evolutionary relationships between 18 alleles that represent the diversity of b alleles in Zea mays. Although most of these alleles have few phenotypic differences, five alleles have very distinct tissue-specific patterns of pigmentation. Superimposing the phenotypes on the molecular phylogeny reveals that the alleles with strong and distinctive patterns of expression are closely related to alleles with weak expression, implying that the distinctive patterns have arisen recently. We have identified apparent insertions in three of the five phenotypically distinct alleles, and the fourth has unique upstream restriction fragment length polymorphisms relative to closely related alleles. The insertion in B-Peru has been shown to be responsible for its unique expression and, in the other two alleles, the presence of the insertion correlates with the phenotype. These results suggest that major changes in gene expression are probably the result of large-scale changes in DNA sequence and/or structure most likely mediated by transposable elements.

GAIP is membrane-anchored by palmitoylation and interacts with the activated (GTP-bound) form of Gαi subunits

GAIP (G Alpha Interacting Protein) is a member of the recently described RGS (Regulators of G-protein Signaling) family that was isolated by interaction cloning with the heterotrimeric G-protein Gαi3 and was recently shown to be a GTPase-activating protein (GAP). In AtT-20 cells stably expressing GAIP, we found that GAIP is membrane-anchored and faces the cytoplasm, because it was not released by sodium carbonate treatment but was digested by proteinase K. When Cos cells were transiently transfected with GAIP and metabolically labeled with [35S]methionine, two pools of GAIP—a soluble and a membrane-anchored pool—were found. Since the N terminus of GAIP contains a cysteine string motif and cysteine string proteins are heavily palmitoylated, we investigated the possibility that membrane-anchored GAIP might be palmitoylated. We found that after labeling with [3H]palmitic acid, the membrane-anchored pool but not the soluble pool was palmitoylated. In the yeast two-hybrid system, GAIP was found to interact specifically with members of the Gαi subfamily, Gαi1, Gαi2, Gαi3, Gαz, and Gαo, but not with members of other Gα subfamilies, Gαs, Gαq, and Gα12/13. The C terminus of Gαi3 is important for binding because a 10-aa C-terminal truncation and a point mutant of Gαi3 showed significantly diminished interaction. GAIP interacted preferentially with the activated (GTP) form of Gαi3, which is in keeping with its GAP activity. We conclude that GAIP is a membrane-anchored GAP with a cysteine string motif. This motif, present in cysteine string proteins found on synaptic vesicles, pancreatic zymogen granules, and chromaffin granules, suggests GAIP’s possible involvement in membrane trafficking.

Spatial memory is related to hippocampal subcellular concentrations of calcium-dependent protein kinase C isoforms in young and aged rats

Relationships were examined between spatial learning and hippocampal concentrations of the α, β2, and γ isoforms of protein kinase C (PKC), an enzyme implicated in neuronal plasticity and memory formation. Concentrations of PKC were determined for individual 6-month-old (n = 13) and 24-month-old (n = 27) male Long–Evans rats trained in the water maze on a standard place-learning task and a transfer task designed for rapid acquisition. The results showed significant relationships between spatial learning and the amount of PKC among individual subjects, and those relationships differed according to age, isoform, and subcellular fraction. Among 6-month-old rats, those with the best spatial memory were those with the highest concentrations of PKCγ in the particulate fraction and of PKCβ2 in the soluble fraction. Aged rats had increased hippocampal PKCγ concentrations in both subcellular fractions in comparison with young rats, and memory impairment was correlated with higher PKCγ concentrations in the soluble fraction. No age difference or correlations with behavior were found for concentrations of PKCγ in a comparison structure, the neostriatum, or for PKCα in the hippocampus. Relationships between spatial learning and hippocampal concentrations of calcium-dependent PKC are isoform-specific. Moreover, age-related spatial memory impairment is associated with altered subcellular concentrations of PKCγ and may be indicative of deficient signal transduction and neuronal plasticity in the hippocampal formation.

PlasmoDB: An integrative database of the Plasmodium falciparum genome. Tools for accessing and analyzing finished and unfinished sequence data

The Plasmodium falciparum Genome Database (http://PlasmoDB.org) integrates sequence information, automated analyses and annotation data emerging from the P.falciparum genome sequencing consortium. To date, raw sequence coverage is available for >90% of the genome, and two chromosomes have been finished and annotated. Data in PlasmoDB are organized by chromosome (1–14), and can be accessed using a variety of tools for graphical and text-based browsing or downloaded in various file formats. The GUS (Genomics Unified Schema) implementation of PlasmoDB provides a multi-species genomic relational database, incorporating data from human and mouse, as well as P.falciparum. The relational schema uses a highly structured format to accommodate diverse data sets related to genomic sequence and gene expression. Tools have been designed to facilitate complex biological queries, including many that are specific to Plasmodium parasites and malaria as a disease. Additional projects seek to integrate genomic information with the rich data sets now becoming available for RNA transcription, protein expression, metabolic pathways, genetic and physical mapping, antigenic and population diversity, and phylogenetic relationships with other apicomplexan parasites. The overall goal of PlasmoDB is to facilitate Internet- and CD-ROM-based access to both finished and unfinished sequence information by the global malaria research community.

Gene genealogies and population variation in plants

Early in the development of plant evolutionary biology, genetic drift, fluctuations in population size, and isolation were identified as critical processes that affect the course of evolution in plant species. Attempts to assess these processes in natural populations became possible only with the development of neutral genetic markers in the 1960s. More recently, the application of historically ordered neutral molecular variation (within the conceptual framework of coalescent theory) has allowed a reevaluation of these microevolutionary processes. Gene genealogies trace the evolutionary relationships among haplotypes (alleles) with populations. Processes such as selection, fluctuation in population size, and population substructuring affect the geographical and genealogical relationships among these alleles. Therefore, examination of these genealogical data can provide insights into the evolutionary history of a species. For example, studies of Arabidopsis thaliana have suggested that this species underwent rapid expansion, with populations showing little genetic differentiation. The new discipline of phylogeography examines the distribution of allele genealogies in an explicit geographical context. Phylogeographic studies of plants have documented the recolonization of European tree species from refugia subsequent to Pleistocene glaciation, and such studies have been instructive in understanding the origin and domestication of the crop cassava. Currently, several technical limitations hinder the widespread application of a genealogical approach to plant evolutionary studies. However, as these technical issues are solved, a genealogical approach holds great promise for understanding these previously elusive processes in plant evolution.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Heterogeneity and Photoinhibition of Photosystem II Studied with Thermoluminescence1

Thermoluminescence (TL) signals were recorded from grana stacks, margins, and stroma lamellae from fractionated, dark-adapted thylakoid membranes of spinach (Spinacia oleracea L.) in the absence and in the presence of 2,6-dichlorphenylindophenol (DCMU). In the absence of DCMU, the TL signal from grana fractions consisted of a homogenous B-band, which originates from recombination of the semi-quinone QB− with the S2 state of the water-splitting complex and reflects active photosystem II (PSII). In the presence of DCMU, the B-band was replaced by the Q-band, which originates from an S2QA− recombination. Margin fractions mainly showed two TL-bands, the B- and C-bands, at approximately 50°C in the absence of DCMU, and Q- and C-bands in the presence of DCMU. The C-band is ascribed to a TyrD+-QA− recombination. In the absence of DCMU, the fractions of stromal lamellae mainly gave rise to a TL emission at 42°C. The intensity of this band was independent of the number of excitation flashes and was shifted to higher temperatures (52°C) after the addition of DCMU. Based on these observations, this band was considered to be a C-band. After photoinhibitory light treatment of uncoupled thylakoid membranes, the TL intensities of the B- and Q-bands decreased, whereas the intensity at 45°C (C-band) slightly increased. It is proposed that the 42 to 52°C band that was observed in marginal and stromal lamellae and in photoinhibited thylakoid membranes reflects inactive PSII centers that are assumed to be equivalent to inactive PSII QB-nonreducing centers.