Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.

α10: A determinant of nicotinic cholinergic receptor function in mammalian vestibular and cochlear mechanosensory hair cells

We report the cloning and characterization of rat α10, a previously unidentified member of the nicotinic acetylcholine receptor (nAChR) subunit gene family. The protein encoded by the α10 nAChR subunit gene is most similar to the rat α9 nAChR, and both α9 and α10 subunit genes are transcribed in adult rat mechanosensory hair cells. Injection of Xenopus laevis oocytes with α10 cRNA alone or in pairwise combinations with either α2-α6 or β2-β4 subunit cRNAs yielded no detectable ACh-gated currents. However, coinjection of α9 and α10 cRNAs resulted in the appearance of an unusual nAChR subtype. Compared with homomeric α9 channels, the α9α10 nAChR subtype displays faster and more extensive agonist-mediated desensitization, a distinct current–voltage relationship, and a biphasic response to changes in extracellular Ca2+ ions. The pharmacological profiles of homomeric α9 and heteromeric α9α10 nAChRs are essentially indistinguishable and closely resemble those reported for endogenous cholinergic eceptors found in vertebrate hair cells. Our data suggest that efferent modulation of hair cell function occurs, at least in part, through heteromeric nAChRs assembled from both α9 and α10 subunits.

Fusicoccin, 14-3-3 Proteins, and Defense Responses in Tomato Plants1

Fusicoccin (FC) is a fungal toxin that activates the plant plasma membrane H+-ATPase by binding with 14-3-3 proteins, causing membrane hyperpolarization. Here we report on the effect of FC on a gene-for-gene pathogen-resistance response and show that FC application induces the expression of several genes involved in plant responses to pathogens. Ten members of the FC-binding 14-3-3 protein gene family were isolated from tomato (Lycopersicon esculentum) to characterize their role in defense responses. Sequence analysis is suggestive of common biochemical functions for these tomato 14-3-3 proteins, but their genes showed different expression patterns in leaves after challenges. Different specific subsets of 14-3-3 genes were induced after treatment with FC and during a gene-for-gene resistance response. Possible roles for the H+-ATPase and 14-3-3 proteins in responses to pathogens are discussed.

Molecular Analysis of (R)-(+)-Mandelonitrile Lyase Microheterogeneity in Black Cherry1

The flavoprotein (R)-(+)-mandelonitrile lyase (MDL; EC 4.1.2.10), which plays a key role in cyanogenesis in rosaceous stone fruits, occurs in black cherry (Prunus serotina Ehrh.) homogenates as several closely related isoforms. Biochemical and molecular biological methods were used to investigate MDL microheterogeneity and function in this species. Three novel MDL cDNAs of high sequence identity (designated MDL2, MDL4, and MDL5) were isolated. Like MDL1 and MDL3 cDNAs (Z. Hu, J.E. Poulton [1997] Plant Physiol 115: 1359–1369), they had open reading frames that predicted a flavin adenine dinucleotide-binding site, multiple N-glycosylation sites, and an N-terminal signal sequence. The N terminus of an MDL isoform purified from seedlings matched the derived amino acid sequence of the MDL4 cDNA. Genomic sequences corresponding to the MDL1, MDL2, and MDL4 cDNAs were obtained by polymerase chain reaction amplification of genomic DNA. Like the previously reported mdl3 gene, these genes are interrupted at identical positions by three short, conserved introns. Given their overall similarity, we conclude that the genes mdl1, mdl2, mdl3, mdl4, and mdl5 are derived from a common ancestral gene and constitute members of a gene family. Genomic Southern-blot analysis showed that this family has approximately eight members. Northern-blot analysis using gene-specific probes revealed differential expression of the genes mdl1, mdl2, mdl3, mdl4, and mdl5.

A Photoperiod-Insensitive Barley Line Contains a Light-Labile Phytochrome B1

Barley (Hordeum vulgare L.) is a long-day plant whose flowering is enhanced when the photoperiod is supplemented with far-red light, and this promotion is mediated by phytochrome. A chemically mutagenized dwarf cultivar of barley was selected for early flowering time (barley maturity daylength response [BMDR]-1) and was made isogenic with the cultivar Shabet (BMDR-8) by backcrossing. BMDR-1 was found to contain higher levels of both phytochrome A and phytochrome B in the dark on immunoblots with monoclonal antibodies from oat (Avena sativa L.) that are specific to different members of the phytochrome gene family. Phytochrome A was light labile in both BMDR-1 and BMDR-8, decreasing to very low levels after 4 d of growth in the light. Phytochrome B was light stable in BMDR-8, being equal in both light and darkness. However, phytochrome B became light labile in BMDR-1 and this destabilization of phytochrome B appeared to make BMDR-1 insensitive to photoperiod. In addition, both the mutant and the wild type lacked any significant promotion of flowering in response to a pulse of far-red light given at the end of day, and the end-of-day, far-red inhibition of tillering is normal in both, suggesting that phytochrome B is not involved with these responses in barley.

Arginase Is Inoperative in Developing Soybean Embryos1

Arginase (EC 3.5.3.1) transcript level and activity were measured in soybean (Glycine max L.) embryos from the reserve deposition stage to postgermination. Using a cDNA probe for a small soybean arginase gene family, no transcript was detected in developing embryos. However, arginase transcripts increased sharply on germination, reaching a maximum at 3 to 5 d after germination. There was low but measurable in vitro arginase specific activity in developing embryos (less than 6% of seedling maximum). During germination arginase specific activity increased in parallel with the sharply increasing arginase transcript level. Seedling arginase activity was largely localized in cotyledons. Arginase activity was assayed in vivo by measuring urea accumulation in a urease-deficient mutant. No urea was detected in developing embryos, whereas accumulated urea paralleled arginase specific activity and transcript level in germinating seedlings. As in planta embryos, cultured cotyledons did not accumulate urea when arginine (Arg) was provided with other amino acids in a “mock” seed-coat exudate. Arg as the sole nitrogen source was converted to urea but did not support cotyledon growth. There appeared to be a lack of recruitment of the low-level arginase activity to hydrolyze free Arg in developing embryos, thus avoiding a futile urea cycle.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.

RNA editing activity is associated with splicing factors in lnRNP particles: The nuclear pre-mRNA processing machinery

Multiple members of the ADAR (adenosine deaminases acting on RNA) gene family are involved in A-to-I RNA editing. It has been speculated that they may form a large multicomponent protein complex. Possible candidates for such complexes are large nuclear ribonucleoprotein (lnRNP) particles. The lnRNP particles consist mainly of four spliceosomal subunits that assemble together with the pre-mRNA to form a large particle and thus are viewed as the naturally assembled pre-mRNA processing machinery. Here we investigated the presence of ADARs in lnRNP particles by Western blot analysis using anti-ADAR antibodies and by indirect immunoprecipitation. Both ADAR1 and ADAR2 were found associated with the spliceosomal components Sm and SR proteins within the lnRNP particles. The two ADARs, associated with lnRNP particles, were enzymatically active in site-selective A-to-I RNA editing. We demonstrate the association of ADAR RNA editing enzymes with physiological supramolecular complexes, the lnRNP particles.

Genomic analysis of orthologous mouse and human olfactory receptor loci

Olfactory receptor (OR) genes represent ≈1% of genomic coding sequence in mammals, and these genes are clustered on multiple chromosomes in both the mouse and human genomes. We have taken a comparative genomics approach to identify features that may be involved in the dynamic evolution of this gene family and in the transcriptional control that results in a single OR gene expressed per olfactory neuron. We sequenced ≈350 kb of the murine P2 OR cluster and used synteny, gene linkage, and phylogenetic analysis to identify and sequence ≈111 kb of an orthologous cluster in the human genome. In total, 18 mouse and 8 human OR genes were identified, including 7 orthologs that appear to be functional in both species. Noncoding homology is evident between orthologs and generally is confined within the transcriptional unit. We find no evidence for common regulatory features shared among paralogs, and promoter regions generally do not contain strong promoter motifs. We discuss these observations, as well as OR clustering, in the context of evolutionary expansion and transcriptional regulation of OR repertoires.

Hexose Transport in Growing Petunia Pollen Tubes and Characterization of a Pollen-Specific, Putative Monosaccharide Transporter1

We investigated the molecular and physiological processes of sugar uptake and metabolism during pollen tube growth and plant fertilization. In vitro germination assays showed that petunia (Petunia hybrida) pollen can germinate and grow not only in medium containing sucrose (Suc) as a carbon source, but also in medium containing the monosaccharides glucose (Glc) or fructose (Fru). Furthermore, high-performance liquid chromatography analysis demonstrated a rapid and complete conversion of Suc into equimolar amounts of Glc and Fru when pollen was cultured in a medium containing 2% Suc. This indicates the presence of wall-bound invertase activity and uptake of sugars in the form of monosaccharides by the growing pollen tube. A cDNA designated pmt1 (petunia monosaccharide transporter 1), which is highly homologous to plant monosaccharide transporters, was isolated from petunia. Pmt1 belongs to a small gene family and is expressed specifically in the male gametophyte, but not in any other vegetative or floral tissues. Pmt1 is activated after the first pollen mitosis, and high levels of mRNA accumulate in mature and germinating pollen. A model describing the transport of sugars to the style, the conversion of Suc into Glc and Fru, and the active uptake by a monosaccharide transporter into the pollen tube is presented.

Formation of a Normal Epidermis Supported by Increased Stability of Keratins 5 and 14 in Keratin 10 Null Mice

The expression of distinct keratin pairs during epidermal differentiation is assumed to fulfill specific and essential cytoskeletal functions. This is supported by a great variety of genodermatoses exhibiting tissue fragility because of keratin mutations. Here, we show that the loss of K10, the most prominent epidermal protein, allowed the formation of a normal epidermis in neonatal mice without signs of fragility or wound-healing response. However, there were profound changes in the composition of suprabasal keratin filaments. K5/14 persisted suprabasally at elevated protein levels, whereas their mRNAs remained restricted to the basal keratinocytes. This indicated a novel mechanism regulating keratin turnover. Moreover, the amount of K1 was reduced. In the absence of its natural partner we observed the formation of a minor amount of novel K1/14/15 filaments as revealed by immunogold electron microscopy. We suggest that these changes maintained epidermal integrity. Furthermore, suprabasal keratinocytes contained larger keratohyalin granules similar to our previous K10T mice. A comparison of profilaggrin processing in K10T and K10−/− mice revealed an accumulation of filaggrin precursors in the former but not in the latter, suggesting a requirement of intact keratin filaments for the processing. The mild phenotype of K10−/− mice suggests that there is a considerable redundancy in the keratin gene family.

Cloning of an arylalkylamine N-acetyltransferase (aaNAT1) from Drosophila melanogaster expressed in the nervous system and the gut.

In insects, neurotransmitter catabolism, melatonin precursor formation, and sclerotization involve arylalkylamine N-acetyltransferase (aaNAT, EC 2.3.1.87) activity. It is not known if one or multiple aaNAT enzymes are responsible for these activities. We recently have purified an aaNAT from Drosophila melanogaster. Here, we report the cloning of the corresponding aaNAT cDNA (aaNAT1) that upon COS cell expression acetylates dopamine, tryptamine, and the immediate melatonin precursor serotonin. aaNAT1 represents a novel gene family unrelated to known acetyl-transferases, except in two weakly conserved amino acid motifs. In situ hybridization studies of aaNAT1 mRNA in embryos reveal hybridization signals in the brain, the ventral cord, the gut, and probably in oenocytes, indicating a broad tissue distribution of aaNAT1 transcripts. Moreover, in day/ night studies we demonstrate a diurnal rhythm of melatonin concentration without a clear-cut change in aaNAT1 mRNA levels. The data suggest that tissue-specific regulation of aaNAT1 may be associated with different enzymatic functions and do not exclude the possibility of additional aaNAT genes.

Higher plants contain homologs of the bacterial celA genes encoding the catalytic subunit of cellulose synthase.

In spite of much effort, no one has succeeded in isolating and characterizing the enzyme(s) responsible for synthesis of cellulose, the major cell wall polymer of plants. We have characterized two cotton (Gossypium hirsutum) cDNA clones and identified one rice (Oryza sativa) cDNA that are homologs of the bacterial celA genes that encode the catalytic subunit of cellulose synthase. Three regions in the deduced amino acid sequences of the plant celA gene products are conserved with respect to the proteins encoded by bacterial celA genes. Within these conserved regions, there are four highly conserved subdomains previously suggested to be critical for catalysis and/or binding of the substrate UDP-glucose (UDP-Glc). An overexpressed DNA segment of the cotton celA1 gene encodes a polypeptide fragment that spans these domains and binds UDP-Glc, while a similar fragment having one of these domains deleted does not. The plant celA genes show little homology at the N- and C-terminal regions and also contain two internal insertions of sequence, one conserved and one hypervariable, that are not found in the bacterial gene sequences. Cotton celA1 and celA2 genes are expressed at high levels during active secondary wall cellulose synthesis in developing cotton fibers. Genomic Southern blot analyses in cotton demonstrate that celA forms a small gene family.

Evolution of alcohol dehydrogenase genes in the palm and grass families.

The alcohol dehydrogenase (Adh; alcohol:NAD+ oxidoreductase, EC 1.1.1.1) gene family has two or three loci in a broad array of angiosperm species. The relative stability in the number of Adh loci led Gottlieb [Gottlieb, L. D. (1982) Science 216, 373-380] to propose that the Adh gene family arose from an ancient gene duplication. In this study, the isolation of three loci from the California fan palm (Washingtonia robusta) is reported. The three loci from palm are highly diverged. One palm Adh gene, referred to here as adhB, has been completely sequenced, including 950 nucleotides of the upstream regulatory region. For the second locus, adhA, 81% of the exon sequence is complete. Both show the same basic structure as grass Adh genes in terms of intron number and intron location. The third locus, adhC, for which only a small amount of sequence is available (12% of exon sequence) appears to be more highly diverged. Comparison of the Adh gene families from palms and grasses shows that the adh1 and adh2 genes of grasses, and the adhA and adhB genes of palms, arose by duplication following the divergence of the two families. This finding suggests that the multiple Adh loci in different monocot lineages are not the result of a single ancestral duplication but, rather, of multiple duplication events.

Association of a novel human FE65-like protein with the cytoplasmic domain of the beta-amyloid precursor protein.

We identified a novel human homologue of the rat FE65 gene, hFE65L, by screening the cytoplasmic domain of beta-amyloid precursor protein (beta PP) with the "interaction trap." The cytoplasmic domains of the beta PP homologues, APLP1 and APLP2 (amyloid precursor-like proteins), were also tested for interaction with hFE65L. APLP2, but not APLP1, was found to interact with hFE65L. We confirmed these interactions in vivo by successfully coimmunoprecipatating endogenous beta PP and APLP2 from mammalian cells overexpressing a hemagglutinin-tagged fusion of the C-terminal region of hFE65L. We report the existence of a human FE65 gene family and evidence supporting specific interactions between members of the beta PP and FE65 protein families. Sequence analysis of the FE65 human gene family reveals the presence of two phosphotyrosine interaction (PI) domains. Our data show that a single PI domain is sufficient for binding of hFE65L to the cytoplasmic domain of beta PP and APLP2. The PI domain of the protein, Shc, is known to interact with the NPXYp motif found in the cytoplasmic domain of a number of different growth factor receptors. Thus, it is likely that the PI domains present in the C-terminal moiety of the hFE65L protein bind the NPXY motif located in the cytoplasmic domain of beta PP and APLP2.