MBE Advance Access originally published online on November 24, 2005
Molecular Biology and Evolution 2006 23(3):608-614; doi:10.1093/molbev/msj067
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Article |
Patterns of Protein Evolution in Tetrahymena thermophila: Implications for Estimates of Effective Population Size


* Department of Biological Sciences, Smith College, Northampton;
Program in Organismic and Evolutionary Biology, University of Massachusetts Amherst; and
Department of Biological, Geological and Environmental Sciences, Cleveland State University
E-mail: lkatz{at}smith.edu.
| Abstract |
|---|
|
|
|---|
High levels of synonymous substitutions among alleles of the surface antigen SerH led to the hypothesis that Tetrahymena thermophila has a tremendously large effective population size, one that is greater than estimated for many prokaryotes (Lynch, M., and J. S. Conery. 2003. Science 302:14011404.). Here we show that SerH is unusual as there are substantially lower levels of synonymous variation at five additional loci (four nuclear and one mitochondrial) characterized from T. thermophila populations. Hence, the effective population size of T. thermophila, a model single-celled eukaryote, is lower and more consistent with estimates from other microbial eukaryotes. Moreover, reanalysis of SerH polymorphism data indicates that this protein evolves through a combination of vertical transmission of alleles and concerted evolution of repeat units within alleles. SerH may be under balancing selection due to a mechanism analogous to the maintenance of antigenic variation in vertebrate immune systems. Finally, the dual nature of ciliate genomes and particularly the amitotic divisions of processed macronuclear genomes may make it difficult to estimate accurately effective population size from synonymous polymorphisms. This is because selection and drift operate on processed chromosomes in macronuclei, where assortment of alleles, disruption of linkage groups, and recombination can alter the genetic landscape relative to more canonical eukaryotic genomes.
Key Words: genome evolution effective population size ciliates protein evolution Tetrahymena thermophila
| Introduction |
|---|
|
|
|---|
Characterizing DNA sequence variation in natural populations of Tetrahymena thermophila is essential for elucidating the evolutionary history of this model organism. Tetrahymena thermophila has been a rich resource for discovery, including self-splicing RNAs and telomeres. Yet, as with many model organisms, only a limited amount of data is available on the ecology and natural history of this model eukaryote. For example, the phylogeography of T. thermophila is relatively unexplored, though populations have been observed along the east coast of the United States (specifically, Massachusetts, Vermont, New Hampshire, Maine, Pennsylvania, Florida) and as far west as Illinois (Nanney and Simon 1999
The elevated numbers of synonymous substitutions among alleles of the surface antigen gene SerH (Gerber et al. 2002
) were interpreted as showing that T. thermophila has a very large effective population size (Lynch and Conery 2003
). In fact, in a study analyzing a diversity of organisms from vertebrates to prokaryotes, T. thermophila has the second highest estimate of effective population size (measured as Neµ)its value falls between estimates from several prokaryotic species with larger census population sizes and geographic distributions (i.e., Procholorococcus sp., Salmonella enterica, and Legionella penumophila) (Lynch and Conery 2003
). However, the estimate of Neµ in T. thermophila was made from only a single locus, SerH, and polymorphisms in any single gene may reflect evolutionary processes (i.e., balancing selection) operating on a single locus rather than population history. Furthermore, based on the estimate of Neµ derived from the SerH data, T. thermophila appears to be an exception to the genome complexity argument that states that genome sizes are largest in lineages with small effective population sizes due to the influence of genetic drift (Lynch and Conery 2003
). In contrast to observations for other eukaryotic taxa, T. thermophila has an estimated effective population size larger than most prokaryotes (Lynch and Conery 2003
) and also a relatively large macronuclear genome (
110 Mb; GenBank accession number AAGF01000000).
The SerH gene is one of a family of Ser genes in Tetrahymena that encodes abundant glycosylphosphatidylinositol-linked surface proteins known as immobilization antigens (Ko and Thompson 1992
; Smith et al. 1992
). As with the variant surface proteins of many parasitic protists (Kusch and Schmidt 2001
), the proteins encoded by these genes are modular, composed of tandem repeats containing even numbers of periodic cysteines. For instance, SerH genes encode H proteins containing 3.54.5 repeats each with eight cysteines, whereas the SerL genes encode L proteins containing 25 repeats each with six cysteines (Tondravi et al. 1990
; Doerder and Gerber 2000
). While the sequencing of SerH alleles from T. thermophila populations revealed considerable sequence variation, the overall structure, hydrophobicity and cysteine periodicity were conserved (Gerber et al. 2002
). At the nucleotide level, 9.7%26.7% of sites were polymorphic in pairwise comparisons, and a high rate of synonymous substitution was observed. Though the exact role of the SerH proteins is unknown, their ubiquitous presence on the surface of Tetrahymena and their tight environmental and genetic regulation suggest that they are ecologically important.
The SerH locus is also associated with elevated levels of recombination within macronuclei when compared to other loci (Deak and Doerder 1998
). In ciliates, macronuclear genomes are the site of all gene expression while the germline micronucleus is transcriptionally inactive and exchanged during conjugation (reviewed in Jahn and Klobutcher 2002
; Yao, Duharcourt, and Chalker 2002
; McGrath and Katz 2004
). For two SerH1 nonsense mutations separated by 726 nt, recombination forming wild-type and double-mutant alleles occurs in virtually 100% of heterozygous macronuclei, a recombination rate significantly higher than that observed for other macronuclear markers (Lovlie, Haller, and Orias 1988
; Longcor et al. 1996
). Whereas only micronuclear recombination can generate new hereditary genotypes, macronuclear recombination can yield new gene combinations within an asexual clonal lineage and therefore may be advantageous during asexual periods of the ciliates life cycle.
Our aim was to distinguish between two hypotheses on protein evolution and effective population size in T. thermophila. The first hypothesis is that the high level of polymorphism at SerH, including numerous synonymous substitutions among alleles, is indicative of a large effective population of this ciliate and is paralleled by high levels of polymorphisms at other loci (Gerber et al. 2002
). Under the second hypothesis, the elevated level of polymorphism in SerH is anomalous and reflects a nonneutral pattern of molecular evolution in this surface antigen gene. To distinguish between these hypotheses, we analyzed variation in four nuclear genes and one mitochondrial gene from a combination of inbred lines and natural isolates of T. thermophila, some of which contain previously characterized divergent SerH alleles (Gerber et al. 2002
).
| Materials and Methods |
|---|
|
|
|---|
Polymerase Chain Reaction, Cloning, and Sequencing
We characterized polymorphisms in three protein-coding and two small-subunit (SSU)-rDNA genes (mitochondrial and nuclear) to assess different regions of the T. thermophila genome. Polymorphisms are expected to accumulate at synonymous sites in protein-coding genes and in regions of low constraint (i.e., some loops) of the rDNAs. Four genes (SSU-rDNA, mt-SSU,
tub, Ef1
) were characterized in the Katz laboratory, where samples were sent "blind" to hide source (natural population vs. clonal line). We relied on previously published primers for amplification of
tub (Israel et al. 2002
(Baldauf 1999
1.6 kb of mt-SSU were designed in our laboratory: SSU mito 159+ (ATAGCGTMHDGGTGMGTAAWA) and SSU mito 1840 (CCACCRRCARKTTCCCMTACC). All products were amplified using Platinum Taq DNA polymerase (Invitrogen, Carlsbad, Calif., #10966026). A 1-kb fragment of the actin gene was characterized in the Doerder laboratory from five natural populations using primers ATGGCTGAAAGTGAATCCC and CTCAGGAGGAGCAACAAC and a high fidelity Z-Taq (Takara, Mirus Bio, Madison, Wisc.) Actin products, analyzed in the Doerder laboratory, were cloned into pGEM-Teasy for sequencing. Resulting sequences were similar to sequences on GenBank (table 1).
|
Sequence Analysis
The sequences obtained in this study were aligned using the multisequence alignment algorithms and ClustalW (Thompson, Higgins, and Gibson 1994
SerH Analysis
Analysis of the evolution of SerH combined published alleles (table 2) with homologous regions form the T. thermophila genome. SerH loci were identified in the completed T. thermophila macronculear genome by Blast analysis, using the Assembly_2 November 2003 scaffolds available at The Institute for Genomic Research http://tigrblast.tigr.org/er-blast. Identification of SerH ESTs was performed by both inspecting the results of Blast analysis of ESTs from several libraries (http://cbr-rbc.nrc-cnrc.gc.ca/reith/tetra/tetra.html) and by searching for similar ESTs to SerH alleles from natural populations using SeqMan (DNAStar, Inc.). SerH alleles and scaffold regions were aligned by ClustalW in Megalign (DNAStar, Inc.). Genealogical analyses of SerH amino acid repeats and of scaffold regions used both neighbor-joining (NJ) and maximum parsimony algorithms. Mean distances were estimated for NJ analysis. Parsimony analyses treated gaps as a 21st state and used 10 random addition sequences in heuristic searches. Bootstraps were calculated using 100 replicates under all models.
|
| Results |
|---|
|
|
|---|
Our data indicate that the degree of DNA polymorphism at the SerH locus (Gerber et al. 2002
tub), elongation factor 1
(Ef1
), actin, and small-subunit ribosomal RNA (SSU-rDNA), as well as the mitochondrial small-subunit rRNA (mt-SSU). Genomic DNAs from natural populations and clonal lines were chosen for this study because they are known to contain divergent SerH alleles (table 1; Gerber et al. 2002
tub, and Ef1
) and one mitochondrial gene (mt-SSU; table 3). Furthermore, there is no evidence of elevated diversity either within the natural population or between inbred lines. In fact, all four of these loci have roughly the same level of polymorphisms, and at least some of the low level of polymorphisms is likely due to either experimental (polymerase chain reaction/sequencing) error or are variants in the polyploid macronucleus. Similarly, a comparison of a 969-bp fragment of actin genes from five different natural isolates (four from Allegheny National Forest and one from New Hampshire) confirms the findings as we uncovered only five polymorphic sites and an average pairwise divergence of 0.27% (table 3). These data are in stark contrast to the high divergence, 20.69% nucleotide and 31.15% amino acid, observed at the SerH locus in the same populations (table 3; Gerber et al. 2002
|
To provide an explanation for the contrast between the high level of SerH polymorphism and lack of variation at other loci, we reanalyzed the published SerH allele sequences based on the patterns of substitutions in three of the 3.54.5 imperfect, cysteine-rich repeat motifs within each allele (table 2; Gerber et al. 2002
|
Genealogical analyses of the first three of the imperfect repeats in SerH reveal a complex pattern of evolution through both vertical descent and concerted evolution. Repeats within one group of alleles, family I (fig. 1c), are very similar to one another, consistent with concerted evolution driving the diversity of these alleles. Within allele family I, repeats in SerH5 are virtually identical at the amino acid level, and repeats within SerH6 are closely related (fig. 1c). Similarly, the first repeat of the SerH1 allele (SerH1_r1) is recently derived from the second repeat of the same allele (SerH1_r2). In contrast, allele families IIIV arose from independent concerted events followed by a period of vertical transmissionrepeats within these families are more similar between alleles than within alleles (fig. 1c). For example, repeats within allele family IV (alleles H7, H9, and H12) have diverged due to vertical transmission following concerted evolution of alleles (fig. 1c). Similarly, allele family II (H8, H10, and H11) also diverged following what appears to be a more ancient concerted eventrepeats within this family do not form a monophyletic group (fig. 1c).
Comparisons of coding and presumed untranslated regions (UTRs) of SerH alleles reveal a pattern of evolution consistent with the mixed ancestry of the SerH alleles as there is more polymorphism in coding domains compared across alleles than in both 5' and 3' UTRs (table 4). For example, among published alleles 5' and 3' UTRs diverge by 12.7% and 6.0%, respectively, while nucleotide divergence in coding domains in the same alleles is 20.2% (table 2). This pattern is consistent with a recombination-based mechanisms driving evolution of repeat regions within SerH alleles and implies that the concerted evolution is regulated so as not to involve UTRs.
|
Comparisons between SerH alleles characterized from several lines with the single macronuclear scaffold available from one inbred strain indicate that there is even a greater diversity of SerH in the macronuclear genomes of natural populations of T. thermophila. Blast analysis of SerH alleles against the completed T. thermophila macronuclear genome identified five loci related to SerH on a single macronuclear scaffold (fig. 2). Only one of these loci, Scaf_5, is identical in sequence to a published allele (SerH3). This macronuclear locus appears to be the only expressed site in the macronucleus as it is the only match found in comparisons with available T. thermophila EST data. Of the remaining four scaffold regions, one (Scaf_1) shows very low levels of similarity to other macronuclear scaffold regions and to published alleles. The other scaffold regions have complex relationships to published alleles, as revealed by genealogical analysis (fig. 3): intriguingly, the five divergent loci on the sequenced macronuclear chromosomes of a laboratory strain of T. thermophila fall in only two of four classes of alleles. This observation suggests that additional variants of SerH will be found on the macronuclear scaffolds of other strains of T. thermophila.
|
|
| Discussion |
|---|
|
|
|---|
Contrary to previous analysis (Lynch and Conery 2003
, SSU-rRNA,
tub, actin, and mt-SSU genes indicate that T. thermophila has a smaller effective population size and that the elevated levels at SerH are aberrant. Based on synonymous differences among alleles in our sampling of three protein-coding genes, we used equations from Lynch and Conery (2003)
0.0008 (table 3). This estimate of Neµ is less than 100 times smaller than the estimate of
0.18 that was calculated based on polymorphisms in the SerH gene. A smaller effective population size for T. thermophila is consistent with the hypothesis of Lynch and Conery (2003)
110 Mb) and large number of genes (
25,000). While a more accurate estimate of Neµ for this ciliate will require a larger, random sample of alleles, the estimate of Neµ in this study (table 3) is in line with estimates for other microbial eukaryotes including Plasmodium falciparum (Neµ = 0.00082) and Encephalitozoon cuniculi (Neµ = 0.003). Both of these eukaryotic microbes are parasiticthey can achieve large census population sizes and have relatively broad geographical distributions (generally tropical and subtropical for P. falciparum; recently commensal with humans for E. cuniculi). Similar estimates of T. thermophila census sizes and current distributions are not available, making it difficult to interpret estimates of effective population size for this model ciliate.
The elevated levels of polymorphism at SerH (table 3) coupled with the mixed patterns of vertical descent and gene conversion of repeats (fig. 1c) indicate that the SerH-coding regions are not evolving neutrally. Instead, the Ser H locus may be evolving under adaptive evolution (e.g., balancing selection) to maintain variation in this surface antigen in populations of T. thermophila. Such a pattern of evolution is analogous to that found in variable surface antigens in Trypanosomes (Donelson 2003
) and the V(D)J system of human immune systems (Lewis and Wu 1997
; Diaz and M.F. 1998
), where balancing selection maintains polymorphisms at loci involved in host-parasite interactions. While the function of the SerH-encoded proteins is unknown, the Ser gene families are regulated by environmental factors, and the resulting proteins are expressed on the surface of the ciliate (Margolin, Loefer, and R.D. 1959
; Allen and Gibson 1973
). Many more candidate Ser gene families with the same type of modular structure are present in the newly sequenced T. thermophila genome (F. P. Doerder, unpublished data), indicating that these antigen genes may constitute a substantial portion of the Tetrahymena genome.
Finally, the rapid amino acid evolution at the SerH locus is consistent with our hypothesis (Katz et al. 2004
) that differential selection on the functional macronucleus and the transcriptionally inactive micronucleus contributes to elevated rates observed in many ciliate proteins (Bhattacharya and Weber 1997
; Moreira, Le Guyader, and Philippe 1999
; Katz et al. 2004
). Ciliates, like animals and a few other eukaryotic lineages, have distinct germline and somatic genomes throughout their life cycles (reviewed in Zufall, Robinson, and Katz 2005
). However, only in ciliate macronuclei can differential amplification of chromosomes, recombination, and breakup of linkage groups create variation among macronuclei during asexual divisions. All of these effects are enhanced during asexual divisions by amitosis, the process by which macronuclei divide without accurately distributing sister chromatids between daughter nuclei. We argue that genome processing combined with amitosis can mask the effects of deleterious mutations in ciliate macronuclei (Katz et al. 2004
). As a result, proteins such as SerH can explore sequence space during iterative asexual cycles by effectively hiding alleles from selection in macronuclei while accumulating potentially compensatory mutations in the micronucleus. Hence, it may be difficult to infer the effective population size of ciliates because selection on processed macronuclear genomes can affect the rates of fixation of synonymous substitutions when compared to other eukaryotes with more canonical genomes.
| Acknowledgements |
|---|
|
|
|---|
This work was support by the National Science Foundation grants DEB-0092908 and DEB-0079325 to L.A.K. and funds provided to F.P.D. by the Graduate College and Department of Biological, Geological and Environmental Sciences at Cleveland State University. Thanks to Elyse Lasser for her work in helping to characterize T. thermophila sequences.
| Footnotes |
|---|
1 Present address: Department of Plant, Soil, and Insect Sciences, Fernald Hall, University of Massachusetts at Amherst, Amherst, MA 01003
| References |
|---|
|
|
|---|
Allen, S. L., and I. Gibson. 1973. Genetics of Tetrahymena. Pp. 307373 in A. M. Elliott, ed. Biology of Tetrahymena. Dowden, Hutchinson and Ross, Stroudsburg, Pa.
Baldauf, S. L. 1999. A search for the origins of animals and fungi: comparing and combining molecular data. Am. Nat. 154:S178S188.[CrossRef][Medline]
Bhattacharya, D., and K. Weber. 1997. The actin gene of the glaucocystophyte Cyanophora paradoxa: analysis of the coding region and introns, and an actin phylogeny of eukaryotes. Curr. Genet. 31:439446.[CrossRef][Web of Science][Medline]
Deak, J. C., and F. P. Doerder. 1998. High frequency intragenic recombination during macronuclear development in Tetrahymena thermophila restores the wild-type SerH1 gene. Genetics 148:11091115.
Diaz, M., and F. M. Flajnik. 1998. Evolution of somatic hypermutation and gene conversion in adaptive immunity. Immunol. Rev. 162:1324.[CrossRef][Web of Science][Medline]
Doerder, F. P., and C. A. Gerber. 2000. Molecular characterization of the SerL paralogs of Tetrahymena thermophila. Biochem. Biophys. Res. Commun. 278:621626.[CrossRef][Web of Science][Medline]
Donelson, J. E. 2003. Antigenic variation and the African trypanosome genome. Acta Trop. 85:391404.[CrossRef][Web of Science][Medline]
Gerber, C. A., A. B. Lopez, S. J. Shook, and F. P. Doerder. 2002. Polymorphism and selection at the SerH immobilization antigen locus in natural populations of Tetrahymena thermophila. Genetics 160:14691479.
Israel, R. L., S. L. Kosakovsky Pond, S. V. Muse, and L. A. Katz. 2002. Evolution of duplicated alpha-tubulin genes in ciliates. Evolution 56:11101122.[CrossRef][Web of Science][Medline]
Jahn, C. L., and L. A. Klobutcher. 2002. Genome remodeling in ciliated protozoa. Ann. Rev. Microbiol. 56:489520.[CrossRef][Web of Science][Medline]
Katz, L. A., E. Lasek-Nesselquist, J. Bornstein, and S. V. Muse. 2004. Dramatic diversity of ciliate histone H4 genes revealed by comparisons of patterns of substitutions and paralog divergences among eukaryotes. Mol. Biol. Evol. 21:555562.
Ko, Y. G., and G. A. Thompson. 1992. Immobilization antigens from Tetrahymena thermophila are glycosyl-phosphatidylinositol-linked proteins. J. Protozool. 39:719723.[Medline]
Kumar, S., K. Tamura, and M. Nei. 2004. MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief. Bioinform. 5:150163.
Kusch, J., and H. J. Schmidt. 2001. Genetically controlled expression of surface variant antigens in free-living protozoa. J. Membr. Biol. 180:101109.[CrossRef][Web of Science][Medline]
Lewis, S. M., and G. E. Wu. 1997. The origins of V(D)J recombination. Cell 88:159162.[CrossRef][Web of Science][Medline]
Longcor, M. A., S. A. Wickert, M.-F. Chau, and E. Orias. 1996. Coassortment of genetic loci during macronuclear division in Tetrahymena thermophila. Eur. J. Protistol. 32(Suppl.):8589.
Lovlie, A., B. L. Haller, and E. Orias. 1988. Molecular evidence for somatic recombination in the ribosomal DNA of Tetrahymena thermophila. Proc. Natl. Acad. Sci. USA 85:51565160.
Lynch, M., and J. S. Conery. 2003. The origins of genome complexity. Science 302:14011404.
Margolin, P., J. B. Loefer, and R. D. Owen. 1959. Immobilizing antigens of Tetrahymena pyriformis. J. Protozool. 6:207215.
McGrath, C. L., and L. A. Katz. 2004. Genome diversity in microbial eukaryotes. Trends Ecol. Evol. 19:3238.
Medlin, L., H. J. Elwood, S. Stickel, and M. L. Sogin. 1988. The characterization of enzymatically amplified eukaryotes 16S-like ribosomal RNA coding regions. Gene 71:491500.[CrossRef][Web of Science][Medline]
Moreira, D., H. Le Guyader, and H. Philippe. 1999. Unusually high evolutionary rate of the elongation factor 1
genes from the ciliorphora and its impact on the phylogeny of eukaryotes. Mol. Biol. Evol. 16:234245.[Abstract]
Nanney, D. L., and E. M. Simon. 1999. Chapter 1. Laboratory and evolutionary history of Tetrahymena. Pp. 125 in D. Asai and J. Forney, eds. Methods in cell biology 62, Tetrahymena. Academic Press, New York.
Smith, D. L., M. S. Berkowitz, D. Potoczak, M. Krause, C. Raab, F. Quinn, and F. P. Doerder. 1992. Characterization of the T-cell, L-cell, I-cell, S-cell, M-cell and P-cell surface (immobilization) antigens of Tetrahymena-Thermophilamolecular-weights, isoforms, and cross-reactivity of antisera. J. Protozool. 39:420428.[Medline]
Swofford, D. 2002. PAUP*: phylogenetic analysis using parsimony (*and other methods). Sinauer Associates, Sunderland, Mass.
Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. Clustal-Wimproving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:46734680.
Tondravi, M. M., R. L. Willis, H. D. Love, and G. A. Bannon. 1990. Molecular characterization of SerH3, a Tetrahymena thermophila gene encoding a temperature regulated surface antigen. Mol. Cell. Biol. 10:60916096.
Yao, M. C., S. Duharcourt, and D. L. Chalker. 2002. Genome-wide rearrangements of DNA in ciliates. Pp. 730758 in N. L. Craig, R. Craigie, M. Gellert, and A. Lambowitz, eds. Mobile DNA II. ASM Press, Washington D.C.
Zufall, R. A., T. Robinson, and L. A. Katz. 2005. Evolution of developmentally regulated genome rearrangements in eukaryotes. J. Exp. Zool. B Mol. Dev. Evol. 304B:448455.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
C.-H. Kuo, N. A. Moran, and H. Ochman The consequences of genetic drift for bacterial genome complexity Genome Res., August 1, 2009; 19(8): 1450 - 1454. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Gentekaki and D. H. Lynn High-Level Genetic Diversity but No Population Structure Inferred from Nuclear and Mitochondrial Markers of the Peritrichous Ciliate Carchesium polypinum in the Grand River Basin (North America) Appl. Envir. Microbiol., May 15, 2009; 75(10): 3187 - 3195. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Catania, F. Wurmser, A. A. Potekhin, E. Przybos, and M. Lynch Genetic Diversity in the Paramecium aurelia Species Complex Mol. Biol. Evol., February 1, 2009; 26(2): 421 - 431. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Snoke, T. U. Berendonk, D. Barth, and M. Lynch Large Global Effective Population Sizes in Paramecium Mol. Biol. Evol., December 1, 2006; 23(12): 2474 - 2479. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Zufall, C. L. McGrath, S. V. Muse, and L. A. Katz Genome Architecture Drives Protein Evolution in Ciliates Mol. Biol. Evol., September 1, 2006; 23(9): 1681 - 1687. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




