Incongruence between test statistics and P values in medical papers

Given an observed test statistic and its degrees of freedom, one may compute the observed P value with most statistical packages. It is unknown to what extent test statistics and P values are congruent in published medical papers. Methods:We checked the congruence of statistical results reported in all the papers of volumes 409–412 of Nature (2001) and a random sample of 63 results from volumes 322–323 of BMJ (2001). We also tested whether the frequencies of the last digit of a sample of 610 test statistics deviated from a uniform distribution (i.e., equally probable digits).Results: 11.6% (21 of 181) and 11.1% (7 of 63) of the statistical results published in Nature and BMJ respectively during 2001 were incongruent, probably mostly due to rounding, transcription, or type-setting errors. At least one such error appeared in 38% and 25% of the papers of Nature and BMJ, respectively. In 12% of the cases, the significance level might change one or more orders of magnitude. The frequencies of the last digit of statistics deviated from the uniform distribution and suggested digit preference in rounding and reporting.Conclusions: this incongruence of test statistics and P values is another example that statistical practice is generally poor, even in the most renowned scientific journals, and that quality of papers should be more controlled and valued

HLA Alleles Influence the Clinical Signature of Amoxicillin-Clavulanate Hepatotoxicity.

BACKGROUND AND AIM The genotype-phenotype interaction in drug-induced liver injury (DILI) is a subject of growing interest. Previous studies have linked amoxicillin-clavulanate (AC) hepatotoxicity susceptibility to specific HLA alleles. In this study we aimed to examine potential associations between HLA class I and II alleles and AC DILI with regards to phenotypic characteristics, severity and time to onset in Spanish AC hepatotoxicity cases. METHODS High resolution genotyping of HLA loci A, B, C, DRB1 and DQB1 was performed in 75 AC DILI cases and 885 controls. RESULTS The distributions of class I alleles A*3002 (P/Pc = 2.6E-6/5E-5, OR 6.7) and B*1801 (P/Pc = 0.008/0.22, OR 2.9) were more frequently found in hepatocellular injury cases compared to controls. In addition, the presence of the class II allele combination DRB1*1501-DQB1*0602 (P/Pc = 5.1E-4/0.014, OR 3.0) was significantly increased in cholestatic/mixed cases. The A*3002 and/or B*1801 carriers were found to be younger (54 vs 65 years, P = 0.019) and were more frequently hospitalized than the DRB1*1501-DQB1*0602 carriers. No additional alleles outside those associated with liver injury patterns were found to affect potential severity as measured by Hy's Law criteria. The phenotype frequencies of B*1801 (P/Pc = 0.015/0.42, OR 5.2) and DRB1*0301-DQB1*0201 (P/Pc = 0.0026/0.07, OR 15) were increased in AC DILI cases with delayed onset compared to those corresponding to patients without delayed onset, while the opposite applied to DRB1*1302-DQB1*0604 (P/Pc = 0.005/0.13, OR 0.07). CONCLUSIONS HLA class I and II alleles influence the AC DILI signature with regards to phenotypic expression, latency presentation and severity in Spanish patients.

Human genetic selection on the MTHFR 677C>T polymorphism.

BACKGROUND The prevalence of genotypes of the 677C>T polymorphism for the MTHFR gene varies among humans. In previous studies, we found changes in the genotypic frequencies of this polymorphism in populations of different ages, suggesting that this could be caused by an increase in the intake of folate and multivitamins by women during the periconceptional period. The aim was to analyze changes in the allelic frequencies of this polymorphism in a Spanish population, including samples from spontaneous abortions (SA). METHODS A total of 1305 subjects born in the 20th century were genotyped for the 677C>T polymorphism using allele specific real-time PCR with Taqman probes. A section of our population (n = 276) born in 1980-1989 was compared with fetal samples (n = 344) from SA of unknown etiology from the same period. RESULTS An increase in the frequency of the T allele (0.38 vs 0.47; p < 0.001) and of the TT genotype (0.14 vs 0.24; p < 0.001) in subjects born in the last quarter of the century was observed. In the 1980-1989 period, the results show that the frequency of the wild type genotype (CC) is about tenfold lower in the SA samples than in the controls (0.03 vs 0.33; p < 0.001) and that the frequency of the TT genotype increases in the controls (0.19 to 0.27) and in the SA samples (0.20 to 0.33 (p < 0.01)); r = 0.98. CONCLUSION Selection in favor of the T allele has been detected. This selection could be due to the increased fetal viability in early stages of embryonic development, as is deduced by the increase of mutants in both living and SA populations.

Polymorphisms in DNA-Repair Genes in a Cohort of Prostate Cancer Patients from Different Areas in Spain: Heterogeneity between Populations as a Confounding Factor in Association Studies.

BACKGROUND Differences in the distribution of genotypes between individuals of the same ethnicity are an important confounder factor commonly undervalued in typical association studies conducted in radiogenomics. OBJECTIVE To evaluate the genotypic distribution of SNPs in a wide set of Spanish prostate cancer patients for determine the homogeneity of the population and to disclose potential bias. DESIGN SETTING AND PARTICIPANTS A total of 601 prostate cancer patients from Andalusia, Basque Country, Canary and Catalonia were genotyped for 10 SNPs located in 6 different genes associated to DNA repair: XRCC1 (rs25487, rs25489, rs1799782), ERCC2 (rs13181), ERCC1 (rs11615), LIG4 (rs1805388, rs1805386), ATM (rs17503908, rs1800057) and P53 (rs1042522). The SNP genotyping was made in a Biotrove OpenArray® NT Cycler. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Comparisons of genotypic and allelic frequencies among populations, as well as haplotype analyses were determined using the web-based environment SNPator. Principal component analysis was made using the SnpMatrix and XSnpMatrix classes and methods implemented as an R package. Non-supervised hierarchical cluster of SNP was made using MultiExperiment Viewer. RESULTS AND LIMITATIONS We observed that genotype distribution of 4 out 10 SNPs was statistically different among the studied populations, showing the greatest differences between Andalusia and Catalonia. These observations were confirmed in cluster analysis, principal component analysis and in the differential distribution of haplotypes among the populations. Because tumor characteristics have not been taken into account, it is possible that some polymorphisms may influence tumor characteristics in the same way that it may pose a risk factor for other disease characteristics. CONCLUSION Differences in distribution of genotypes within different populations of the same ethnicity could be an important confounding factor responsible for the lack of validation of SNPs associated with radiation-induced toxicity, especially when extensive meta-analysis with subjects from different countries are carried out.

Universal complex structures in written language

Quantitative linguistics has provided us with a number of empirical laws that characterise the evolution of languages and competition amongst them. In terms of language usage, one of the most influential results is Zipf’s law of word frequencies. Zipf’s law appears to be universal, and may not even be unique to human language. However, there is ongoing controversy over whether Zipf’s law is a good indicator of complexity. Here we present an alternative approach that puts Zipf’s law in the context of critical phenomena (the cornerstone of complexity in physics) and establishes the presence of a large-scale “attraction” between successive repetitions of words. Moreover, this phenomenon is scale-invariant and universal – the pattern is independent of word frequency and is observed in texts by different authors and written in different languages. There is evidence, however, that the shape of the scaling relation changes for words that play a key role in the text, implying the existence of different “universality classes” in the repetition of words. These behaviours exhibit striking parallels with complex catastrophic phenomena.

KAM theory for conformally symplectic systems

We present a KAM theory for some dissipative systems (geometrically, these are conformally symplectic systems, i.e. systems that transform a symplectic form into a multiple of itself). For systems with n degrees of freedom depending on n parameters we show that it is possible to find solutions with n-dimensional (Diophantine) frequencies by adjusting the parameters. We do not assume that the system is close to integrable, but we use an a-posteriori format. Our unknowns are a parameterization of the solution and a parameter. We show that if there is a sufficiently approximate solution of the invariance equation, which also satisfies some explicit non–degeneracy conditions, then there is a true solution nearby. We present results both in Sobolev norms and in analytic norms. The a–posteriori format has several consequences: A) smooth dependence on the parameters, including the singular limit of zero dissipation; B) estimates on the measure of parameters covered by quasi–periodic solutions; C) convergence of perturbative expansions in analytic systems; D) bootstrap of regularity (i.e., that all tori which are smooth enough are analytic if the map is analytic); E) a numerically efficient criterion for the break–down of the quasi–periodic solutions. The proof is based on an iterative quadratically convergent method and on suitable estimates on the (analytical and Sobolev) norms of the approximate solution. The iterative step takes advantage of some geometric identities, which give a very useful coordinate system in the neighborhood of invariant (or approximately invariant) tori. This system of coordinates has several other uses: A) it shows that for dissipative conformally symplectic systems the quasi–periodic solutions are attractors, B) it leads to efficient algorithms, which have been implemented elsewhere. Details of the proof are given mainly for maps, but we also explain the slight modifications needed for flows and we devote the appendix to present explicit algorithms for flows.

Intracellular excision and reintegration dynamics of the ICEclc genomic island of Pseudomonas knackmussii sp. strain B13.

Genomic islands are DNA elements acquired by horizontal gene transfer that are common to a large number of bacterial genomes, which can contribute specific adaptive functions, e.g. virulence, metabolic capacities or antibiotic resistances. Some genomic islands are still self-transferable and display an intricate life-style, reminiscent of both bacteriophages and conjugative plasmids. Here we studied the dynamical process of genomic island excision and intracellular reintegration using the integrative and conjugative element ICEclc from Pseudomonas knackmussii B13 as model. By using self-transfer of ICEclc from strain B13 to Pseudomonas putida and Cupriavidus necator as recipients, we show that ICEclc can target a number of different tRNA(Gly) genes in a bacterial genome, but only those which carry the GCC anticodon. Two conditional traps were designed for ICEclc based on the attR sequence, and we could show that ICEclc will insert with different frequencies in such traps producing brightly fluorescent cells. Starting from clonal primary transconjugants we demonstrate that ICEclc is excising and reintegrating at detectable frequencies, even in the absence of recipient. Recombination site analysis provided evidence to explain the characteristics of a larger number of genomic island insertions observed in a variety of strains, including Bordetella petri, Pseudomonas aeruginosa and Burkholderia.

Sequence polymorphisms in Pvs48/45 and Pvs47 gametocyte and gamete surface proteins in Plasmodium vivax isolated in Korea

Nucleotide sequence analyses of the Pvs48/45 and Pvs47 genes were conducted in 46 malaria patients from the Republic of Korea (ROK) (n = 40) and returning travellers from India (n = 3) and Indonesia (n = 3). The domain structures, which were based on cysteine residue position and secondary protein structure, were similar between Plasmodium vivax (Pvs48/45 and Pvs47) and Plasmodium falciparum (Pfs48/45 and Pfs47). In comparison to the Sal-1 reference strain (Pvs48/45, PVX_083235 and Pvs47, PVX_083240), Korean isolates revealed seven polymorphisms (E35K, H211N, K250N, D335Y, A376T, I380T and K418R) in Pvs48/45. These isolates could be divided into five haplotypes with the two major types having frequencies of 47.5% and 20%, respectivelfy. In Pvs47, 10 polymorphisms (F22L, F24L, K27E, D31N, V230I, M233I, E240D, I262T, I273M and A373V) were found and they could be divided into four haplotypes with one major type having a frequency of 75%. The Pvs48/45 isolates from India showed a unique amino acid substitution site (K26R). Compared to the Sal-1 and ROK isolates, the Pvs47 isolates from travellers returning from India and Indonesia had amino acid substitutions (S57T and I262K). The current data may contribute to the development of the malaria transmission-blocking vaccine in future clinical trials.

Significant genetic and phenotypic changes arising from clonal growth of a single spore of an arbuscular mycorrhizal fungus over multiple generations.

Arbuscular mycorrhizal fungi (AMF) are highly successful plant symbionts. They reproduce clonally producing multinucleate spores. It has been suggested that some AMF harbor genetically different nuclei. However, recent advances in sequencing the Glomus irregulare genome have indicated very low within-fungus polymorphism. We tested the null hypothesis that, with no genetic differences among nuclei, no significant genetic or phenotypic variation would occur among clonal single spore lines generated from one initial AMF spore. Furthermore, no additional variation would be expected in the following generations of single spore lines. Genetic diversity contained in one initial spore repeatedly gave rise to genetically different variants of the fungus with novel phenotypes. The genetic changes represented quantitative changes in allele frequencies, most probably as a result of changes in the frequency of genetic variation partitioned on different nuclei. The genetic and phenotypic variation is remarkable, given that it arose repeatedly from one clonal individual. Our results highlight the dynamic nature of AMF genetics. Even though within-fungus genetic variation is low, some is probably partitioned among nuclei and potentially causes changes in the phenotype. Our results are important for understanding AMF genetics, as well as for researchers and biotechnologists hoping to use AMF genetic diversity for the improvement of AMF inoculum.

Genetic polymorphisms in the glutamate-rich protein of Plasmodium falciparum field isolates from a malaria-endemic area of Brazil

The genetic diversity displayed by Plasmodium falciparum, the most deadly Plasmodium species, is a significant obstacle for effective malaria vaccine development. In this study, we identified genetic polymorphisms in P. falciparum glutamate-rich protein (GLURP), which is currently being tested in clinical trials as a malaria vaccine candidate, from isolates found circulating in the Brazilian Amazon at variable transmission levels. The study was performed using samples collected in 1993 and 2008 from rural villages situated near Porto Velho, in the state of Rondônia. DNA was extracted from 126 P. falciparum-positive thick blood smears using the phenol-chloroform method and subjected to a nested polymerase chain reaction protocol with specific primers against two immunodominant regions of GLURP, R0 and R2. Only one R0 fragment and four variants of the R2 fragment were detected. No differences were observed between the two time points with regard to the frequencies of the fragment variants. Mixed infections were uncommon. Our results demonstrate conservation of GLURP-R0 and limited polymorphic variation of GLURP-R2 in P. falciparum isolates from individuals living in Porto Velho. This is an important finding, as genetic polymorphisms in B and T-cell epitopes could have implications for the immunological properties of the antigen.

The phlebotomine fauna (Diptera: Psychodidae) of Guaraí, state of Tocantins, with an emphasis on the putative vectors of American cutaneous leishmaniasis in rural settlement and periurban areas

Phlebotomine sandflies were captured in rural settlement and periurban areas of the municipality of Guaraí in the state of Tocantins (TO), an endemic area of American cutaneous leishmaniasis (ACL). Forty-three phlebotomine species were identified, nine of which have already been recognised as ACL vectors. Eleven species were recorded for the first time in TO. Nyssomyia whitmani was the most abundant species, followed by Evandromyia bourrouli, Nyssomyia antunesi and Psychodopygus complexus. The Shannon-Wiener diversity index and the evenness index were higher in the rural settlement area than in the periurban area. The evaluation of different ecotopes within the rural area showed the highest frequencies of Ev. bourrouli and Ny. antunesi in chicken coops, whereas Ny. whitmani predominated in this ecotope in the periurban area. In the rural settlement area, Ev. bourrouli was the most frequently captured species in automatic light traps and Ps. complexus was the most prevalent in Shannon trap captures. The rural settlement environment exhibited greater phlebotomine biodiversity than the periurban area. Ps. complexus and Psychodopygus ayrozai naturally infected with Leishmania (Viannia) braziliensis were identified. The data identified Ny. whitmani as a potential ACL vector in the periurban area, whereas Ps. complexus was more prevalent in the rural environment associated with settlements.

N-cadherin is essential for retinoic acid-mediated cardiomyogenic differentiation in mouse embryonic stem cells.

Contraction forces developed by cardiomyocytes are transmitted across the plasma membrane through end-to-end connections between the myocytes, called intercalated disks, which enable the coordinated contraction of heart muscle. A component of the intercalated disk, the adherens junction, consists of the cell adhesion molecule, N-cadherin. Embryos lacking N-cadherin die at mid-gestation from cardiovascular abnormalities. We have evaluated the role of N-cadherin in cardiomyogenesis using N-cadherin-null mouse embryonic stem (ES) cells grown as embryoid bodies (EBs) in vitro. Myofibrillogenesis, the spatial orientation of myofibers, and intercellular contacts including desmosomes were normal in N-cadherin-null ES cell-derived cardiomyocytes. The effect of retinoic acid (RA), a stage and dose-dependent cardiogenic factor, was assessed in differentiating ES cells. all-trans (at) RA increased the number of ES cell-derived cardiomyocytes by approximately 3-fold (at 3 x 10(-9) M) in wt EBs. However, this effect was lost in N-cadherin-null EBs. In the presence of supplemented at-RA, the emergence of spontaneously beating cardiomyocytes appeared to be delayed and slightly less efficient in N-cadherin-null compared with wt and heterozygous EBs (frequencies of EBs with beating activity at 5 days: 54+/-18% vs. 96+/-0.5%, and 93+/-7%, respectively; peak frequencies of EBs with beating activity: 83+/-8% vs. 96+/-0.5% and 100%, respectively). In conclusion, cardiomyoyctes differentiating from N-cadherin-null ES cells in vitro show normal myofibrillogenesis and intercellular contacts, but impaired responses to early cardiogenic effects mediated by at-RA. These results suggest that N-cadherin may be essential for RA-induced cardiomyogenesis in mouse ES cells in vitro.

Referrals of cancer versus non-cancer patients to a palliative care consult team: do they differ?

This retrospective study compared 100 consecutive non-cancer (NC) patients referred to a palliative care consult team (PCT) in a Swiss university hospital to 506 cancer (C) patients referred during the same period. The frequencies of reported symptoms were similar in both groups. The main reasons for referral in the NC group were symptom control, global evaluation, and assistance with discharge. Requests for symptom control predominated in the C group. Prior to the first visit, 50% of NC patients were on opioids, compared to 58% of C patients. After the first visit, the proportion of NC patients on opioids increased to 64% and the proportion of C patients to 73%. The median daily oral morphine equivalent dose for NC patients taking opioids prior to the first PCT visit was higher than that for C patients (60 mg versus 45 mg). At the time of death or discharge, the percentage of NC patients on opioids was 64%, while that of C patients was 76%. Moreover, NC patients were on significantly lower median doses of opioids than C patients (31 mg versus 60 mg). Over half the NC patients died during hospitalization, as compared to 33% of C patients. Only 6% of NC patients were discharged to palliative care units, as compared to 22% of C patients.

Parenteral is more efficient than mucosal immunization to induce regression of human papillomavirus-associated genital tumors.

Cervical cancer is a public health concern as it represents the second cause of cancer death in women worldwide. High-risk human papillomaviruses (HPV) are the etiologic agents, and HPV E6 and/or E7 oncogene-specific therapeutic vaccines are under development to treat HPV-related lesions in women. Whether the use of mucosal routes of immunization may be preferable for inducing cell-mediated immune responses able to eradicate genital tumors is still debated because of the uniqueness of the female genital mucosa (GM) and the limited experimentation. Here, we compared the protective activity resulting from immunization of mice via intranasal (i.n.), intravaginal (IVAG) or subcutaneous (s.c.) routes with an adjuvanted HPV type 16 E7 polypeptide vaccine. Our data show that s.c. and i.n. immunizations elicited similar frequencies and avidity of TetE71CD81 and E7-specific Interferon-gamma-secreting cells in the GM, whereas slightly lower immune responses were induced by IVAG immunization. In a novel orthotopic murine model, both s.c. and i.n. immunizations allowed for complete long-term protection against genital E7-expressing tumor challenge. However, only s.c. immunization induced complete regression of already established genital tumors. This suggests that the higher E7-specific systemic response observed after s.c. immunization may contribute to the regression of growing genital tumors, whereas local immune responses may be sufficient to impede genital challenges. Thus, our data show that for an efficiently adjuvanted protein-based vaccine, parenteral vaccination route is superior to mucosal vaccination route for inducing regression of established genital tumors in a murine model of HPV-associated genital cancer.

Cerebrospinal fluid biomarkers of central nervous system dysfunction in HIV-infected patients

Background :¦In addition to opportunistic infections of the central nervous system (CNS), which are due to immunosuppression related to HIV, the latter virus, itself, can cause neuropathological abnormalities which are located mainly in the basal ganglia and are characterized by microglial giant cells, reactive astrocytosis and perivascular monocytes. This HIV encephalopathy is characterized, clinically, by psycho-motor slowing, memory loss, difficulties in complex tasks requiring executive functions, as well as motor disorders .These cognitive deficits are grouped under the acronym of HIV-associated neurocognitive disorders (HAND). In fact, HANDs are subdivided in three groups in accordance with the severity of the cognitive impairment: Asymptomatic Neurocognitive Impairment (ANI), Mild/moderate Neurocognitive Disorders (MND) and HIV Associated Dementia (HAD).¦While the incidence of HAD has significantly decreased in the era of combined antiretrobiral therapy (cART), the prevalence of milder forms of HIV-associated neurocognitive disorders HAND seem to have increased. There are many potential reasons to explain this state of facts.¦An important question is to understand how soon the brain may be affected by HIV. Since performing a biopsy in these patients is not an issue, the study of the CSF represents the best available way to look at putative biomarkers of inflammation/neurodegeneration in the CNS. Here, we wanted to examined the putative usefulness of different biomarkers as early indicators of anti-retroviral failure at the level of the CNS. We chose to study the CSF levels of:¦Amyloid-β 1-42 (Aβ42), Tau total (tTau), phosphorylated Tau (pTau), Neopterin and S100-β.¦Indeed, these molecules are representative biomarkers of the major cells of the CNS, i.e. neurons,¦macrophages/microglia and astrocytes.¦To examine how sensitive were these CSF biomarkers to indicate CNS insults caused by HIV, we proposed to take advantage of the MOST (Monotherapy Switzerland/Thailand study) study, recently published in AIDS. Thus, we collaborated with Prof. Pietro Vernazza in St-Gall. In MOST study, monotherapy (MT) consisting in ritonavir-boosted lopinavir (LPV/r) was compared to continuous conventional antiretroviral therapy including several molecules, hereafter referred as CT¦Methods :We tested 61 cerebrospinal fluid (CSF) samples from 52 patients enrolled in MOST, including 34 CSF samples of CT and 27 of MT (mean duration on MT: 47+20 weeks) in patients who maintained full VL suppression in blood (<50cps/ml). Using enzyme-linked immunosorbent assay (ELISA), we determined the CSF concentration of S100-beta (astrocytosis), neopterin (microglia, inflammation), total Tau (tTau), phosphorylated Tau (pTau), and amyloid-beta 1-42 (Abeta), the latter three markers indicating neuronal damages. The CSF samples of 37 HIV-negative patients with Alzheimer dementia (AD) served as controls. Results are expressed in pg/ml and reported as median ± interquartile range. Mann Whitney-U test was used to compare the results of a given biomarker between two groups and the Fisher test to compare frequencies.¦Results: We found a higher concentration of S100-beta (570±1132) and neopterin (2.5±2.9) in the CSF of MT versus CT (0±532, p=0.002 and 1.2±2.5, p=0.058, respectively). A cutoff of 940 pg/ml for S100-beta allowed to discriminate MT (11 above versus 16 below) from CT (1 vs 33, p=0.0003). At a lesser extent, a cutoff of 11 pg/ml for neopterin separated MT (4 above versus 23) from CT (0 vs 34, p=0.034) (Figure).¦In AD, tTau was higher (270±414) and Abeta lower (234±328) than in CT (150±153, p=0.0078, and 466±489, p=0.007, respectively). Such as for CT, Abeta was lower in AD than in MT (390±412, p=0.01). However, contrasting with CT, the levels of tTau were not different between AD and MT (199±177, p=0.11). S100b (173±214; p=0.0006) and neopterin (1.1±0.9; p=0.0014) were lower in AD than MT.¦Conclusions: Despite full VL-suppression in blood, HIV monotherapy is sufficient to trigger inflammation and, especially, astrocytosis. CSF markers of patients on CT have the same profile as reported for healthy subjects, suggesting that CT permits a good control of HIV in the brain. Finally, the levels of tTau, which are relatively similar between AD and MT patients, suggest that neurons are damaged during monotherapy.