Protein-coding genes are epigenetically regulated in Arabidopsis polyploids

The fate of redundant genes resulting from genome duplication is poorly understood. Previous studies indicated that ribosomal RNA genes from one parental origin are epigenetically silenced during interspecific hybridization or polyploidization. Regulatory mechanisms for protein-coding genes in polyploid genomes are unknown, partly because of difficulty in studying expression patterns of homologous genes. Here we apply amplified fragment length polymorphism (AFLP)–cDNA display to perform a genome-wide screen for orthologous genes silenced in Arabidopsis suecica, an allotetraploid derived from Arabidopsis thaliana and Cardaminopsis arenosa. We identified ten genes that are silenced from either A. thaliana or C. arenosa origin in A. suecica and located in four of the five A. thaliana chromosomes. These genes represent a variety of RNA and predicted proteins including four transcription factors such as TCP3. The silenced genes in the vicinity of TCP3 are hypermethylated and reactivated by blocking DNA methylation, suggesting epigenetic regulation is involved in the expression of orthologous genes in polyploid genomes. Compared with classic genetic mutations, epigenetic regulation may be advantageous for selection and adaptation of polyploid species during evolution and development.

Measuring genome evolution

The determination of complete genome sequences provides us with an opportunity to describe and analyze evolution at the comprehensive level of genomes. Here we compare nine genomes with respect to their protein coding genes at two levels: (i) we compare genomes as “bags of genes” and measure the fraction of orthologs shared between genomes and (ii) we quantify correlations between genes with respect to their relative positions in genomes. Distances between the genomes are related to their divergence times, measured as the number of amino acid substitutions per site in a set of 34 orthologous genes that are shared among all the genomes compared. We establish a hierarchy of rates at which genomes have changed during evolution. Protein sequence identity is the most conserved, followed by the complement of genes within the genome. Next is the degree of conservation of the order of genes, whereas gene regulation appears to evolve at the highest rate. Finally, we show that some genomes are more highly organized than others: they show a higher degree of the clustering of genes that have orthologs in other genomes.

Evolution of C4 Photosynthesis in Flaveria Species1 : Isoforms of NADP-Malic Enzyme

NADP-malic enzyme (NADP-ME, EC 1.1.1.40), a key enzyme in C4 photosynthesis, provides CO2 to the bundle-sheath chloroplasts, where it is fixed by ribulose-1,5-bisphosphate carboxylase/oxygenase. We characterized the isoform pattern of NADP-ME in different photosynthetic species of Flaveria (C3, C3-C4 intermediate, C4-like, C4) based on sucrose density gradient centrifugation and isoelectric focusing of the native protein, western-blot analysis of the denatured protein, and in situ immunolocalization with antibody against the 62-kD C4 isoform of maize. A 72-kD isoform, present to varying degrees in all species examined, is predominant in leaves of C3 Flaveria spp. and is also present in stem and root tissue. By immunolabeling, NADP-ME was found to be mostly localized in the upper palisade mesophyll chloroplasts of C3 photosynthetic tissue. Two other isoforms of the enzyme, with molecular masses of 62 and 64 kD, occur in leaves of certain intermediates having C4 cycle activity. The 62-kD isoform, which is the predominant highly active form in the C4 species, is localized in bundle-sheath chloroplasts. Among Flaveria spp. there is a 72-kD constitutive form, a 64-kD form that may have appeared during evolution of C4 metabolism, and a 62-kD form that is necessary for the complete functioning of C4 photosynthesis.

Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.

Translocation events in the evolution of aminoacyl-tRNA synthetases.

We have characterized hisS, the gene encoding the histidyl-tRNA synthetase (HisRS) from the tetraodontoid fish Fugu rubripes. The hisS gene is about 3.5 kbp long and contains 13 exons and 12 introns of 172 bp, on average. The Fugu hisS gene encodes a putative protein of 519 amino acids with the three motifs identified as signatures of class 2 aminoacyl-tRNA synthetases. A model for the shifting of intron 8 between Fugu and hamster is proposed based on the successive appearance of a cryptic splicing site followed by an insertion mutation that created a new acceptor site. In addition, sequence comparisons suggest that the hisS gene has undergone a translocation through the first intron. As a result, the Fugu HisRS has an N-terminal sequence markedly different from that in the human and hamster enzymes. We propose that similar events have been responsible for variations at the N-terminal end of other aminoacyl-tRNA synthetases. Our analysis suggests that this involves exchanges through introns of two exons encoding an ancestral 32-amino acid motif.

Identification of functional domains and evolution of Tc1-like transposable elements.

Tc1-like transposable elements from teleost fish have been phylogenetically examined to determine the mechanisms involved in their evolution and conserved domains of function. We identified two new functional domains in these elements. The first is a bipartite nuclear localization signal, indicating that transposons can take advantage of the transport machinery of host cells for nuclear uptake of their transposases. The second is a novel combination of a paired domain-related protein motif juxtaposed to a leucine zipper-like domain located in the putative DNA-binding regions of the transposases. This domain coexists with a special inverted repeat structure in certain transposons in such phylogenetically distant hosts as fish and insects. Our data indicate that reassortment of functional domains and horizontal transmission between species are involved in the formation and spread of new types of transposable elements.

Integration of multiple repeats of geminiviral DNA into the nuclear genome of tobacco during evolution.

Integration of viral DNA into the host nuclear genome, although not unusual in bacterial and animal systems, has surprisingly not been reported for plants. We have discovered geminvirus-related DNA (GRD) sequences, in the form of distinct sets of multiple direct repeats comprising three related repeat classes, situated in a unique locus in the Nicotiana tabacum (tobacco) nuclear genome. The organization of these sequences is similar or identical in eight different tobacco cultivars we have examined. DNA sequence analysis reveals that each repeat has sequences most resembling those of the New World geminiviral DNA replication origin plus the adjacent AL1 gene, encoding the viral replication protein. We believe these GRD sequences originated quite recently in Nicotiana evolution through integration of geminiviral DNA by some combination of the processes of illegitimate recombination, amplification, deletions, and rearrangements. These events must have occurred in plant tissue that was subsequently able to contribute to meristematic tissue yielding gametes. GRD may have been retained in tobacco by selection or by random fixation in a small evolving population. Although we cannot detect transcription of these sequences, this does not exclude the possibility that they may originally have been expressed.

Using in vitro selection to direct the covalent attachment of human immunodeficiency virus type 1 Rev protein to high-affinity RNA ligands.

We have used an in vitro selection procedure called crosslinking SELEX (SELEX = systematic evolution of ligands by exponential enrichment) to identify RNA sequences that bind with high affinity and crosslink to the Rev protein from human immunodeficiency virus type 1 (HIV-1). A randomized RNA library substituted with the photoreactive chromophore 5-iodouracil was irradiated with monochromatic UV light in the presence of Rev. Those sequences with the ability to photocrosslink to Rev were partitioned from the rest of the RNA pool, amplified, and used for the next round of selection. Rounds of photocrosslinking selection were alternated with rounds of selection for RNA sequences with high affinity to Rev. This iterative, dual-selection method yielded RNA molecules with subnanomolar dissociation constants and high efficiency photocrosslinking to Rev. Some of the RNA molecules isolated by this procedure form a stable complex with Rev that is resistant to denaturing gel electrophoresis in the absence of UV irradiation. In vitro selection of nucleic acids by using modified nucleotides allows the isolation of nucleic acid molecules with potentially limitless chemical capacities to covalently attack a target molecule.

The histone fold: a ubiquitous architectural motif utilized in DNA compaction and protein dimerization.

The histones of all eukaryotes show only a low degree of primary structure homology, but our earlier crystallographic results defined a three-dimensional structural motif, the histone fold, common to all core histones. We now examine the specific architectural patterns within the fold and analyze the nature of the amino acid residues within its functional segments. The histone fold emerges as a fundamental protein dimerization motif while the differentiations of the tips of the histone dimers appear to provide the rules of core octamer assembly and the basis for nucleosome regulation. We present evidence for the occurrence of the fold from archaebacteria to mammals and propose the use of this structural motif to define a distinct family of proteins, the histone fold superfamily. It appears that evolution has conserved the conformation of the fold even through variations in primary structure and among proteins with various functional roles.

Myc and Max: molecular evolution of a family of proto-oncogene products and their dimerization partner.

The myc gene family encodes a group of transcription factors that regulate cell proliferation and differentiation. These genes are widely studied because of their importance as proto-oncogenes. Phylogenetic analyses are described here for 45 Myc protein sequences representing c-, N-, L-, S-, and B-myc genes. A gene duplication early in vertebrate evolution produced the c-myc lineage and another lineage that later gave rise to the N- and L-myc lineages by another gene duplication. Evolutionary divergence in the myc gene family corresponds closely to the known branching order of the major vertebrate groups. The patterns of sequence evolution are described for five separate highly conserved regions, and these analyses show that differential rates of sequence divergence (= mosaic evolution) have occurred among conserved motifs. Further, the closely related dimerization partner protein Max exhibits significantly less sequence variability than Myc. It is suggested that the reduced variability in max stems from natural selection acting to preserve dimerization capability with products of myc and related genes.

Shk1, a homolog of the Saccharomyces cerevisiae Ste20 and mammalian p65PAK protein kinases, is a component of a Ras/Cdc42 signaling module in the fission yeast Schizosaccharomyces pombe.

We describe a protein kinase, Shk1, from the fission yeast Schizosaccharomyces pombe, which is structurally related to the Saccharomyces cerevisiae Ste20 and mammalian p65PAK protein kinases. We provide genetic evidence for physical and functional interaction between Shk1 and the Cdc42 GTP-binding protein required for normal cell morphology and mating in S. pombe. We further show that expression of the STE20 gene complements the shk1 null mutation and that Shk1 is capable of signaling to the pheromone-responsive mitogen-activated protein kinase cascade in S. cerevisiae. Our results lead us to propose that signaling modules composed of small GTP-binding proteins and protein kinases related to Shk1, Ste20, and p65PAK, are highly conserved in evolution and participate in both cytoskeletal functions and mitogen-activated protein kinase signaling pathways.

The alpha and gamma 1 isoforms of protein phosphatase 1 are highly and specifically concentrated in dendritic spines.

Protein phosphatase 1 (PP1) is a highly conserved enzyme that has been implicated in diverse biological processes in the brain as well as in nonneuronal tissues. The present study used light and electron microscopic immunocytochemistry to characterize the distribution of two PP1 isoforms, PP1 alpha and PP1 gamma 1, in the rat neostriatum. Both isoforms are heterogeneously distributed in brain with the highest immunoreactivity being found in the neostriatum and hippocampal formation. Further, both isoforms are highly and specifically concentrated in dendritic spines. Weak immunoreactivity is present in dendrites, axons, and some axon terminals. Immunoreactivity for PP1 alpha is also present in the perikaryal cytoplasm and nuclei of most medium- and large-sized neostriatal neurons. The specific localization of PP1 in dendritic spines is consistent with a central role for this enzyme in signal transduction. The data support the concept that, in the course of evolution, spines developed as specialized signal transduction organelles enabling neurons to integrate diverse inputs from multiple afferent nerve terminals.