Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58 Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates. (C) 2009 Elsevier B.V. All rights reserved.

Recurrent Parasitemias and Population Dynamics of Plasmodium vivax Polymorphisms in Rural Amazonia

Clinical trials documented alarming post-treatment Plasmodium vivax recurrence rates caused by recrudescence of surviving asexual blood stages, relapse from hypnozoites, or new infections. Here we describe high rates of P vivax recurrence (26-40% 180 days after treatment) in two cohorts of rural Amazonians exposed to low levels of malaria transmission after a vivax malaria episode treated with chloroquine-primaquine. Microsatellite analysis of 28 paired acute infection and recurrence parasites showed only two pairs of identical haplotypes (consistent with recrudescences or reactivation of homologous hypnozoites) and four pairs of related haplotypes (sharing alleles at 11-13 of 14 microsatellites analyzed). Local isolates of P vivax were extraordinarily diverse and rarely shared the same haplotype, indicating that frequent recurrences did not favor the persistence or reappearance of clonal lineages of parasites in the Population. This fast haplotype replacement rate may represent the typical population dynamics Of neutral polymorphisms in parasites from low-endemicity areas.

Worldwide sequence conservation of transmission-blocking vaccine candidate Pvs230 in Plasmodium vivax

Pfs230, surface protein of gametocyte/gamete of the human malaria parasite, Plasmodium falciparum, is a prime candidate of malaria transmission-blocking vaccine. Plasmodium vivax has an ortholog of Pfs230 (Pvs230), however, there has been no study in any aspects on Pvs230 to date. To investigate whether Pvs230 can be a vivax malaria transmission-blocking vaccine, we performed evolutionary and population genetic analysis of the Pvs230 gene (pvs230: PVX_003905). Our analysis of Pvs230 and its orthologs in eight Plasmodium species revealed two distinctive parts: an interspecies variable part (IVP) containing species-specific oligopeptide repeats at the N-terminus and a 7.5 kb interspecies conserved part (ICP) containing 14 cysteine-rich domains. Pvs230 was closely related to its orthologs, Pks230 and Pcys230, in monkey malaria parasites. Analysis of 113 pvs230 sequences obtained from worldwide, showed that nucleotide diversity is remarkably low in the non-repeat 8-kb region of pvs230 (theta pi = 0.00118) with 77 polymorphic nucleotide sites, 40 of which results in amino acid replacements. A signature of purifying selection but not of balancing selection was seen on pvs230. Functional and/or structural constraints may limit the level of polymorphism in pvs230. The observed limited polymorphism in pvs230 should ground for utilization of Pvs230 as an effective transmission-blocking vaccine. (C) 2011 Elsevier Ltd. All rights reserved.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7 Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5` and 3` segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition. we uncovered two surprising patterns of evolution in repetitive sequence. Firsts in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi. (C) 2009 Elsevier B.V. All rights reserved.

Sequence-specific reconstruction from fragmentary databases using seed sequences: implementation and validation on SAGE, proteome and generic sequencing data

Motivation: DNA assembly programs classically perform an all-against-all comparison of reads to identify overlaps, followed by a multiple sequence alignment and generation of a consensus sequence. If the aim is to assemble a particular segment, instead of a whole genome or transcriptome, a target-specific assembly is a more sensible approach. GenSeed is a Perl program that implements a seed-driven recursive assembly consisting of cycles comprising a similarity search, read selection and assembly. The iterative process results in a progressive extension of the original seed sequence. GenSeed was tested and validated on many applications, including the reconstruction of nuclear genes or segments, full-length transcripts, and extrachromosomal genomes. The robustness of the method was confirmed through the use of a variety of DNA and protein seeds, including short sequences derived from SAGE and proteome projects.

Intraspecific sequence variation in 16S rRNA gene of Ureaplasma diversum isolates

Ureaplasma diversum infection in bulls may result in seminal vesiculitis, balanoposthitis and alterations in spermatozoids. In cows, it can cause placentitis, fetal alveolitis, abortion and the birth of weak calves. U. diversum ATCC 49782 (serogroups A), ATCC 49783 (serogroup C) and 34 field isolates were used for this study. These microorganisms were submitted to Polymerase Chain Reaction for 16S gene sequence determination using Tact High Fidelity and the products were purified and bi-directionally sequenced. Using the sequence obtained, a fragment containing four hypervariable regions was selected and nucleotide polymorphisms were identified based on their position within the 16S rRNA gene. Forty-four single nucleotide polymorphisms (SNP) were detected. The genotypic variability of the 16S rRNA gene of U. diversum isolates shows that the taxonomy classification of these organisms is likely much more complex than previously described and that 16S rRNA gene sequencing may be used to suggest an epidemiologic pattern of different origin strains. (c) 2011 Elsevier B.V. All rights reserved.

Sequence and Structural Analysis of the 5 ` Noncoding Region of Hepatitis C Virus in Patients With Chronic Infection

Hepatitis C virus (HCV), exhibits considerable genetic diversity, but presents a relatively well conserved 5 ` noncoding region (5 ` NCR) among all genotypes. In this study, the structural features and translational efficiency of the HCV 5 ` NCR sequences were analyzed using the programs RNAfold, RNAshapes and RNApdist and with a bicistronic dual luciferase expression system, respectively. RNA structure prediction software indicated that base substitutions will alter potentially the 5 ` NCR structure. The heterogeneous sequence observed on 5 ` NCR led to important changes in their translation efficiency in different cell culture lines. Interactions of the viral RNA with cellular transacting factors may vary according to the cell type and viral genome polymorphisms that may result in the translational efficiency observed. J. Med. Virol. 81: 1212-1219, 2009. (C) 2009 Wiley-Liss, Inc.

Braid groups of non-orientable surfaces and the Fadell-Neuwirth short exact sequence

Let M be a compact, connected non-orientable surface without boundary and of genus g >= 3. We investigate the pure braid groups P,(M) of M, and in particular the possible splitting of the Fadell-Neuwirth short exact sequence 1 -> P(m)(M \ {x(1), ..., x(n)}) hooked right arrow P(n+m)(M) (P*) under right arrow P(n)(M) -> 1, where m, n >= 1, and p* is the homomorphism which corresponds geometrically to forgetting the last m strings. This problem is equivalent to that of the existence of a section for the associated fibration p: F(n+m)(M) -> F(n)(M) of configuration spaces, defined by p((x(1), ..., x(n), x(n+1), ..., x(n+m))) = (x(1), ..., x(n)). We show that p and p* admit a section if and only if n = 1. Together with previous results, this completes the resolution of the splitting problem for surface pure braid groups. (C) 2009 Elsevier B.V. All rights reserved.

Acanthoscurrin Fragment 101-132: Total Synthesis at 60 degrees C of a Novel Difficult Sequence

Glycine-rich proteins (GRP), serve a variety of biological functions. Acanthoscurrin is an antimicrobial GRP isolated front hemocytes-of the Brazilian spider Acanthoscurria gomesiana. Aiming to contribute to the knowledge of the secondary structure and stepwise solid-phase synthesis of GRPs` glycine-rich domains, we attempted to prepare G(101)GGLGGGRGGGYG(113) GGGGYGGGYG(123)GGy(126)GGGKYK(132)-NH(2), acanthoscurrin C-terminal amidated fragment. Although a theoretical prediction did not indicate high aggregation potential for this peptide, repetitive incomplete aminoacylations were observed after incorporating Tyr(126) to the growing peptide-MBHA resin (Boc chemistry) at 60 degrees C. The problem was not solved by varying the coupling reagents or solvents, adding chaotropic salts to the reaction media or changing the resin/chemistry (Rink amide resin/Fmoc chemistry). Some improvement was mode when CLEAR amide resin (Fmoc chemistry) was 32 used, as it allowed for obtaining fragment (G(113)-K(132) NIR-FT-Raman spectra collected for samples of the growing peptide-MBHA, -Rink amide resin and -CLEAR amide resin revealed the presence of beta-sheet structures. Only the combination of CLEAR-amide resin, 60 degrees C, Fmoc-(Fmoc-Hmb)Gly-OH and LiCl (the last two used alternately) was able to inhibit the phenomenon, as proven by NIR-FT-Raman analysis of the growing peptide-resin, allowing the total synthesis of desired 132 fragment Gly(101)-K(132). In summary, this work describes a new difficult sequence, contributes to understanding stepwise solid-phase synthesis of this type of peptide and shows that, at least while protected and linked to a resin, this GRPs glycine-rich motif presents all early tendency to assume beta-sheet structures. (c) 2008 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 92: 65-75, 2009.

Lithium butylchalcogenolate induced Michael-aldol tandem sequence: easy and rapid access to highly functionalized organochalcogenides and unsaturated compounds

Lithium ""butylchalcogenolates are generated in situ by reacting the elements (S, Se, and Te) with (n)butyl-lithium at 0 degrees C. Reaction of the lithium alkylchalcogenolates with activated alkenes and aldehydes gives the corresponding aldol adducts. The selenium-containing products give Morita-Baylis-Hillman adducts after the oxidation/elimination of the selenoxide. The whole sequence can be performed in a one-pot procedure. (C) 2009 Elsevier Ltd. All rights reserved.

Fast Acquisition of (13)C NMR Spectra using the Steady-state Free Precession Sequence

In this paper, we show that the steady-state free precession sequence can be used to acquire (13)C high-resolution nuclear magnetic resonance spectra and applied to qualitative analysis. The analysis of brucine sample using this sequence with 60 degrees flip angle and time interval between pulses equal to 300 ms (acquisition time, 299.7 ms; recycle delay, 300 ms) resulted in spectrum with twofold enhancement in signal-to-noise ratio, when compared to standard (13)C sequence. This gain was better when a much shorter time interval between pulses (100 ms) was applied. The result obtained was more than fivefold enhancement in signal-to-noise ratio, equivalent to more than 20-fold reduction in total data recording time. However, this short time interval between pulses produces a spectrum with severe phase and truncation anomalies. We demonstrated that these anomalies can be minimized by applying an appropriate apodization function and plotting the spectrum in the magnitude mode.

Risk of recurrent stroke and death after first stroke in long-distance ski race participants

Background: Physical activity is of benefit for primary prevention of cardiovascular diseases, but it appears to increase the risk for atrial fibrillation. We aimed to study a cohort of patients following a first stroke in individuals with previous high physical activity, compare them to the general population with respect to recurrent stroke and death, and relate these to atrial fibrillation. Methods and results: From the participants of the Vasaloppet, the world's largest ski-race, and matched individuals from the general population (n=708 604), we identified 5964 patients hospitalized with a first-time stroke between 1994 and 2010. Individuals with severe diseases were excluded. One half percent of skiers and 1% of nonskiers were hospitalized due to stroke. The incidence rate was 8.3 per 100 person-years among skiers and 11.1 among nonskiers. The hazard ratio (HR) for recurrent stroke or death between the 2 groups was 0.76 (95% CI 0.67 to 0.86). The result was consistent in subgroups. The HR for death was 0.66 (95% CI 0.56 to 0.78) and for recurrent stroke 0.82 (95% CI 0.70 to 0.96). After adjustment for smoking and socioeconomic factors, the HR for death was consistent at 0.70 (95% CI 0.56 to 0.87) while the HR for recurrent stroke was not statistically significant. Outcomes for skiers with atrial fibrillation tended to show a lower risk than for nonskiers. Conclusions: This large cohort study supports the hypothesis that patients with a stroke and with prior regular physical activity have a lower risk of death, while their risk for recurrent stroke is similar to that of nonskiers. The skiers had a higher incidence of atrial fibrillation, but still no increased risk of recurring stroke.

All maximal-pairs in step-leap representation of melodic sequence

This paper proposes an efficient pattern extraction algorithm that can be applied on melodic sequences that are represented as strings of abstract intervallic symbols; the melodic representation introduces special “binary don’t care” symbols for intervals that may belong to two partially overlapping intervallic categories. As a special case the well established “step–leap” representation is examined. In the step–leap representation, each melodic diatonic interval is classified as a step (±s), a leap (±l) or a unison (u). Binary don’t care symbols are used to represent the possible overlapping between the various abstract categories e.g. *=s, *=l and #=-s, #=-l. We propose an O(n+d(n-d)+z)-time algorithm for computing all maximal-pairs in a given sequence x=x[1..n], where x contains d occurrences of binary don’t cares and z is the number of reported maximal-pairs.