Structure and evolution of 5S rRNA genes and pseudogenes in the genusAspergillus

Summary We have cloned and determined the nucleotide sequence of 18 DNA fragments hybridizing to 5S rRNA from twoAspergillus species-A. wentii andA. awamori. Four of the analyzed sequences were pseudogenes. The gene sequences of these two species were very similar and differed fromAspergillus nidulans at both constant and microheterogeneous sites.     

The legumin gene family: structure and evolutionary implications of Vicia faba B-type genes and pseudogenes

We have characterized several Vicia faba genes encoding methionine residue-free group B subunits of the 11S or legumin storage proteins. The respective gene subfamily consists of 10 to 15 members, six of them having been studied by DNA sequence analysis. Four functional genes (LeB2, LeB4, LeB6, LeB7) are highly homologous in their coding region and 0.3 kb of their 3 flanking sequences. On the other hand, two pseudogenes (LeB1, LeB5) have accumulated a large number of mutations including an identical 0.7 kb internal deletion; they are both flanked by a repetitive element. Analysis of sequence changes show that transitions are nearly double as frequent as transversions. CpG is the most infrequent dinucleotide whereas TpA is significantly underrepresented in exon sequences. End points of deletions are correlated with short direct repeats and preferentially found in the two introns. Our studies indicate that the Vicia faba legumin B gene subfamily contains a group of expressed, highly homologous genes as well as more diverged pseudogenes.     

Nucleotide sequence of the 16S-23S spacer region in the rrnA operon from a blue-green alga, Anacystis nidulans

Summary The nucleotide sequence of an entire spacer region between the 16S and 23S rRNA genes of the rrnA operon from a blue-green alga, Anacystis nidulans, has been determined. The spacer region is 545 base pairs long and encodes tRNA^fle and tRNA^Ala in the order of 16S rRNA-tRNA^fle-tRNA^Ala-23S rRNA. A striking feature is that the A. nidulans tRNA^fle gene contains no 3-CCA sequence while the tRNA^Ala gene does. These spacer tRNA genes show strong sequence homology with those of chloroplasts and bacteria.     

ThePenicillium chrysogenum andAspergillus nidulans wetA developmental regulatory genes are functionally equivalent

Aspergillus nidulans andPenicillium chrysogenum are related fungi that reproduce asexually by forming multicellular conidiophores and uninucleate conidia. InA. nidulans, spore maturation is controlled by thewetA (AwetA) regulatory gene. We cloned a homologous gene (PwetA) fromP. chrysogenum to determine if spore maturation is regulated by a similar mechanism in this species. ThePwetA andAwetA genes are similar in structure and functional organization. The inferred polypeptides share 77% overall amino acid sequence similarity, with several regions having > 85% similarity. The genes also had significant, local sequence similarities in their 5 flanking regions, including conserved binding sites for the product of the regulatory geneabaA.PwetA fully complemented anA. nidulans wetA deletion mutation, demonstrating thatPwetA and its 5 regulatory sequences function normally inA. nidulans. These results indicate that the mechanisms controlling sporulation inA. nidulans andP. chrysogenum are evolutionarily conserved.     

Nucleotide sequence of the rrnB 16S ribosomal RNA gene from Mycoplasma capricolum

Summary The nucleotide sequences of the rrnB 16S ribosomal RNA gene and its 5-and 3-flanking regions from Mycoplasma capricolum have been determined. The coding sequence is 1521 base pairs long, being 21 base pairs shorter than that of the Scherichia coli 16S rRNA gene. The 16S rRNA sequence of M. capricolum reveals 74% and 76% identity with that of E. coli and Anacystis nidulans, respectively. The secondary structure model constructed from the M. capricolum 16S rRNA.gene sequence resembles that proposed for E. coli 16S rRNA. A large stem structure can be constructed between the 5- and 3-flanking sequences of the 16S rRNA gene. The flanking regions are extremely rich in AT.     

Parallel origins of duplications and the formation of pseudogenes in mitochondrial DNA from parthenogenetic lizards (Heteronotia binoei; Gekkonidae)

Summary Analysis of mitochondrial DNAs (mtDNAs) from parthenogenetic lizards of theHeteronotia binoei complex with restriction enzymes revealed an 5-kb addition present in all 77 individuals. Cleavage site mapping suggested the presence of a direct tandem duplication spanning the 16S and 12S rRNA genes, the control region and most, if not all, of the gene for the subunit 1 of NADH dehydrogenase (ND1). The location of the duplication was confirmed by Southern hybridization. A restriction enzyme survey provided evidence for modifications to each copy of the duplicated sequence, including four large deletions. Each gene affected by a deletion was complemented by an intact version in the other copy of the sequence, although for one gene the functional copy was heteroplasmic for another deletion. Sequencing of a fragment from one copy of the duplication which encompassed the tRNA^leu(UUR) and parts of the 16S rRNA and ND1 genes, revealed mutations expected to disrupt function. Thus, evolution subsequent to the duplication event has resulted in mitochondrial pseudogenes. The presence of duplications in all of these parthenogens, but not among representatives of their maternal sexual ancestors, suggests that the duplications arose in the parthenogenetic form. This provides the second instance inH. binoei of mtDNA duplication associated with the transition from sexual to parthenogenetic reproduction. The increased incidence of duplications in parthenogenetic lizards may be caused by errors in mtDNA replication due to either polyploidy or hybridity of their nuclear genomes.     

A computational approach for identifying pseudogenes in the ENCODE regions

Background

Pseudogenes are inheritable genetic elements showing sequence similarity to functional genes but with deleterious mutations. We describe a computational pipeline for identifying them, which in contrast to previous work explicitly uses intron-exon structure in parent genes to classify pseudogenes. We require alignments between duplicated pseudogenes and their parents to span intron-exon junctions, and this can be used to distinguish between true duplicated and processed pseudogenes (with insertions).

Results

Applying our approach to the ENCODE regions, we identify about 160 pseudogenes, 10% of which have clear 'intron-exon' structure and are thus likely generated from recent duplications.

Conclusion

Detailed examination of our results and comparison of our annotation with the GENCODE reference annotation demonstrate that our computation pipeline provides a good balance between identifying all pseudogenes and delineating the precise structure of duplicated genes.

     

Evolution of the plastid ribosomal RNA operon in a nongreen parasitic plant: Accelerated sequence evolution,altered promoter structure,and tRNA pseudogenes

The nucleotide sequence of a 7.4 kb region containing the entire plastid ribosomal RNA operon of the nongreen parasitic plant Epifagus virginiana has been determined. Analysis of the sequence indicates that all four rRNA genes are intact and almost certainly functional. In contrast, the split genes for tRNA^Ile and tRNA^Ala present in the 16S-23S rRNA spacer region have become pseudogenes, and deletion upstream of the 16S rRNA gene has removed a tRNA^Val gene and most of the promoter region for the rRNA operon. The rate of nucleotide substitution in 16S and 23S rRNAs is several times higher in Epifagus than in tobacco, a related photosynthetic plant. Possible reasons for this, including relaxed translational constraints, are discussed.     

Organization of the ribosomal RNA genes in Mycoplasma hyopneumoniae: the 5S rRNA gene is separated from the 16S and 23S rRNA genes

Summary In order to study the organization of the ribosomal RNA genes of Mycoplasma hyopneumoniae the rRNA genes were cloned in phage vectors EMBL3 and EMBL4. By subcloning the restriction fragments into various plasmids and analysing the resulting clones by Southern and Northern blot hybridization, a restriction map of the rRNA genes was generated and the organization of the rRNA genes was determined. The results show that the genes for the 16S and 23S rRNAs are closely spaced and occur only once in the genome, whereas the 5S rRNA gene is separated from the other two genes by more than 4 kb.     

Estimating the intensity of male-driven evolution in rodents by using X-linked and Y-linked Ube 1 genes and pseudogenes

Using sequence data from the last introns of ZFX and ZFY genes, we previously estimated the male-to-female ratio () of mutation rate to be close to 6 in higher primates and 1.8 in rodents. As the mutation rate may vary among different regions of the mammalian genome, it is interesting to see whether sequence data from other regions will give similar estimates. In this study, we have determined the partial genomic sequences of the ubiquitin-activating enzyme El genes (Ube 1x and Ube 1y for the X-linked and Y-linked homologues, respectively) of mice and rats and two mouse Ube 1y pseudogenes. From the intron sequences of the Ube 1 genes, we calculated the divergence of the Y-linked genes (Y = 0.161) and that of the X-linked genes (X = 0.107) between mouse and rat, and found the Y/X ratio to be 1.50. This ratio led to an estimate of = 2.0 with a 95% confidence interval of (1.0, 3.9). Similar estimates of were obtained if mouse Ube 1y pseudogenes were used instead of the mouse Ube 1y functional gene. These estimates are consistent with our previous estimate for rodents and suggest that the sex ratio of mutation rate in rodents is approximately only one-third of that in higher primates. Our estimate of the divergence time between Ube 1x and Ube 1y supports the view that the two genes separated before the eutherian radiation.Correspondence to: W.-H. Li     

Cloning and molecular characterisation of the amdR controlled gatA gene of Aspergillus nidulans

Summary The gamma-amino-n-butyrate transaminase gene (gatA) of Aspergillus nidulans is one of several genes under positive control by the regulatory gene amdR (also called intA). The gatA gene has been cloned from a cosmid library by complementation of a gatA mutation. The sequence of a 2.6 kb genomic fragment containing gatA has been determined. An open reading frame of 1497 bp within this sequences is interrupted by three putative introns and predicts a protein of 55 kDa. Northern analysis confirms control of gatA RNA levels by amdR and also indicates that gatA is not strongly regulated by areA-mediated nitrogen metabolite repression. A. nidulans transformants containing multiple copies of a plasmid carrying an 88 bp fragment from the 5 untranscribed region of gatA grew poorly on substrates whose utilisation is dependent on genes controlled by amdR. This indicated titration of limiting amounts of the amdR gene product by this 88 bp fragment. Comparison of this sequence with the 5 region of the coregulated gene, amdS, reveals probable sites of action for the amdR protein.     

Restriction endonuclease mapping of ribosomal RNA genes: Sequence divergence and the origin of the tetraploid treefrog Hyla versicolor

Hyla chrysoscelis (2n=24) and H. versicolor (2n=48) are a diploid-tetraploid species pair of treefrogs. Restriction endonuclease mapping of ribosomal RNA (rRNA) gene repeat units of diploids collected from eastern and western populations reveals no differences within rRNA gene coding regions but distinctive differences within the nontranscribed spacers. A minimum of two physical maps is required to construct an rRNA gene map for the tetraploid, whose repeat units appear to be a composite, with about 50% of the elements resembling the western diploid population and about 50% resembling the eastern population. These results imply that this population of the tetraploid species may have arisen from a genetically hybrid diploid. Alternatively, the dual level of sequence heterogeneity in H. versicolor may reflect some type of gene flow between the two species. The coding region of the rRNA genes in the tetraploid differs from that in either diploid in about 20% of all repeat units, as exemplified by a BamHI site located near the 5 terminus of the 28 S rRNA gene. If the 20% variant class of 28 S rRNA gene coding sequences is expressed, then there must be two structural classes of ribosomes; if only the 80% sequence class is expressed, then a genetic control mechanism must be capable of distinguishing between the two different sequence variants. It is postulated that the 20% variant sequence class may be correlated with a partial functional diploidization of rRNA genes in the tetraploid species.This research was supported, in part, by NSF Grants CDP-8002341 and PRM-8106947 and by faculty research grants from Miami University to J.C.V.     

Cloning, characterization of the acyl-CoA : 6-amino penicillanic acid acyltransferase gene of Aspergillus nidulans and linkage to the isopenicillin N synthase gene

Summary The penDE gene encoding acyl-CoA:6-amino penicillanic acid acyltransferase (AAT), the last enzyme of the penicillin biosynthetic pathway, has been cloned from the DNA of Aspergillus nidulans. The gene contains three introns which are located in the 5 region of the open reading frame. It encodes a protein of 357 amino acids with a molecular weight of 39 240 Da. The penDE gene of A. nidulans shows 73% similarity at the nucleotide level with the penDE gene of Penicillium chrysogenum. The A. nidulans gene was expressed in P. chrysogenum and complemented the AAT deficiency of the non-producer mutants of P. chrysogenum, npe6 and npe8. The penDE gene of A. nidulans is linked to the pcbC gene, which encodes the isopenicillin N synthase, as also occurs in P. chrysogenum. Both genes show the same orientation and are separated by an intergenic region of 822 nucleotides.     

The ribosomal RNA genes from chloroplasts of mustard (Sinapis alba L.): mapping and sequencing of the leader region

Summary The genes coding for rRNAs from mustard chloroplasts were mapped within the inverted repeat regions of intact ctDNA and on ctDNA fragments cloned in pBR322. R-loop analysis and restriction endonuclease mapping show that the genes for 16S rRNA map at distances of 17 kb from the junctions of the repeat regions with the large unique region. The genes for 23S rRNA are located at distances of 2.8 kb from the junctions with the small unique region. Genes for 4.5S and 5S rRNA are located in close proximity to the 23S rRNA genes towards the small unique region. DNA sequencing of portions of the 5 terminal third from the mustard 16S rRNA gene shows 96–99% homology with the corresponding regions of the maize, tobacco and spinach chloroplast genes. Sequencing of the region proximal to the 16S rRNA gene reveals the presence of a tRNA^Val gene in nearly the same position and with identical sequence as in maize, tobacco and spinach. Somewhat less but still strong homology is also observed for the tDNA ^Val/16S rDNA intercistronic regions and for the regions upstream of the tRNA^Val gene. However, due to many small and also a few larger deletions and insertions in the leader region, common reading frames coding for homologous peptides larger than 44 amino acids can not be detected; it is therefore unlikely that this region contains a protein coding gene.     

Primate DRB6 pseudogenes: clue to the evolutionary origin of the HLA-DR2 haplotype

The HLA-DR2 haplotype contains three \-chain encoding DRB genes and one -chain encoding DRA gene. Of the three DRB genes, two are presumably functional (HLA-DRB1 and HLA-DRB5), whereas the third (HLA-DRBV1) is a pseudogene. A pseudogene closely related to HLA-DRBVI is present in the chimpanzee (Patr-DRB6) and in the gorilla (Gogo-DRB6). We sequenced the HLA-DRBVI and Patr-DRB6 pseudogenes (all exons and most of the introns), and compared the sequence to that of the Gogo-DRB6 gene (of which only the exon sequence is available). All three pseudogenes seem to lack exon 1 and contain other deletions responsible for shifts in the translational reading frame. At least the HLA-DRBVI pseudogene, however, seems to be transcribed nevertheless. The chimpanzee pseudogene contains two inserts in intron 2, one of which is an Alu repeat belonging to the Sb subfamily, while the other remains unidentified. These inserts are lacking in the human gene. A comparison with sequences published by other investigators revealed the presence of the HLA-DRBVI pseudogene also in the DRI and DRw10 haplotypes. Measurements of genetic distances indicate DRB6 to be closely related to the DRB2 pseudogene and to the HLA-DRB4 functional gene. In humans, gorillas, and chimpanzees, the DRB6 pseudogene is associated with the same functional gene (DRB5) indicating that this linkage disequilibrium is at least six million years old and that DR2 is one of the oldest DR haplotypes in higher primates.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M77284-M77295. Address correspondence and offprint requests to: J. Klein.     

Multiple nuclear pseudogenes of mitochondrial DNA exist in the canine genome

Many copies of nuclear counterparts of mitochondrial DNA (mtDNA) were found in nuclear DNA from sperm heads of the domestic dog, Canis familiaris, by DNA-DNA hybridization and DNA sequencing. Nuclear counterparts homologous to the mtDNA D-loop region were cloned into lambda phage vectors (EMBL4 and lambda gt11), and nucleotide sequences of seven different mtDNA pseudogenes were then determined. The seven pseudogenes were E3 (474 bp; 82% homology with canine mtDNA), E13 (1867 bp; 67%), 8B (2375 bp; 78%), 12A (2650 bp; 79%), 33 (4131 bp; 86%), 47 (4251 bp; 86%), and E17 (5721 bp; 71%). These seven mtDNA pseudogenes corresponded to portions of cytoplasmic mtDNA containing the genes ile, ND1, leu, 16S rRNA, val, 12S rRNA, phe, D-loop, pro, thr, cytb, and glu. A neighbor-joining phylogenetic tree constructed from 12S rRNA sequences in mtDNA pseudogenes 8B, 33, 47, and E17 and in 10 mtDNA fragments from other species showed that these four pseudogenes form a monophyletic clade with canine mtDNA. A neighbor-joining phylogenetic tree based on the 318-bp cytb region showed that the canine pseudogenes existed before the divergence of 17 related canids, and their divergence dates were calculated at around 4.4 to 8.6 million years ago.     

Isolation and nucleotide sequence of an autonomously replicating sequence (ARS) element functional in Candida albicans and Saccharomyces cerevisiae

Summary An 8.6-kb fragment was isolated from an EcoRI digest of Candida albicans ATCC 10261 genomic DNA which conferred the property of autonomous replication in Saccharomyces cervisiae on the otherwise non-replicative plasmid pMK155 (5.6 kb). The DNA responsible for the replicative function was subcloned as a 1.2-kb fragment onto a non-replicative plasmid (pRC3915) containing the C. albicans URA3 and LEU2 genes to form plasmid pRC3920. This plasmid was capable of autonomous replication in both S. cerevisiae and C. albicans and transformed S. cerevisiae AH22 (leu2 ^–) to Leu⁺ at a frequency of 2.15 × 10³ transformants per pg DNA, and transformed C. albicans SGY-243 (ura3) to Ura⁺ at a frequency of 1.91 × 10³ transformants per g DNA. Sequence analysis of the cloned DNA revealed the presence of two identical regions of eleven base pairs (5TTTTATGTTTT3) which agreed with the consensus of autonomously replicating sequence (ARS) cores functional in S. cerevisiae. In addition there were two 10/11 and numerous 9/11 matches to the core consensus. The two 11/11 matches to the consensus, CaARS1 and CaARS2, were located on opposite strands in a non-coding AT-rich region and were separated by 107 bp. Also present on the C. albicans DNA, 538 by from the ARS cores, was a gene for 5S rRNA which showed sequence homology with several other yeast 5S rRNA genes. A sub-fragment (494 bp) containing the 5S rRNA gene (but not the region containing the ARS cores) hybridized to genomic DNAs from a number of yeast species, including S. cerevisiae, C. tropicalis, C. pseudotropicalis, C. parapsilosis, C. kruseii, C. (Torulopsis) glabrata and Neurospora crassa. The 709-bp ARS element (but not the 5S rRNA gene) was necessary for high-frequency transformation and autonomous plasmid replication in both S. cerevisiae and C. albicans.EMBL/GenBank database accession number: X16634 (5S rRNA)     

ThehapC gene ofAspergillus nidulans is involved in the expression of CCAAT-containing promoters

The 5 regulatory region of theamdS gene ofAspergillus nidulans, which encodes an acetamidase required for growth on acetamide as a carbon and nitrogen source, contains a CCAAT sequence which is required for setting the basal level ofamdS expression. Mobility shift studies have identified a factor inA. nidulans nuclear extracts which binds to this CCAAT sequence. InSaccharomyces cerevisiae theHAP3 gene encodes one component of a multisubunit complex that binds CCAAT sequences. A search of the EMBL and SwissProt databases has revealed anA. nidulans sequence with significant homology to theHAP3 gene adjacent to the previously cloned regulatory geneamdR. Sequencing of the remainder of this region has confirmed the presence of a gene, designatedhapC, with extensive homology toHAP3. The predicted amino acid sequence of HapC shows extensive identity to HAP3 in the central conserved domain, but shows little conservation in the flanking sequences. A haploid carrying ahapC deletion has been created and is viable, but grows poorly on all media tested. This null mutant grows especially slowly on acetamide as a sole carbon and nitrogen source, indicating thathapC plays a role inamdS expression. In agreement with this notion, it has been shown that thehapC deletion results in reduced levels of expression of anamdS::lacZ reporter gene and this effect is particularly evident under conditions of carbon limitation. Nuclear extracts prepared from thehapC deletion mutant show no CCAAT binding activity to theamdS orgatA promoters, indicating thathapC may encode a component of the complex binding at this sequence.     

A novel class of <Emphasis Type="Italic">Helitron</Emphasis>- related transposable elements in maize contain portions of multiple pseudogenes

We recently described a maize mutant caused by an insertion of a Helitron type transposable element (Lal, S.K., Giroux, M.J., Brendel, V., Vallejos, E. and Hannah, L.C., 2003, Plant Cell, 15: 381–391). Here we describe another Helitron insertion in the barren stalk1 gene of maize. The termini of a 6525 bp insertion in the proximal promoter region of the mutant reference allele of maize barren stalk1 gene (ba1-ref) shares striking similarity to the Helitron insertion we reported in the Shrunken-2 gene. This insertion is embedded with pseudogenes that differ from the pseudogenes discovered in the mutant Shrunken-2 insertion. Using the common terminal ends of the mutant insertions as a query, we discovered other Helitron insertions in maize BAC clones. Based on the comparison of the insertion site and PCR amplified genomic sequences, these elements inserted between AT dinucleotides. These putative non-autonomous Helitroninsertions completely lacked sequences similar to RPA (replication protein A) and DNA Helicases reported in other species. A blastn analysis indicated that both the 5 and 3 termini of Helitrons are repeated in the maize genome. These data provide strong evidence that Helitron type transposable elements are active and may have played an essential role in the evolution and expansion of the maize genome.     

Cloning and characterization of the rDNA repeat unit of Podospora anserina

Summary DNA coding for ribosomal RNA in Podospora anserina has been cloned and was found as a tandemly repeated 8.3 kb sequence. The cloned rDNA was characterized by restriction endonuclease mapping. The location of 5.8S, 18S and 28S rRNA coding regions was established by DNA-RNA hybridization and S1 nuclease mapping. The organization of P. anserina rRNA genes is similar to that of Neurospora crassa and Aspergillus nidulans. The rDNA unit does not contain the sequence coding for 5S RNA.