The likelihood approach to compare populations : a study on DNA evidence and pitfalls of intuitions

The paper follows on from earlier work [Taroni F and Aitken CGG. Probabilistic reasoning in the law, Part 1: assessment of probabilities and explanation of the value of DNA evidence. Science & Justice 1998; 38: 165-177]. Different explanations of the value of DNA evidence were presented to students from two schools of forensic science and to members of fifteen laboratories all around the world. The responses were divided into two groups; those which came from a school or laboratory identified as Bayesian and those which came from a school or laboratory identified as non-Bayesian. The paper analyses these responses using a likelihood approach. This approach is more consistent with a Bayesian analysis than one based on a frequentist approach, as was reported by Taroni F and Aitken CGG. [Probabilistic reasoning in the law, Part 1: assessment of probabilities and explanation of the value of DNA evidence] in Science & Justice 1998.

Cost-effective experimental design to support modeling of concentration-response functions

Modeling concentration-response function became extremely popular in ecotoxicology during the last decade. Indeed, modeling allows determining the total response pattern of a given substance. However, reliable modeling is consuming in term of data, which is in contradiction with the current trend in ecotoxicology, which aims to reduce, for cost and ethical reasons, the number of data produced during an experiment. It is therefore crucial to determine experimental design in a cost-effective manner. In this paper, we propose to use the theory of locally D-optimal designs to determine the set of concentrations to be tested so that the parameters of the concentration-response function can be estimated with high precision. We illustrated this approach by determining the locally D-optimal designs to estimate the toxicity of the herbicide dinoseb on daphnids and algae. The results show that the number of concentrations to be tested is often equal to the number of parameters and often related to the their meaning, i.e. they are located close to the parameters. Furthermore, the results show that the locally D-optimal design often has the minimal number of support points and is not much sensitive to small changes in nominal values of the parameters. In order to reduce the experimental cost and the use of test organisms, especially in case of long-term studies, reliable nominal values may therefore be fixed based on prior knowledge and literature research instead of on preliminary experiments

A score test for non-nested hypotheses with applications to discrete data models

This paper suggests that a convenient score test against non-nested alternatives can be constructed from the linear combination of the likelihood functions of the competing models. It is shown that this procedure is essentially a test for the correct specification of the conditional distribution of the variable of interest.

Use of DNA profiles for investigation using a simulated national DNA database: Part I. Partial SGM Plus profiles.

In traditional criminal investigation, uncertainties are often dealt with using a combination of common sense, practical considerations and experience, but rarely with tailored statistical models. For example, in some countries, in order to search for a given profile in the national DNA database, it must have allelic information for six or more of the ten SGM Plus loci for a simple trace. If the profile does not have this amount of information then it cannot be searched in the national DNA database (NDNAD). This requirement (of a result at six or more loci) is not based on a statistical approach, but rather on the feeling that six or more would be sufficient. A statistical approach, however, could be more rigorous and objective and would take into consideration factors such as the probability of adventitious matches relative to the actual database size and/or investigator's requirements in a sensible way. Therefore, this research was undertaken to establish scientific foundations pertaining to the use of partial SGM Plus loci profiles (or similar) for investigation.

Oligocene CO2 decline promoted C4 photosynthesis in grasses.

C4 photosynthesis is an adaptation derived from the more common C3 photosynthetic pathway that confers a higher productivity under warm temperature and low atmospheric CO2 concentration [1, 2]. C4 evolution has been seen as a consequence of past atmospheric CO2 decline, such as the abrupt CO2 fall 32-25 million years ago (Mya) [3-6]. This relationship has never been tested rigorously, mainly because of a lack of accurate estimates of divergence times for the different C4 lineages [3]. In this study, we inferred a large phylogenetic tree for the grass family and estimated, through Bayesian molecular dating, the ages of the 17 to 18 independent grass C4 lineages. The first transition from C3 to C4 photosynthesis occurred in the Chloridoideae subfamily, 32.0-25.0 Mya. The link between CO2 decrease and transition to C4 photosynthesis was tested by a novel maximum likelihood approach. We showed that the model incorporating the atmospheric CO2 levels was significantly better than the null model, supporting the importance of CO2 decline on C4 photosynthesis evolvability. This finding is relevant for understanding the origin of C4 photosynthesis in grasses, which is one of the most successful ecological and evolutionary innovations in plant history.

DNA evidence, probabilistic evaluation and collaborative tests.

Forensic scientists working in 12 state or private laboratories participated in collaborative tests to improve the reliability of the presentation of DNA data at trial. These tests were motivated in response to the growing criticism of the power of DNA evidence. The experts' conclusions in the tests are presented and discussed in the context of the Bayesian approach to interpretation. The use of a Bayesian approach and subjective probabilities in trace evaluation permits, in an easy and intuitive manner, the integration into the decision procedure of any revision of the measure of uncertainty in the light of new information. Such an integration is especially useful with forensic evidence. Furthermore, we believe that this probabilistic model is a useful tool (a) to assist scientists in the assessment of the value of scientific evidence, (b) to help jurists in the interpretation of judicial facts and (c) to clarify the respective roles of scientists and of members of the court. Respondents to the survey were reluctant to apply this methodology in the assessment of DNA evidence.

Access to specialized pediatric cancer care in Switzerland.

BACKGROUND: Specialized pediatric cancer centers (PCCs) are thought to be essential to obtain state-of-the-art care for children and adolescents. We determined the proportion of childhood cancer patients not treated in a PCC, and described their characteristics and place of treatment. PROCEDURE: The Swiss Childhood Cancer Registry (SCCR) registers all children treated in Swiss PCCs. The regional cancer registries (covering 14/26 cantons) register all cancer patients of a region. The children of the SCCR with data from 7 regions (11 cantons) were compared, using specialized software for record linkage. All children <16 years of age at diagnosis with primary malignant tumors, diagnosed between 1990 and 2004, and living in one of these regions were included in the analysis. RESULTS: 22.1% (238/1,077) of patients recorded in regional registries were not registered in the SCCR. Of these, 15.7% (169/1,077) had never been in a PCC while 6.4% (69/1,077) had been in a PCC but were not registered in the SCCR, due to incomplete data flow. In all diagnostic groups and in all age groups, a certain proportion of children was treated outside a PCC, but this proportion was largest in children suffering from malignant bone tumors/soft tissue sarcomas and from malignant epithelial neoplasms, and in older children. The proportion of patients treated in a PCC increased over the study period (P < 0.0001). CONCLUSIONS: One in six childhood cancer patients in Switzerland was not treated in a PCC. Whether these patients have different treatment outcomes remained unclear.

Ecological components and evolution of selfing in the freshwater snail Galba truncatula.

The reproductive assurance hypothesis emphasizes that self-fertilization should evolve in species with reduced dispersal capability, low population size or experiencing recurrent bottlenecks. Our work investigates the ecological components of the habitats colonized by the snail, Galba truncatula, that may influence the evolution of selfing. Galba truncatula is a preferential selfer inhabiting freshwater habitats, which vary with respect to the degree of permanence. We considered with a population genetic approach the spatial and the temporal degree of isolation of populations of G. truncatula. We showed that patches at distances of only a few meters are highly structured. The effective population sizes appear quite low, in the order of 10 individuals or less. This study indicates that individuals of the species G. truncatula are likely to be alone in a site and have a low probability of finding a partner from a nearby site to reproduce. These results emphasize the advantage of selfing in this species.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

Emergence of young human genes after a burst of retroposition in primates.

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past approximately 63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome-derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.

False phylogenies on wood mice due to cryptic cytochrome-b pseudogene.

The phylogeny and phylogeography of the Old World wood mice (subgenus Sylvaemus, genus Apodemus, Muridae) are well-documented. Nevertheless, the distributions of species, such as A. fulvipectus and A. ponticus remain dubious, as well as their phylogenetic relationships with A. sylvaticus. We analysed samples of Apodemus spp. across Europe using the mitochondrial cytochrome-b gene (cyt-b) and compared the DNA and amino-acid compositions of previously published sequences. The main result stemming from this study is the presence of a well-differentiated lineage of Sylvaemus including samples of various species (A. sylvaticus, A. fulvipectus, A. ponticus) from distant locations, which were revealed to be nuclear copies of the mitochondrial cyt-b. The presence of this cryptic pseudogene in published sequences is supported by different pathways. This has led to important errors in previous molecular trees and hence to partial misinterpretations in the phylogeny of Apodemus.

Effective size of two feral domestic cat populations (Felis catus L): effect of the mating system.

A variety of behavioural traits have substantial effects on the gene dynamics and genetic structure of local populations. The mating system is a plastic trait that varies with environmental conditions in the domestic cat (Felis catus) allowing an intraspecific comparison of the impact of this feature on genetic characteristics of the population. To assess the potential effect of the heterogenity of males' contribution to the next generation on variance effective size, we applied the ecological approach of Nunney & Elam (1994) based upon a demographic and behavioural study, and the genetic 'temporal methods' of Waples (1989) and Berthier et al. (2002) using microsatellite markers. The two cat populations studied were nearly closed, similar in size and survival parameters, but differed in their mating system. Immigration appeared extremely restricted in both cases due to environmental and social constraints. As expected, the ratio of effective size to census number (Ne/N) was higher in the promiscuous cat population (harmonic mean = 42%) than in the polygynous one (33%), when Ne was calculated from the ecological method. Only the genetic results based on Waples' estimator were consistent with the ecological results, but failed to evidence an effect of the mating system. Results based on the estimation of Berthier et al. (2002) were extremely variable, with Ne sometimes exceeding census size. Such low reliability in the genetic results should retain attention for conservation purposes.

Serum sCD163 levels are associated with type 2 diabetes mellitus and are influenced by coffee and wine consumption: results of the Di@bet.es study.

OBJECTIVE Serum levels of soluble TNF-like weak inducer of apoptosis (sTWEAK) and its scavenger receptor CD163 (sCD163) have been linked to insulin resistance. We analysed the usefulness of these cytokines as biomarkers of type 2 diabetes in a Spanish cohort, together with their relationship to food consumption in the setting of the Di@bet.es study. RESEARCH DESIGN AND METHODS This is a cross-sectional, matched case-control study of 514 type 2 diabetes subjects and 517 controls with a Normal Oral Glucose Tolerance Test (NOGTT), using data from the Di@bet.es study. Study variables included clinical and demographic structured survey, food frequency questionnaire and physical examination. Serum concentrations of sTWEAK and sCD163 were measured by ELISA. Linear regression analysis determined which variables were related to sTWEAK and sCD163 levels. Logistic regression analysis was used to estimate odd ratios of presenting type 2 diabetes. RESULTS sCD163 concentrations and sCD163/sTWEAK ratio were 11.0% and 15.0% higher, respectively, (P<0.001) in type 2 diabetes than in controls. Following adjustment for various confounders, the OR for presenting type 2 diabetes in subjects in the highest vs the lowest tertile of sCD163 was [(OR), 2,01 (95%CI, 1,46-2,97); P for trend <0.001]. Coffee and red wine consumption was negatively associated with serum levels of sCD163 (P = 0.0001 and; P = 0.002 for coffee and red wine intake, respectively). CONCLUSIONS High circulating levels of sCD163 are associated with type 2 diabetes in the Spanish population. The association between coffee and red wine intake and these biomarkers deserves further study to confirm its potential role in type 2 diabetes.

Undertreatment of acute pain (oligoanalgesia) and medical practice variation in prehospital analgesia of adult trauma patients: a 10 yr retrospective study.

BACKGROUND: Prehospital oligoanalgesia is prevalent among trauma victims, even when the emergency medical services team includes a physician. We investigated if not only patients' characteristics but physicians' practice variations contributed to prehospital oligoanalgesia. METHODS: Patient records of conscious adult trauma victims transported by our air rescue helicopter service over 10 yr were reviewed retrospectively. Oligoanalgesia was defined as a numeric rating scale (NRS) >3 at hospital admission. Multilevel logistic regression analysis was used to predict oligoanalgesia, accounting first for patient case-mix, and then physician-level clustering. The intraclass correlation was expressed as the median odds ratio (MOR). RESULTS: A total of 1202 patients and 77 physicians were included in the study. NRS at the scene was 6.9 (1.9). The prevalence of oligoanalgesia was 43%. Physicians had a median of 5.7 yr (inter-quartile range: 4.2-7.5) of post-graduate training and 27% were female. In our multilevel analysis, significant predictors of oligoanalgesia were: no analgesia [odds ratio (OR) 8.8], National Advisory Committee for Aeronautics V on site (OR 4.4), NRS on site (OR 1.5 per additional NRS unit >4), female physician (OR 2.0), and years of post-graduate experience [>4.0 to ≤5.0 (OR 1.3), >3.0 to ≤4.0 (OR 1.6), >2.0 to ≤3.0 (OR 2.6), and ≤2.0 yr (OR 16.7)]. The MOR was 2.6, and was statistically significant. CONCLUSIONS: Physicians' practice variations contributed to oligoanalgesia, a factor often overlooked in analyses of prehospital pain management. Further exploration of the sources of these variations may provide innovative targets for quality improvement programmes to achieve consistent pain relief for trauma victims.