The early adaptive evolution of calmodulin

Interaction between gene duplication and natural selection in molecular evolution was investigated utilizing a phylogenetic tree constructed by the parsimony procedure from amino acid sequences of 50 calmodulin- family protein members. The 50 sequences, belonging to seven protein lineages related by gene duplication (calmodulin itself, troponin-C, alkali and regulatory light chains of myosin, parvalbumin, intestinal calcium-binding protein, and glial S-100 phenylalanine-rich protein), came from a wide range of eukaryotic taxa and yielded a denser tree (more branch points within each lineage) than in earlier studies. Evidence obtained from the reconstructed pattern of base substitutions and deletions in these ancestral loci suggests that, during the early history of the family, selection acted as a transforming force on expressed genes among the duplicates to encode molecular sites with new or modified functions. In later stages of descent, however, selection was a conserving force that preserved the structures of many coadapted functional sites. Each branch of the family was found to have a unique average tempo of evolutionary change, apparently regulated through functional constraints. Proteins whose functions dictate multiple interaction with several other macromolecules evolved more slowly than those which display fewer protein-protein and protein-ion interactions, e.g., calmodulin and next troponin-C evolved at the slowest average rates, whereas parvalbumin evolved at the fastest. The history of all lineages, however, appears to be characterized by rapid rates of evolutionary change in earlier periods, followed by slower rates in more recent periods. A particularly sharp contrast between such fast and slow rates is found in the evolution of calmodulin, whose rate of change in earlier eukaryotes was manyfold faster than the average rate over the past 1 billion years. In fact, the amino acid replacements in the nascent calmodulin lineage occurred at residue positions that in extant metazoans are largely invariable, lending further support to the Darwinian hypothesis that natural selection is both a creative and a conserving force in molecular evolution.      

Two soluble glycosyltransferases glycosylate less efficiently in vivo than their membrane bound counterparts

Many Golgi glycosyltransferases are type II membrane proteins which are cleaved to produce soluble forms that are released from cells. Cho and Cummings recently reported that a soluble form of alpha1, 3- galactosyltransferase was comparable to its membrane bound counterpart in its ability to galactosylate newly synthesized glycoproteins (Cho,S.K. and Cummings,R.D. (1997) J. Biol. Chem., 272, 13622-13628). To test the generality of their findings, we compared the activities of the full length and soluble forms of two such glycosyltransferases, ss1,4 N-Acetylgalactosaminyltransferase (GM2/GD2/ GA2 synthase; GalNAcT) and beta galactoside alpha2,6 sialyltransferase (alpha2,6-ST; ST6Gal I), for production of their glycoconjugate products in vivo . Unlike the full length form of GalNAcT which produced ganglioside GM2 in transfected cells, soluble GalNAcT did not produce detectable GM2 in vivo even though it possessed in vitro GalNAcT activity comparable to that of full length GalNAcT. When compared with cells expressing full length alpha2,6-ST, cells expressing a soluble form of alpha2,6-ST contained 3-fold higher alpha2,6-ST mRNA levels and secreted 7-fold greater alpha2,6-ST activity as measured in vitro , but in striking contrast contained 2- to 4-fold less of the alpha2,6-linked sialic acid moiety in cellular glycoproteins in vivo . In summary these results suggest that unlike alpha1,3-galactosyltransferase the soluble forms of these two glycosyltransferases are less efficient at glycosylation of membrane proteins and lipids in vivo than their membrane bound counterparts.      

Intercalation of proflavine and a platinum derivative of proflavine into double-helical Poly(A)

The equilibria and kinetics of the interactions of proflavine (PR) and its platinum-containing derivative [PtCl(tmen)(2)HNC(13)H(7)(NHCH(2)CH(2))(2)](+) (PRPt) with double-stranded poly(A) have been investigated by spectrophotometry and Joule temperature-jump relaxation at ionic strength 0.1 M, 25 degrees C, and pH 5.2. Spectrophotometric measurements indicate that base-dye interactions are prevailing. T-jump experiments with polarized light showed that effects due to field-induced alignment could be neglected. Both of the investigated systems display two relaxation effects. The kinetic features of the reaction are discussed in terms of a two-step series mechanism in which a precursor complex DS(I) is formed in the fast step, which is then converted to a final complex in the slow step. The rate constants of the fast step are k(1) = (2.5 +/- 0.4) x 10(6) M(-1) s(-1), k(-1) = (2.4 +/- 0.1) x 10(3) s(-1) for poly(A)-PR and k(1) = (2.3 +/- 0.1) x 10(6) M(-1) s(-1), k(-1) = (1.6 +/- 0.2) x 10(3) s(-1) for poly(A)-PRPt. The rate constants for the slow step are k(2) = (4.5 +/- 0.5) x 10(2) s(-1), k(-2) = (1.7 +/- 0.1) x 10(2) s(-1) for poly(A)-PR and k(2) = 9.7 +/- 1.2 s(-1), k(-2) = 10.6 +/- 0.2 s(-1) for poly(A)-PRPt. Spectrophotometric measurements yield for the equilibrium constants and site size the values K = (4.5 +/- 0.1) x 10(3) M(-1), n = 1.3 +/- 0.5 for poly(A)-PR and K = (2.9 +/- 0.1) x 10(3) M(-1), n = 2.3 +/- 0.6 for poly(A)-PRPt. The values of k(1) are similar and lower than expected for diffusion-limited reactions. The values of k(-1) are similar as well. It is suggested that the formation of DS(I) involves only the proflavine residues in both systems. In contrast, the values of k(2) and k(-2) in poly(A)-PRPt are much lower than in poly(A)-PR. The results suggest that in the complex DS(II) of poly(A)-PRPt both proflavine and platinum residues are intercalated. In addition, a very slow process was detected and ascribed to the covalent binding of Pt(II) to the adenine.     

Ceramide inhibits axonal growth and nerve growth factor uptake without compromising the viability of sympathetic neurons.

Ceramide inhibits axonal growth of cultured rat sympathetic neurons when the ceramide content of distal axons, but not cell bodies, is increased (Posse de Chaves, E. I., Bussiere, M. Vance, D. E., Campenot, R. B., and Vance, J.E. (1997) J. Biol. Chem. 272, 3028-3035). We now report that inhibition of growth does not result from cell death since although ceramide is a known apoptotic agent, C(6)-ceramide given to the neurons for 24 h did not cause cell death but instead protected the neurons from death induced by deprivation of nerve growth factor (NGF). We also find that a pool of ceramide generated from sphingomyelin in distal axons, but not cell bodies, inhibits axonal growth. Analysis of endogenous sphingomyelinase activities demonstrated that distal axons are rich in neutral sphingomyelinase activity but contain almost no acidic sphingomyelinase, which is concentrated in cell bodies/proximal axons. Together, these observations are consistent with the idea that generation of ceramide from sphingomyelin by a neutral sphingomyelinase in axons inhibits axonal growth. Furthermore, we demonstrate that treatment of distal axons with ceramide inhibits the uptake of NGF and low density lipoproteins by distal axons by approximately 70 and 40%, respectively, suggesting that the inhibition of axonal growth by ceramide might be due, at least in part, to impaired endocytosis of NGF. However, inhibition of endocytosis of NGF by ceramide could not be ascribed to decreased phosphorylation of TrkA.     

A comparison of adiposity measures as predictors of all-cause mortality: the Melbourne Collaborative Cohort Study

Objective: Our goal was to examine five different measures of adiposity as predictors of all‐cause mortality. Research Methods and Procedures: Subjects were 16,969 men and 24,344 women enrolled between 1990 and 1994 in the Melbourne Collaborative Cohort Study (27 to 75 years of age). There were 2822 deaths over a median follow‐up period of 11 years. BMI, waist circumference, and waist‐to‐hip ratio were obtained from direct anthropometric measurements. Fat mass and percentage fat were estimated by bioelectric impedance analysis. Results: Comparing the top quintile with the second quintile, for men there was an increased risk of between 20% and 30% for all‐cause mortality associated with each of the anthropometric measures. For women, there was an increased risk of 30% (95% confidence interval for hazard ratio, 1.1–1.6) observed for waist circumference and 50% (1.2–1.8) for waist‐to‐hip ratio, but little or no increased risk for BMI, fat mass, and percentage fat. Waist‐to‐hip ratio was positively and monotonically associated with all‐cause mortality for both men and women. There was a linear association between waist circumference and all‐cause mortality for men, whereas a U‐shaped association was observed for women. Discussion: Measures of central adiposity were better predictors of mortality in women in the Melbourne Collaborative Cohort Study compared with measures of overall adiposity. We recommend measuring waist and hip circumferences in population studies investigating the risk of all‐cause mortality associated with obesity. The use of additional measures such as bioelectric impedance is not justified for this outcome.     

Comparison of genomic signatures of selection on Plasmodium falciparum between different regions of a country with high malaria endemicity

Background

Genome wide sequence analyses of malaria parasites from widely separated areas of the world have identified contrasting population structures and signatures of selection. To compare relatively closely situated but ecologically contrasting regions within an endemic African country, population samples of Plasmodium falciparum clinical isolates were collected in Ghana from Kintampo in the central forest-savannah area, and Navrongo in a drier savannah area ~350 km to the north with more seasonally-restricted transmission. Parasite DNA was sequenced and paired-end reads mapped to the P. falciparum reference genome.

Results

High coverage genome wide sequence data for 85 different clinical isolates enabled analysis of 121,712 single nucleotide polymorphisms (SNPs). The local populations had similar proportions of mixed genotype infections, similar SNP allele frequency distributions, and eleven chromosomal regions had elevated integrated haplotype scores (|iHS|) in both. A between-population Rsb metric comparing extended haplotype homozygosity indicated a stronger signal within Kintampo for one of these regions (on chromosome 14) and in Navrongo for two of these regions (on chromosomes 10 and 13). At least one gene in each of these identified regions is a potential target of locally varying selection. The candidates include genes involved in parasite development in mosquitoes, members of variant-expressed multigene families, and a leading vaccine-candidate target of immunity.

Conclusions

Against a background of very similar population structure and selection signatures in the P. falciparum populations of Ghana, three narrow genomic regions showed evidence indicating local differences in historical timing or intensity of selection. Sampling of closely situated populations across heterogeneous environments has potential to refine the mapping of important loci under temporally or spatially varying selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1746-3) contains supplementary material, which is available to authorized users.     

A Genome Wide Association Study of Plasmodium falciparum Susceptibility to 22 Antimalarial Drugs in Kenya

Background

Drug resistance remains a chief concern for malaria control. In order to determine the genetic markers of drug resistant parasites, we tested the genome-wide associations (GWA) of sequence-based genotypes from 35 Kenyan P. falciparum parasites with the activities of 22 antimalarial drugs.

Methods and Principal Findings

Parasites isolated from children with acute febrile malaria were adapted to culture, and sensitivity was determined by in vitro growth in the presence of anti-malarial drugs. Parasites were genotyped using whole genome sequencing techniques. Associations between 6250 single nucleotide polymorphisms (SNPs) and resistance to individual anti-malarial agents were determined, with false discovery rate adjustment for multiple hypothesis testing. We identified expected associations in the pfcrt region with chloroquine (CQ) activity, and other novel loci associated with amodiaquine, quinazoline, and quinine activities. Signals for CQ and primaquine (PQ) overlap in and around pfcrt, and interestingly the phenotypes are inversely related for these two drugs. We catalog the variation in dhfr, dhps, mdr1, nhe, and crt, including novel SNPs, and confirm the presence of a dhfr-164L quadruple mutant in coastal Kenya. Mutations implicated in sulfadoxine-pyrimethamine resistance are at or near fixation in this sample set.

Conclusions/Significance

Sequence-based GWA studies are powerful tools for phenotypic association tests. Using this approach on falciparum parasites from coastal Kenya we identified known and previously unreported genes associated with phenotypic resistance to anti-malarial drugs, and observe in high-resolution haplotype visualizations a possible signature of an inverse selective relationship between CQ and PQ.     

Drug-resistant genotypes and multi-clonality in Plasmodium falciparum analysed by direct genome sequencing from peripheral blood of malaria patients

Naturally acquired blood-stage infections of the malaria parasite Plasmodium falciparum typically harbour multiple haploid clones. The apparent number of clones observed in any single infection depends on the diversity of the polymorphic markers used for the analysis, and the relative abundance of rare clones, which frequently fail to be detected among PCR products derived from numerically dominant clones. However, minority clones are of clinical interest as they may harbour genes conferring drug resistance, leading to enhanced survival after treatment and the possibility of subsequent therapeutic failure. We deployed new generation sequencing to derive genome data for five non-propagated parasite isolates taken directly from 4 different patients treated for clinical malaria in a UK hospital. Analysis of depth of coverage and length of sequence intervals between paired reads identified both previously described and novel gene deletions and amplifications. Full-length sequence data was extracted for 6 loci considered to be under selection by antimalarial drugs, and both known and previously unknown amino acid substitutions were identified. Full mitochondrial genomes were extracted from the sequencing data for each isolate, and these are compared against a panel of polymorphic sites derived from published or unpublished but publicly available data. Finally, genome-wide analysis of clone multiplicity was performed, and the number of infecting parasite clones estimated for each isolate. Each patient harboured at least 3 clones of P. falciparum by this analysis, consistent with results obtained with conventional PCR analysis of polymorphic merozoite antigen loci. We conclude that genome sequencing of peripheral blood P. falciparum taken directly from malaria patients provides high quality data useful for drug resistance studies, genomic structural analyses and population genetics, and also robustly represents clonal multiplicity.     

Comparison of genomic signatures of selection on <Emphasis Type="Italic">Plasmodium falciparum</Emphasis> between different regions of a country with high malaria endemicity

Background

Genome wide sequence analyses of malaria parasites from widely separated areas of the world have identified contrasting population structures and signatures of selection. To compare relatively closely situated but ecologically contrasting regions within an endemic African country, population samples of Plasmodium falciparum clinical isolates were collected in Ghana from Kintampo in the central forest-savannah area, and Navrongo in a drier savannah area ~350 km to the north with more seasonally-restricted transmission. Parasite DNA was sequenced and paired-end reads mapped to the P. falciparum reference genome.

Results

High coverage genome wide sequence data for 85 different clinical isolates enabled analysis of 121,712 single nucleotide polymorphisms (SNPs). The local populations had similar proportions of mixed genotype infections, similar SNP allele frequency distributions, and eleven chromosomal regions had elevated integrated haplotype scores (|iHS|) in both. A between-population Rsb metric comparing extended haplotype homozygosity indicated a stronger signal within Kintampo for one of these regions (on chromosome 14) and in Navrongo for two of these regions (on chromosomes 10 and 13). At least one gene in each of these identified regions is a potential target of locally varying selection. The candidates include genes involved in parasite development in mosquitoes, members of variant-expressed multigene families, and a leading vaccine-candidate target of immunity.

Conclusions

Against a background of very similar population structure and selection signatures in the P. falciparum populations of Ghana, three narrow genomic regions showed evidence indicating local differences in historical timing or intensity of selection. Sampling of closely situated populations across heterogeneous environments has potential to refine the mapping of important loci under temporally or spatially varying selection.