Multigenomic Delineation of Plasmodium Species of the Laverania Subgenus Infecting Wild-living Chimpanzees and Gorillas

Plasmodium falciparum, the major cause of malaria morbidity and mortality worldwide, is only distantly related to other human malaria parasites and has thus been placed in a separate subgenus, termed Laverania. Parasites morphologically similar to P. falciparum have been identified in African apes, but only one other Laverania species, Plasmodium reichenowi from chimpanzees, has been formally described. Although recent studies have pointed to the existence of additional Laverania species, their precise number and host associations remain uncertain, primarily because of limited sampling and a paucity of parasite sequences other than from mitochondrial DNA. To address this, we used limiting dilution polymerase chain reaction to amplify additional parasite sequences from a large number of chimpanzee and gorilla blood and fecal samples collected at two sanctuaries and 30 field sites across equatorial Africa. Phylogenetic analyses of more than 2,000 new sequences derived from the mitochondrial, nuclear, and apicoplast genomes revealed six divergent and well-supported clades within the Laverania parasite group. Although two of these clades exhibited deep subdivisions in phylogenies estimated from organelle gene sequences, these sublineages were geographically defined and not present in trees from four unlinked nuclear loci. This greatly expanded sequence data set thus confirms six, and not seven or more, ape Laverania species, of which P. reichenowi, Plasmodium gaboni, and Plasmodium billcollinsi only infect chimpanzees, whereas Plasmodium praefalciparum, Plasmodium adleri, and Pladmodium blacklocki only infect gorillas. The new sequence data also confirm the P. praefalciparum origin of human P. falciparum

Adaptive Evolution of RH5 in <i>Ape Plasmodium</i> species of the <i>Laverania</i> Subgenus

Plasmodium falciparum, the major cause of malaria morbidity and mortality in humans, has been shown to have emerged after cross-species transmission of one of six host-specific parasites (subgenus Laverania) infecting wild chimpanzees (Pan troglodytes) and western gorillas (Gorilla gorilla). Binding of the parasite-encoded ligand RH5 to the host protein basigin is essential for erythrocyte invasion and has been implicated in host specificity. A recent study claimed to have found two amino acid changes in RH5 that “drove the host shift leading to the emergence of P. falciparum as a human pathogen.” However, the ape Laverania data available at that time, which included only a single distantly related chimpanzee parasite sequence, were inadequate to justify any such conclusion. Here, we have investigated Laverania Rh5 gene evolution using sequences from all six ape parasite species. Searching for gene-wide episodic selection across the entire Laverania phylogeny, we found eight codons to be under positive selection, including three that correspond to contact residues known to form hydrogen bonds between P. falciparum RH5 and human basigin. One of these sites (residue 197) has changed subsequent to the transmission from apes to humans that gave rise to P. falciparum, suggesting a possible role in the adaptation of the gorilla parasite to the human host. We also found evidence that the patterns of nucleotide polymorphisms in P. falciparum are not typical of Laverania species and likely reflect the recent demographic history of the human parasite

A Proteomic Approach to Study the Effect of Thiotaurine on Human Neutrophil Activation

Thiotaurine, a thiosulfonate related to taurine and hypotaurine, is formed by a metabolic process from cystine and generated by a transulfuration reaction between hypotaurine and thiocysteine. Thiotaurine can produce hydrogen sulfide (H2S) from its sulfane sulfur moiety. H2S is a gaseous signaling molecule which can have regulatory roles in inflammatory process. In addition, sulfane sulfur displays the capacity to reversibly bind to other sulfur atoms. Thiotaurine inhibits PMA-induced activation of human neutrophils, and hinders neutrophil spontaneous apoptosis. Here, we present the results of a proteomic approach to study the possible effects of thiotaurine at protein expression level. Proteome analysis of human neutrophils has been performed comparing protein extracts of resting or PMA-activated neutrophils in presence or in absence of thiotaurine. In particular, PMA-stimulated neutrophils showed high level of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) expression compared to the level of the same glycolytic enzyme in the resting neutrophils. Conversely, decreased expression of GAPDH has been observed when human neutrophils were incubated with 1 mM thiotaurine before activation with PMA. This result, confirmed by Western blot analysis, suggests again that thiotaurine shows a bioactive role in the mechanisms underlying the inflammatory process, influencing the energy metabolism of activated leukocytes and raises the possibility that thiotaurine, acting as a sulfur donor, could modulate neutrophil activation via persulfidation of target proteins, such as GAPDH

European American Stratification in Ovarian Cancer Case Control Data: The Utility of Genome-Wide Data for Inferring Ancestry

We investigated the ability of several principal components analysis (PCA)-based strategies to detect and control for population stratification using data from a multi-center study of epithelial ovarian cancer among women of European-American ethnicity. These include a correction based on an ancestry informative markers (AIMs) panel designed to capture European ancestral variation and corrections utilizing un-thinned genome-wide SNP data; case-control samples were drawn from four geographically distinct North-American sites. The AIMs-only and genome-wide first principal components (PC1) both corresponded to the previously described North or Northwest-Southeast axis of European variation. We found that the genome-wide PCA captured this primary dimension of variation more precisely and identified additional axes of genome-wide variation of relevance to epithelial ovarian cancer. Associations evident between the genome-wide PCs and study site corroborate North American immigration history and suggest that undiscovered dimensions of variation lie within Northern Europe. The structure captured by the genome-wide PCA was also found within control individuals and did not reflect the case-control variation present in the data. The genome-wide PCA highlighted three regions of local LD, corresponding to the lactase (LCT) gene on chromosome 2, the human leukocyte antigen system (HLA) on chromosome 6 and to a common inversion polymorphism on chromosome 8. These features did not compromise the efficacy of PCs from this analysis for ancestry control. This study concludes that although AIMs panels are a cost-effective way of capturing population structure, genome-wide data should preferably be used when available

Inherited variation in immune genes and pathways and glioblastoma risk

To determine whether inherited variations in immune function single-nucleotide polymorphisms (SNPs), genes or pathways affect glioblastoma risk, we analyzed data from recent genome-wide association studies in conjunction with predefined immune function genes and pathways. Gene and pathway analyses were conducted on two independent data sets using 6629 SNPs in 911 genes on 17 immune pathways from 525 glioblastoma cases and 602 controls from the University of California, San Francisco (UCSF) and a subset of 6029 SNPs in 893 genes from 531 cases and 1782 controls from MD Anderson (MDA). To further assess consistency of SNP-level associations, we also compared data from the UK (266 cases and 2482 controls) and the Mayo Clinic (114 cases and 111 controls). Although three correlated epidermal growth factor receptor (EGFR) SNPs were consistently associated with glioblastoma in all four data sets (Mantel–Haenzel P values = 1 × 10−5 to 4 × 10−3), independent replication is required as genome-wide significance was not attained. In gene-level analyses, eight immune function genes were significantly (minP < 0.05) associated with glioblastoma; the IL-2RA (CD25) cytokine gene had the smallest minP values in both UCSF (minP = 0.01) and MDA (minP = 0.001) data sets. The IL-2RA receptor is found on the surface of regulatory T cells potentially contributing to immunosuppression characteristic of the glioblastoma microenvironment. In pathway correlation analyses, cytokine signaling and adhesion–extravasation–migration pathways showed similar associations with glioblastoma risk in both MDA and UCSF data sets. Our findings represent the first systematic description of immune genes and pathways that characterize glioblastoma risk

Transcriptome Profiling of Whole Blood Cells Identifies PLEK2 and C1QB in Human Melanoma

Developing analytical methodologies to identify biomarkers in easily accessible body fluids is highly valuable for the early diagnosis and management of cancer patients. Peripheral whole blood is a "nucleic acid-rich" and "inflammatory cell-rich" information reservoir and represents systemic processes altered by the presence of cancer cells.We conducted transcriptome profiling of whole blood cells from melanoma patients. To overcome challenges associated with blood-based transcriptome analysis, we used a PAXgene™ tube and NuGEN Ovation™ globin reduction system. The combined use of these systems in microarray resulted in the identification of 78 unique genes differentially expressed in the blood of melanoma patients. Of these, 68 genes were further analyzed by quantitative reverse transcriptase PCR using blood samples from 45 newly diagnosed melanoma patients (stage I to IV) and 50 healthy control individuals. Thirty-nine genes were verified to be differentially expressed in blood samples from melanoma patients. A stepwise logit analysis selected eighteen 2-gene signatures that distinguish melanoma from healthy controls. Of these, a 2-gene signature consisting of PLEK2 and C1QB led to the best result that correctly classified 93.3% melanoma patients and 90% healthy controls. Both genes were upregulated in blood samples of melanoma patients from all stages. Further analysis using blood fractionation showed that CD45(-) and CD45(+) populations were responsible for the altered expression levels of PLEK2 and C1QB, respectively.The current study provides the first analysis of whole blood-based transcriptome biomarkers for malignant melanoma. The expression of PLEK2, the strongest gene to classify melanoma patients, in CD45(-) subsets illustrates the importance of analyzing whole blood cells for biomarker studies. The study suggests that transcriptome profiling of blood cells could be used for both early detection of melanoma and monitoring of patients for residual disease

De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum)

<p>Abstract</p> <p>Background</p> <p>Transcriptome sequencing data has become an integral component of modern genetics, genomics and evolutionary biology. However, despite advances in the technologies of DNA sequencing, such data are lacking for many groups of living organisms, in particular, many plant taxa. We present here the results of transcriptome sequencing for two closely related plant species. These species, <it>Fagopyrum esculentum </it>and <it>F. tataricum</it>, belong to the order Caryophyllales - a large group of flowering plants with uncertain evolutionary relationships. <it>F. esculentum </it>(common buckwheat) is also an important food crop. Despite these practical and evolutionary considerations <it>Fagopyrum </it>species have not been the subject of large-scale sequencing projects.</p> <p>Results</p> <p>Normalized cDNA corresponding to genes expressed in flowers and inflorescences of <it>F. esculentum </it>and <it>F. tataricum </it>was sequenced using the 454 pyrosequencing technology. This resulted in 267 (for <it>F. esculentum</it>) and 229 (<it>F. tataricum</it>) thousands of reads with average length of 341-349 nucleotides. <it>De novo </it>assembly of the reads produced about 25 thousands of contigs for each species, with 7.5-8.2× coverage. Comparative analysis of two transcriptomes demonstrated their overall similarity but also revealed genes that are presumably differentially expressed. Among them are retrotransposon genes and genes involved in sugar biosynthesis and metabolism. Thirteen single-copy genes were used for phylogenetic analysis; the resulting trees are largely consistent with those inferred from multigenic plastid datasets. The sister relationships of the Caryophyllales and asterids now gained high support from nuclear gene sequences.</p> <p>Conclusions</p> <p>454 transcriptome sequencing and <it>de novo </it>assembly was performed for two congeneric flowering plant species, <it>F. esculentum </it>and <it>F. tataricum</it>. As a result, a large set of cDNA sequences that represent orthologs of known plant genes as well as potential new genes was generated.</p

Obtaining EQ-5D-5L utilities from the disease specific quality of life Alzheimer’s disease scale: development and results from a mapping study

Purpose The Quality of Life Alzheimer’s Disease Scale (QoL-AD) is commonly used to assess disease specific health-related quality of life (HRQoL) as rated by patients and their carers. For cost-effectiveness analyses, utilities based on the EQ-5D are often required. We report a new mapping algorithm to obtain EQ-5D indices when only QoL-AD data are available. Methods Different statistical models to estimate utility directly, or responses to individual EQ-5D questions (response mapping) from QoL-AD, were trialled for patient-rated and proxy-rated questionnaires. Model performance was assessed by root mean square error and mean absolute error. Results The response model using multinomial regression including age and sex, performed best in both the estimation dataset and an independent dataset. Conclusions The recommended mapping algorithm allows researchers for the first time to estimate EQ-5D values from QoL-AD data, enabling cost-utility analyses using datasets where the QoL-AD but no utility measures were collected