The Small Ribosomal Subunit RNA Isoforms in Plasmodium Cynomolgi
Total Page:16
File Type:pdf, Size:1020Kb
Copyright 0 1994 by the Genetics Society of America The Small Ribosomal Subunit RNA Isoforms in Plasmodium cynomolgi Vladimir Corredor' and Vincenzo Enea Department of Medical and Molecular Parasitology, New York University Medical Center, New York, New York 1001 0 Manuscript received May 8, 1993 Accepted for publication November 3, 1993 ABSTRACT We report the isolation, characterization and analysis of the small subunit rRNA genes in PZasmodium cynomolgi (Ceylon). As in other Plasmodium species, these genes are present in low copy number, are unlinked and form two types that are distinct in sequence and are expressed stage specifically. The asexually expressed(type A) genes are present in four copies inthe Ceylon- and in five copies inthe Berok' strain. Surprisingly, the sexually expressed (typeB) gene is present in a single copy. The vast majority of the differences between gene types is confined to the variable regions. The pattern of divergence is different from that observed in Plasmodium berghei or in Plasmodium falciparum.Analysis of the small subunit rRNA sequences of P. cynomolgi, P. berghei and P. falciparum,indicates that the two gene types do not evolve independently but rather interact (through gene conversion or some form of recombi- nation) to such an extent as to erase whatever stage-specific sequence signatures they may have had in the last common ancestor. HE organization of the ribosomal RNA genes in Here we report the isolation, characterization and T Plasmodium (the etiological agent of malaria) is analysis ofthe SSU rRNA isoforms in P. cynomolga (Cey- unusual on several counts: the genes arefew in number lon) and theresults of some comparative work on other (4-8 copies per haploid genome), scattered on differ- strains and on other species. ent chromosomes, and, most notably, are organized in two distinct subgroups whose expression is developmen- MATERIALSAND METHODS tally regulated (DAMEand MCCUTCHAN1983; LANGSLEY Parasites: P. cynomolgi parasites werepropagated in rhesus et al. 1983; VANDER PLOEGet al. 1985; GUNDERSONet al. monkeys (Macaca mulatta).Infected red blood cells (IRBC) 1987; WATERSet al. 1989). Type A genes are expressed containing parasites of the Berok and Ceylon strains of P. cy- in the asexual stages (which take place in the vertebrate nomolgi used for pulse field gel (PFG) analysis were kindly host), and type B genes are expressed in the sexual provided by Dr. J. W. Barnwell. stages (which initiate in the vertebrate host and con- Nucleic acids and restriction enzyme digestion analysis: DNA preparations from plasmid, phage and parasites as well tinue in the insect vector) (GUNDERSONet al. 1987; as restriction enzyme digestions, agarose gel fractionation, WATERSet al. 1989). The origin and significance of blotting and hybridizations wereperformed as previously de- this condition is not understood. The small subunit scribed (ENFAet al. 1984; SAMBROOKet al. 1989). Following ribosomal RNA (SSU rRNA) sequence of both gene hybridization with oligonucleotides, filters were washed se- types has been determined only in Plasmodium berghei quentially in 2 X SSC/O.2% SDS for 30 min, 1 X SSC/O.2% SDS for 30 min and in 0.2 X SSC/O.2% SDS for 5 min at room and Plasmodium falciparum (GUNDERSONet al. 1987; temperature. Total RNA was extracted as previously described MCCUTCHANet al. 1988). Analysis of these sequences (ENFA et al. 1984) and electrophoresed on agarose/ showed that: (1)the isoforms of a given species are more formaldehyde gels (SAMBROOKet al. 1989). similar to one another thanto the corresponding genes A P. cynomolgi (Ceylon strain) EcoRI genomic library was in the otherspecies and (2) the pattern of conservation constructed by standard methods in lambda EMBL4 (SA" BROOK et al. 1989) and screened with oligonucleotides span- and divergence between isoforms is quite different in ning various conservedor semiconserved regions of the SSU the two species (ENEAand CORREDOR1991). To extend rRNA segment. the database we chose to analyze a lineage, Plasmodium Oligonucleotides: Oligonucleotides weresynthesized by cynomolgi, which being relatively distantly related to standard phosphoramidite chemistry in an Applied Biosystems P. berghei and to P. falciparum (DIGIOVANNI et al. 1990), 381A DNAsynthesizer according to the manufacturers instruc- tions and purified on acrylamide gels. Oligonucleotides were and being composed of a numberof closelyrelated but 5' end labeled with [y3*P]ATP(5000 Ci/mmol) (Amersham) distinct strains (COATNEYet al. 1971; GALINSKIet al. and polynucleotide kinase (SAMBROOKet aL 1989). Oligonuclee 1987),would allow for analyses involvinglarge phyletic tides representing sequenceswith different degreeof consem- distances as well as short ones. tion were used for screening genomic libraries, clone character- ization and sequencing in both orientations (Table 1). 'Present address: Department of Cell Biology and Anatomy, Box 1007, PCR amplification and cloning: The SSU rDNA comple- Mount Sinai Medical Center, OneGustave L. Levy Place, New York, NY 10029- ment was amplified using 1-2 ng of total parasite genomic 6574. DNA and oligonucleotides 0009 and 2134.1. Polymerasechain Genetics 136: 857-865 (March, 1994) 858 V. Corredor and V. Enea TABLE 1 Southern transfer was performed by standard methods ex- cept that the gel was soaked in 0.15 N HCl for 15 min before Catalog of oligonucleotides used in this study denaturation and neutralization. ~~ ~~~~ Primer Sequence" Primer Polarity Position Sequence determination and analysis: he sequence of the V7 variable region (NEEFSet al. 1990) of Berok and Ceylon 0001 actggcaagatcaaccaggtt 1-21 blood stage SSU rRNA was determined by primer extension of 0009 tgatcttgccagtagtcatat 9 -29 the RNA (GHOSH1980), using a 5' end-labeled oligonucleotide 0096 tataactgttttaatgagccg 97-1 17 (1880) and M-MLV reverse transcriptase. (BRL)and the product 0382 caggctccctctccggaatcg 389-409 sequenced by the "GILBERT (1980) method as described. 0401 atgacgggtaacggggaatta 363-383 0416 ctgctgccttccttagatgt 423-443 Sequence of genomic clones, pUC subclones and cloned 0508 gagaggtagtgacaagaaata 470 - 490 amplification products was determined by the chain termina- 0651 taactgcaacaattttaatat 611-631 tion dideoxy method (SANGERet al. 1977) using either Taq 1037 acctctgacatctgaatacaa 969-989 polymerase and 5' end-labeled primers as recommended by 1126 gaaagttaagggagtgaagac 1064-1084 the manufacturer (Promega) or T7 polymerase sequencing in 621 1 gtcttcactcccttaactttc 1064-1084 the presence of [35S]dCTPusing 7deaza dGTP and 7-deaza 1162 atagtttatggttaagatta 1099-1119 1171 aagactttgatttctcataag 1202-1222 dITP analogs as substitutes for dGTP as recommended by the 1273 gcttaatttgactcaacacgg 1303-1323 manufacturer (Pharmacia LKB). Oligonucleotides repre- 1508 tggttaattccgataacgaac 1432-1452 sented in Table 1 were used for sequencing reactions in both 1623 tctaacacaaggaagtttaag 1715-1735 orientations. 1880 gacatcacagacctgttgttg 1737-1757 Sequence alignments were generated with the PILEUP/ 1863 cttactaggcattcctcgttc 1944-1964 high road program available inthe Genetic Computer Group 1900 tgccctttgtacacaccgccc 1997-2017 2134 tcacctacggaaaccttgtta 2133-2153 (DEVEREUXet al. 1984) with gap opening penalties ranging 2134.1 gtaacaaggtttccgtaggtga 2133-2153 from 0.5 to 5.0 and gap extension penalties equal to 0.06 times 2401 aagaaagagctattaatctgt 1358-1378 or 0.1 times the gap opening penalty. Different alignments 2401.1 tgacagattaatagctctttc 1356-1376 were collated with the help of the program MALIGNED 3553 ctttaaaagataggatttacg 2099-2119 (CWUC 1991). V7-84 ccgcaatttagcaaatcatac 1562-1582 Parsimony and distance trees were generated with the PHY- V2-114 185-205 caagacattcgtcgagagtaa LIP package, version 3.4 (FELSENSTEIN1985). Bootstrappings v7-174 acaatcaascaaatacattct 1668-1688 were performed with the PHYLIP programs DNABOOT for * Primers display different degrees of conservation varying from parsimony-based trees, and SEQBOOT for distance-based universally conserved to strain specific. Primers used for screeningof trees (FELSENSTEIN1985). Alternative topologies were gener- libraries only hybridize to Plasmodium DNA as determined hy align- ing the corresponding region to the available SSU rRNA genes con- ated with the program DNAMOVE for parsimony (FELSENSTEIN sidered as more closely relatedto the parasite's hosts. These se- 1985) and TREE for distance trees (OLSEN1991). quences include the human (Homo sapiens), and mosquito (Aedes The three pairs of A and B sequences were analyzed in de- albopictus) genes. All primer sequences are oriented5' to 3' and the tail. A number of alignments were generated, including pair- polarity given with respect to the rRNA sequence. wise alignments, separate alignments of the A and B se- Primer positions refer to the sequence alignment in Figure 1. quences, alignments of the A and B sequences of pairs of reactions (PCR) were performed in 50-p1 reaction volumes, as species and alignments of the full set. For each of the above, described by SAIKI(1990) with the following parameters: 90"/1 several gap opening and gap extension penalties were used. min denaturation, 55"/1 min annealing and 71"/2 min poly- These were used as inputs to generate trees by distance meth- merization, for 30 cycles. ods, parsimony and evolutionary parsimony. PCR products were extracted with phenol,phenol/ chloroform and chloroform followed by two to threecycles of RESULTS ethanol precipitation in the presence of ammonium acetate. DNA was blunt ended with Klenow fragment (New England