Phylogenetic Relationship of Locally Isolated Paramecium Species Inferred from Histone H4 Genes
Phylogenetic Relationship of Locally Isolated Paramecium Species Inferred from Histone H4 Genes
Fareeda Tasneem1, Farah Rauf Shakoori1 and Abdul Rauf Shakoori2,*
1Department of Zoology, University of the Punjab, Quaid-i-Azam Campus, Lahore, 54590, Pakistan
2School of Biological Sciences, University of the Punjab, Quaid-i-Azam Campus, Lahore, 54590, Pakistan
The histone H4 gene is supposed to be highly conserved among all organisms with the exception of variations found in the ciliate species. Thus a fragment of histone H4 gene (160bp) was sequenced in ten locally isolated strains of Paramecium species including a standing-alone on the basis of 18SrDNA gene sequence strain FT8. In order to readdress the relationships of FT8 strain with other strains of Paramecium species, a molecular phylogenetic analysis was performed on the basis of H4 gene sequences. A phylogenetic tree was constructed using the aligned sequences of our study with 22 entries of closely related species from GenBank Data. Highly polymorphic sites were observed in FT8 strain sequence as compared to other species. Analysis of H4 gene sequences in our study showed they are closely related and behaved as a good substitute for phylogenetic analysis. Phylogenetic tree constructed by maximum likelihood method placed the FT8 strain very distinctly from other species, giving it a clearly separate position. Phylogeny results based on the H4 gene sequences corresponded to the results obtained by 18SrDNA sequences. Both markers proved the FT8 strain is a highly divergent species based on the phylogenetic relationships, opening new realms for the researchers of this field.
Received 16 February 2017
Revised 28 April 2017
Accepted 03 July 2017
Available online 08 September 2017
FRS and ARS conceived and supervised the study. FT executed the experimental work and wrote the manuscirpt.
Phylogenetic relationships of ciliates, Ribotyping, 18SrDNA analysis, FT8, Paramecium.
* Corresponding author: firstname.lastname@example.org
0030-9923/2017/0005-1767 $ 9.00/0
Copyright 2017 Zoological Society of Pakistan
Ciliates are characterized as monophyletic group of the unicellular eukaryote. Their distinctive feature is the presence of cilia and other infraciliature structures including kinetosomes (; ; ; ; ). All species of ciliates harbor two different kinds of nuclei; micronucleus that is required for sexual reproduction and other one is transcriptionally active macronucleus (; ; ; ; ;; ; ; ). Among them genus Paramecium of the oligohymenophorea ciliates is well studied in terms of its diversity, ultrastructure, physiology and genetics. Especially, the sibling species of the Paramecium aurelia complex have been intensively studied because of their complex genomics and as an example of species radiation (; ; ; ; ). It, therefore, becomes imperative to study the evolution and phylogenetic relationships of these popularly used eukaryotic model organisms. However, the systematics and phylogeny of Paramecium has been well developed in last decade (; ; ; ; ) but at the same time some relations within the genus remain dubious and questionable leaving many described species as invalid (; ; ; ; ; ).
Like many other molecular markers (rDNA, HSP70, COI, COII), genes for the histone proteins (mainly H3 and H4) are being used as phylogenetic markers because of their considerable variation found in ciliates (). It would be pertinent to mention that because of the highly conserved nature of these proteins in vast majority of organisms, histones are not considered to be appropriate markers for phylogenetic studies (; ; ). These proteins are very conservative also in ciliates; however, there were several attempts to use the fragments of their genes as phylogenetic markers to infer the relationships of closely related groups–at species or even within-species level. Recently, investigations on H4 proteins for understanding of the evolutionary relationship of Paramecium species have increased tremendously. For example, investigated the H4 gene of a number of species of the P. aurelia complex (P. aurelia: P. primaurelia, P. biaurelia, P. tetraurelia, P. pentaurelia, P. septaurelia, P. octaurelia, P. decaurelia, P. undecurelia, P. dodecaurelia, P. tredecaurelia and P. quadecaurelia) in order to untangle their phylogenetic connections. Similarly, H4 gene sequence have been studied for the phylogentic relationship of P. jeningsi strains by .
Independent evolution of histone proteins and their prominent sequence variation in ciliates as compared to other animals have made them potential candidate for cross check analysis. made a comparative analysis of H4 and H3 genes including the amino acid sequences of corresponding proteins, and an intergenic DNA fragment between them. According to them, phylogenetic relationships of P. tetraurelia and P. bursaria based on H4 and H3 sequences were similar to the rDNA phylogeny, thus making it valid for the studies on evolutionary relationships. compared thirteen ciliates including Paramecium with other eukaryotes. They found large variations among ciliates as compared to homologous fragments of other animals based on parsimony and maximum likelihood methods. This gives an idea of functional constraints of these proteins and their adaptive evolution probably due to nuclear dualism and peculiar organization of chromatin in macronucleus of ciliates. So the aim of the present study was to resolve the phylogenetic relationship of Paramecium species on the basis of analysis of H4 gene sequences. FT8 strain (Paramecium caudatum pakistanicus) being a unique species on the basis of 18SrDNA sequences already published by was also included in our study. This strain has been included in this study to make evolutionary comparison of this species on the basis of 18S rDNA and H4 gene sequences and to find its position in phylogenetic tree with reference to other Paramecium species.
Materials and Methods
Sampling and maintenance of Paramecium strains
Nine isolates of Paramecium sp. were collected from different regions of Punjab Province, Pakistan: two strains (FT2.1 and F3.1) from Kasur, three (FT4.1, FT10.1 and FT11.1) from Lahore, two (FT5.1, FT6.1) from Sheikhupura and three (FT7.1, FT9.1) from Mansehra. All of them were purified and acclimatized according to .
Genomic DNA extraction
Genomic DNA extraction was followed by incubation of Paramecium cells with 10mM Tris-HCl (pH 7.5) for 10-12h at 27°C. Incubated cells were harvested at low speed centrifugation of 6741x g at 7°C for 10 min. The cells were lysed in the lysis buffer (42% urea, 0.30M NaCl, 10Mm Tris-HCl pH 7.5, 10mM EDTA, 10% Nonidet P40 and 1% SDS). The lysis mixture was extracted two times with phenol:chloroform (1:1) following .
Amplification of H4 genes
A 160bp fragment of the H4 gene was ampliﬁed by using the primers H4-F02 (5’GGT ATT ACT AAG CCC GCT ATC AGA AGA3’) and H4-R02 (5’GGT CTT TCT TCT GGC GTG TTC AGT GTA3’) as used by , who applied H4 histone gene for phylogenetic analysis of P. jenningsi. PCR was performed in a final volume of 50µl containing 10mM dNTPs, 5µl MgCl, 3µl of 10x buffer (Thermo Fisher Scientific), 1µl 10mM of each primer and 1µl of 2.5 U Taq-polymerase (Thermo Fisher Scientific). PCR amplification protocol composed of 1 cycle at 94°C for 5 min, followed by 35 cycles, each of 2 min denaturation at 94°C, 1 min annealing at 54°C, 2 min extension at 72°C and final extension at 72°C for 20 min. Amplification was performed in Gene Amp PCR System 9700 (Applied Biosystem).
Sequencing of amplified products
PCR products, after their appropriate size was confirmed on 0.8% agarose gel, were cut and purified using Nucleospin Extract 11 (Macherey-Nagel Germany). Purified products were sequenced by Macrogen (South Korea).
The nucleotide sequences of all analyzed strains were compared with the sequences present in database using online website program (). All the sequences were aligned by online Muscle alignment program () with the available sequences from GenBank/EMBL databases under the following accession numbers: P. tetraurelia (XM001425872), P. tetraurelia (XM001455606), P. tetraurelia (XM001452606), P. tetraurelia (XM001452073), P. tetraurelia (XM001459068). P. decaurelia (DQ067622), P.tetraurelia (AJ004699), P. tetraurelia (XM001442554), 19 P. tredecaurelia (DQ067629), P. pentaurelia (DQ067623), P. primaurelia (DQ067620), P. quadecuarelia (DQ0676630), P. decaurelia (DQ067626). P. jenningsi (DQ001056), P. jenningsi (DQ001064), P. jenningsi (DQ001062), P. jenningsi (DQ001059), P. jenningsi (DQ001061), P. jenningsi (DQ001057). P. undecaurelia (DQ067627), P. septaurelia (DQ067624), P. tertraurelia (AJ004700), P. caudatum (AB670962) and P. bursaria (AJ004702). Percentage and distance matrix of Paramecium species of present study with the above mentioned Paramecium species from GenBank were compared. The phylogenetic tree was constructed to get the final position of all isolated Paramecium species along with FT8 strain by using online website program () by . Maximum likelihood (neighbor joining) method was used to construct phylogenetic trees describing the relationships of the examined strains. Accuracy of inferred topologies was assessed via bootstrap analysis of 1000 replicates.
PCR amplification reactions of all nine strains yielded a fragment of approximately 160bp in size. Sequences of the histone H4 gene fragment of all strains were obtained. Sequencing was based on both the forward and reverse primers used for amplification. Sequences were verified for identity by comparing with the H4 gene sequences of Paramecium species already available in GenBank.
Careful alignment of the nucleotide sequences derived from the H4 genes of the ten analyzed strains of Paramecium species revealed sixteen polymorphic sites. These sites are indicating both the transition (A to G, C to T) and transversion (A to T, T to G) substitutions, out of which four transitions revealed by only FT8 strain. All of these strain sequences were also aligned to homologous H4 gene sequences of Paramecium species obtained from GenBank. This comparison with other species displayed genetic diversity with 25 polymorphic sites involving both types of mutations (). The extents of divergence between the studied isolates are displayed in the distance
Phylogeny constructed by bootstrap analysis revealed a similar pattern of species emplacement within the tree. Most of the species collected from GenBank as a result of blast analysis of our strain sequences were belonging to the species of the P. aurelia complex. The tree represents nine entries of Paramecium tetraurelia under variety of accession numbers submitted by different authors along with the sequence of fifteen other species. All of the species in phylogenetic tree are not making any definite group or clades, rather are dispersed throughout the tree. However, most of the strains under our study did not show close relationship with any of the species. Rather all of them fell into same group by making their own clade. Whereas, FT8 was the only strain found to be at very distinct position from the rests of the species.
With the early studies on protein it was found out that histone H4 was strongly conserved throughout the plants and animal kingdom. Despite the fact that H4 protein sequences are invariant among vertebrates, considerable differences have been observed within ciliates allowing them to systematize their position in phylogenetic tree. H4 sequence variation was first time described in Tetrahymena pyriformis by . In spite of the fact that variations in the H4 gene sequences of different species of ciliates do not show the direct effect of the reproductive affinity, they can contribute the differences at genetic level and play an alternative role in the evolution of whole genome duplication that may lead to complete reproductive isolation .
Recently, the trend is shifting to confirm the phylogenetic results based on mitochondrial and rRNA genes with histone H4 sequences. Likewise, many things left undecided by mtDNA and rRNA genome sequences are needed to be conﬁrmed with some other data. Histone H4 gene has been selected as a reliable marker for the testing of phylogenies based on rRNA or mtDNA sequence. Mostly it has been observed that variations in the nucleotide sequences often occur as synonymous substitutions that tend to arrange the remote species close to each other in the phylogenetic tree (). So in order to confirm the exact position of FT8 strain () in phylogenetic tree, H4 gene sequence was selected as a better substitute along with some other strains. This has also been reported that species of Paramecium aurelia complex contains a single gene of H4 (). This single gene variations have been extensively studied which is the source of genetic isolation of the P. aurelia sibling species, exists from one or more than one site changes “but constitute a dense evolutionary cluster” .
Species of our study mostly fall under the complex of Paramecium aurelia species except FT8 strain that has already been reported as a new species based on the 18SrRNA phylogeny. Present results based on Histone H4 gene sequences confirmed the FT8 strain has a peculiar phylogenetic position in Paramecium. This species is showing its distinct position representing it as a highly divergent from other species. Rests of the strains belonging to Paramecium aurelia species are showing their positions a little different from each other. The reason could be the isolation of strains from different regions that caused them to make the geographical groups. All of the strains are falling together in same position making their own clade.
According to the tree, FT8 as well as P. caudatum sequence, are interspersed in the mass of the P. aurelia entries that is revealing not a very good picture. This is the weak point of this study which has a risk to create biases in the phylogeny. Although, it can be true with this set of data, as the fragment of the sequence analyzed was very short. However, this drawback need to be resolved by amplifying the bigger portion of H4 gene. The overall picture of the strains under present study in the tree is representing their close relationship with P. aurelia species except FT8 strain which is showing its completely different position as compared to others, indicating it as a new species.
Statement of conflict of interest
Authors have declared no conflict of interest.
Aury, J.M., Jaillon, O., Duret, L., Noel, B., Jubin, C., Porcel, B.M., Segurens, B., Daubin, V., Anthouard, V., Aiach, N., Arnaiz, O., Billaut, A., Beisson, J., Blanc, I., Bouhouche, K., Camara, K., Duharcourt , S., Guigo, R., Gogendeau, D., Katinka, M., Keller, A.M., Kissmehl, R., Klotz, C., Koll, F., Mouel , A.L., Lepere, G., Malinsky, S., Nowacki, M., Nowak, J.K., Plattner, H., Poulain, H., Ruiz, F., Serrano, F., Zagulski, M., Dessen, P., Betermier, M., Weissenbach, J., Scarpelli, C., Schachter, V., Sperling, L., Meyer, E., Cohen, J and Wincker, P., 2006. Global trends of whole-genome duplications revealed by the ciliate Paramecium tetraurelia. Nature, 444: 171-178.
Baroin-Tourancheau, A., Delgado, P., Perasso, R. and Adoutte, A., 1992. A broad molecular phylogeny of ciliates: Identification of major evolutionary trends and radiations within the phylum. Proc. natl. Acad. Sci. USA., 89: 9764-9768.
Barth, D., Krenek, S., Fokin, S.I., and Berendonk, T.U., 2006. Intraspecific genetic variation in Paramecium revealed by mitochondrial cytochrome c oxidase I sequences. J. Eukary. Microbiol., 53: 20-25.
Boscaro, V., Fokin, S., Verni, F and Petroni, G., 2012. Survey of Paramecium duboscqui using three markers and assessment of the molecular variability in the genus Paramecium. Mol. Phylogen. Evolut., 65: 1004–1013.
Dereeper, A., Guignon, V., Blanc, G., Audic, S., Buffet, S., Chevenet, F., Dufayard, J.F., Guindon, S., Lefort, V., Lescot, M., Claverie, J.M. and Gascuel, O., 2008. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucl. Acids Res., 36(Web Server issue): W465-9.
Fokin, S.I., Przyboś, E., Chivilev, S.M., Beier, C.L,, Horn, M., Skotarczak, B., Wodecka, B. and Fujishima, M., 2004. Morphological and molecular investigations of Paramecium schewiakoffi sp. nov. (Ciliophora, Oligohymenophorea) and current status of distribution and taxonomy of Paramecium spp. Eur. J. Protistol., 40:225–243.
Gould, S.B., Kraft, L.G.K., Dooren, G.G.V., Goodman, C.D., Ford, K.L., Cassin, A.M., Bacic, A., McFadden, G., and Waller, R.F., 2011. Ciliate pellicu lar proteome identiﬁes novel protein families with characteristic repeat motifs that are common to alveolates.. Mol. Biol. Evolut., 28: 1319-1331.
Hampl, V., Hug, L., Leigh, J.W., Dacks, J.B., Lang, B.F., Simpson, A.G. and Roger, A.J., 2009. Phylogenomic analyses support the monophyly of Excavata and resolve relationships among eukaryotic ‘‘supergroups’’. Proc. natl. Acad. Sci. U.S.A., 106: 3859-3864.
Katz, L.A., Snoeyenbos-West, O. and Doerder, F.P., 2006. Patterns of protein evolution in Tetrahymena thermophila: implications for estimates of effective population size. Mol. Biol. Evolut., 23: 608-614.
Lanzoni, O., Fokin1, S., Lebedeva, N., Migunova, A., Petroni1, G. and Potekhin, A., 2016. Rare freshwater ciliate Paramecium chlorelligerum Kahl, 1935 and its macronuclear symbiotic bacterium “Candidatus Holospora parva” Olivia Lanzoni. PLoS One, 11: e0167928.
McGrath, C.L., Gout, J.F., Johri, P., Doak, T.G. and Lynch, M., 2014. Differential retention and divergent resolution of duplicate genes following whole-genome duplication. Genome Res., 24: 1665-1675.
Pineau, P., Henry, M., Suspene, R., Marchio, A., Dettai, A., Debruyne, R., Petit, T., Lecu, K.L., Moisson, P., Dejean, A., Wain-Hobson, S. and Vartanian, J.P.A., 2005. Universal primer set for PCR amplification of nuclear Histone H4 genes from all animal species. Mol. Biol. Evolut., 22: 582-588.
Przybos, E., Greczek-Stachura, M., Prajer, M., Potekhin, A. and Cotsinian, A., 2008. Two species of the Paramecium aurelia complex (Ciliophora, Protista) from the Black Sea region (Russia) with their RAPD-PCR fingerprints characteristics. Protistology, 5: 207-212.
Przybos, E., Maciejewska, A. and Skotarczak, B., 2006. Relationships of species of the Paramecium aurelia complex (Protozoa, Ph. Ciliophora, Cl. Oligohymenophorea) based on sequences of the histone H4 gene fragment. Folia boil. (Kraków), 54: 1-2.
Shakoori, F.R., Tasneem, F., Al-Ghanim, K., Mahboob, S., Al-Misned, F., Jahan, N. and Shakoori, A.R., 2014. Variability in secondary structure of 18S Ribosomal RNA as topological marker for identification of Paramecium species. J. cell. Biochem., 115: 2077-2088.
Struder-Kypke, M.C., Wright, A.D.G., Fokin, S. and Lynn, D.H., 2000. Phylogenetic relationships of the genus Paramecium inferred from small subunit rRNA gene sequences. Mol. Phylogenet. Evolut., 14: 122-130.
Tarcz, S., Potekhin, A., Rautian, M. and Przybos, E., 2012. Variation in ribosomal and mitochondrial DNA sequences demonstrates the existence of intraspeciﬁc groups in Paramecium multimicronucleatum (Ciliophora, Oligohymenophorea). Mol. Phylogenet. Evolut., 63: 500-509.
Wang, L., Zhao, Y., Qian, C., Wei, G., Zhu, B. and Liu, C., 2016 Gene expressions of heat shock proteins in Bombyx mori egg parasitized by a parasitoid wasp, Telenomus theophilae. Pakistan J. Zool., 48: 507-512.