Improved statistical analysis of low abundance phenomena in bimodal bacterial populations.

Accurate detection of subpopulation size determinations in bimodal populations remains problematic yet it represents a powerful way by which cellular heterogeneity under different environmental conditions can be compared. So far, most studies have relied on qualitative descriptions of population distribution patterns, on population-independent descriptors, or on arbitrary placement of thresholds distinguishing biological ON from OFF states. We found that all these methods fall short of accurately describing small population sizes in bimodal populations. Here we propose a simple, statistics-based method for the analysis of small subpopulation sizes for use in the free software environment R and test this method on real as well as simulated data. Four so-called population splitting methods were designed with different algorithms that can estimate subpopulation sizes from bimodal populations. All four methods proved more precise than previously used methods when analyzing subpopulation sizes of transfer competent cells arising in populations of the bacterium Pseudomonas knackmussii B13. The methods' resolving powers were further explored by bootstrapping and simulations. Two of the methods were not severely limited by the proportions of subpopulations they could estimate correctly, but the two others only allowed accurate subpopulation quantification when this amounted to less than 25% of the total population. In contrast, only one method was still sufficiently accurate with subpopulations smaller than 1% of the total population. This study proposes a number of rational approximations to quantifying small subpopulations and offers an easy-to-use protocol for their implementation in the open source statistical software environment R.

Evolutionary branching in deme-structured populations.

Adaptive dynamics shows that a continuous trait under frequency dependent selection may first converge to a singular point followed by spontaneous transition from a unimodal trait distribution into a bimodal one, which is called "evolutionary branching". Here, we study evolutionary branching in a deme-structured population by constructing a quantitative genetic model for the trait variance dynamics, which allows us to obtain an analytic condition for evolutionary branching. This is first shown to agree with previous conditions for branching expressed in terms of relatedness between interacting individuals within demes and obtained from mutant-resident systems. We then show this branching condition can be markedly simplified when the evolving trait affect fecundity and/or survival, as opposed to affecting population structure, which would occur in the case of the evolution of dispersal. As an application of our model, we evaluate the threshold migration rate below which evolutionary branching cannot occur in a pairwise interaction game. This agrees very well with the individual-based simulation results.

Four functionally distinct populations of human effector-memory CD8+ T lymphocytes.

In humans, the pathways of memory and effector T cell differentiation remain poorly defined. We have dissected the functional properties of ex vivo effector-memory (EM) CD45RA-CCR7- T lymphocytes present within the circulating CD8+ T cell pool of healthy individuals. Our studies show that EM T cells are heterogeneous and are subdivided based on differential CD27 and CD28 expression into four subsets. EM(1) (CD27+CD28+) and EM(4) (CD27-CD28+) T cells express low levels of effector mediators such as granzyme B and perforin and high levels of CD127/IL-7Ralpha. EM(1) cells also have a relatively short replicative history and display strong ex vivo telomerase activity. Therefore, these cells are closely related to central-memory (CD45RA-CCR7+) cells. In contrast, EM(2) (CD27+CD28-) and EM(3) (CD27-CD28-) cells express mediators characteristic of effector cells, whereby EM(3) cells display stronger ex vivo cytolytic activity and have experienced larger numbers of cell divisions, thus resembling differentiated effector (CD45RA+CCR7-) cells. These data indicate that progressive up-regulation of cytolytic activity and stepwise loss of CCR7, CD28, and CD27 both characterize CD8+ T cell differentiation. Finally, memory CD8+ T cells not only include central-memory cells but also EM(1) cells, which differ in CCR7 expression and may therefore confer memory functions in lymphoid and peripheral tissues, respectively.

Treatment of acute uncomplicated falciparum malaria with artemether-lumefantrine in nonimmune populations: a safety, efficacy, and pharmacokinetic study.

The efficacy and safety of artemether-lumefantrine for the treatment of malaria in nonimmune populations are not well defined. In this study, 165 nonimmune patients from Europe and non-malarious areas of Colombia with acute, uncomplicated falciparum malaria or mixed infection including P. falciparum were treated with the six-dose regimen of artemether-lumefantrine. The parasitologic cure rate at 28 days was 96.0% for the per protocol population (119/124 patients). Median times to parasite clearance and fever clearance were 41.5 and 36.8 hours, respectively. No patient had gametocytes after Day 7. Treatment was well tolerated; most adverse events were mild to moderate and seemed to be related to malaria. There were few serious adverse events, none of which were considered to be drug-related. No significant effects on ECG or laboratory parameters were observed. In conclusion, the six-dose regimen of artemether-lumefantrine was effective and well tolerated in the treatment of acute uncomplicated falciparum malaria in nonimmune patients.

Spatial analysis of corresponding fingerprint features from match and close non-match populations

The development of statistical models for forensic fingerprint identification purposes has been the subject of increasing research attention in recent years. This can be partly seen as a response to a number of commentators who claim that the scientific basis for fingerprint identification has not been adequately demonstrated. In addition, key forensic identification bodies such as ENFSI [1] and IAI [2] have recently endorsed and acknowledged the potential benefits of using statistical models as an important tool in support of the fingerprint identification process within the ACE-V framework. In this paper, we introduce a new Likelihood Ratio (LR) model based on Support Vector Machines (SVMs) trained with features discovered via morphometric and spatial analyses of corresponding minutiae configurations for both match and close non-match populations often found in AFIS candidate lists. Computed LR values are derived from a probabilistic framework based on SVMs that discover the intrinsic spatial differences of match and close non-match populations. Lastly, experimentation performed on a set of over 120,000 publicly available fingerprint images (mostly sourced from the National Institute of Standards and Technology (NIST) datasets) and a distortion set of approximately 40,000 images, is presented, illustrating that the proposed LR model is reliably guiding towards the right proposition in the identification assessment of match and close non-match populations. Results further indicate that the proposed model is a promising tool for fingerprint practitioners to use for analysing the spatial consistency of corresponding minutiae configurations.

Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe.

Ten microsatellite loci and a partial sequence of the COII mitochondrial gene were used to investigate genetic differentiation in B. terrestris, a bumble bee of interest for its high-value crop pollination. The analysis included eight populations from the European continent, five from Mediterranean islands (six subspecies altogether) and one from Tenerife (initially described as a colour form of B. terrestris but recently considered as a separate species, B. canariensis). Eight of the 10 microsatellite loci displayed high levels of polymorphism in most populations. In B. terrestris populations, the total number of alleles detected per polymorphic locus ranged from 3 to 16, with observed allelic diversity from 3.8 +/- 0.5 to 6.5 +/- 1.4 and average calculated heterozygosities from 0.41 +/- 0.09 to 0.65 +/- 0.07. B. canariensis showed a significantly lower average calculated heterozygosity (0.12 +/- 0.08) and observed allelic diversity (1.5 +/- 0.04) as compared to both continental and island populations of B. terrestris. No significant differentiation was found among populations of B. terrestris from the European continent. In contrast, island populations were all significantly and most of them strongly differentiated from continental populations. B. terrestris mitochondrial DNA is characterized by a low nucleotide diversity: 0.18% +/- 0.07%, 0.20% +/- 0.04% and 0.27% +/- 0.04% for the continental populations, the island populations and all populations together, respectively. The only haplotype found in the Tenerife population differs by a single nucleotide substitution from the most common continental haplotype of B. terrestris. This situation, identical to that of Tyrrhenian islands populations and quite different from that of B. lucorum (15 substitutions between terrestris and lucorum mtDNA) casts doubts on the species status of B. canariensis. The large genetic distance between the Tenerife and B. terrestris populations estimated from microsatellite data result, most probably, from a severe bottleneck in the Canary island population. Microsatellite and mitochondrial DNA data call for the protection of the island populations of B. terrestris against importation of bumble bees of foreign origin which are used as crop pollinators.

Reproductive and post-reproductive life history of wild-caught Drosophila melanogaster under laboratory conditions.

The life history of the fruit fly (Drosophila melanogaster) is well understood, but fitness components are rarely measured by following single individuals over their lifetime, thereby limiting insights into lifetime reproductive success, reproductive senescence and post-reproductive lifespan. Moreover, most studies have examined long-established laboratory strains rather than freshly caught individuals and may thus be confounded by adaptation to laboratory culture, inbreeding or mutation accumulation. Here, we have followed the life histories of individual females from three recently caught, non-laboratory-adapted wild populations of D. melanogaster. Populations varied in a number of life-history traits, including ovariole number, fecundity, hatchability and lifespan. To describe individual patterns of age-specific fecundity, we developed a new model that allowed us to distinguish four phases during a female's life: a phase of reproductive maturation, followed by a period of linear and then exponential decline in fecundity and, finally, a post-ovipository period. Individual females exhibited clear-cut fecundity peaks, which contrasts with previous analyses, and post-peak levels of fecundity declined independently of how long females lived. Notably, females had a pronounced post-reproductive lifespan, which on average made up 40% of total lifespan. Post-reproductive lifespan did not differ among populations and was not correlated with reproductive fitness components, supporting the hypothesis that this period is a highly variable, random 'add-on' at the end of reproductive life rather than a correlate of selection on reproductive fitness. Most life-history traits were positively correlated, a pattern that might be due to genotype by environment interactions when wild flies are brought into a novel laboratory environment but that is unlikely explained by inbreeding or positive mutational covariance caused by mutation accumulation.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Determination of testosterone misuse in sport, an international comparative study

Introduction: Urinary steroid profiling is used in doping controls to detect testosterone abuse. A testosterone over epitestosterone (T/E) ratio exceeding 4.0 is considered as suspicious of testosterone administration, irrespectively of individual heterogeneous factors such as the athlete's ethnicity. A deletion polymorphism in the UGT2B17 gene was demonstrated to account for a significant part of the inter-individual variability in the T/E between Caucasians and Asians. However, the anti-doping strategy includes the determination of carbon isotope ratio on androgen metabolites which has been demonstrated to be reliable for the direct detection of testosterone misuse. Herein, we examined the profiles and the variability in the 13C/12Cratios of urinary steroids in a widely heterogeneous cohort of professional soccer players residing in different world countries (Argentina, Italy, Japan, South-Africa, Switzerland and Uganda). Aim: The determination of threshold values based on genotype information and diet specific of the ethnicity is expected to enhance significantly the detection of testosterone misuse. Methods: The steroid profile of 57 Africans, 32 Asians, 50 Caucasians and 32 Hispanics was determined by gas chromatography-mass spectrometry. The carbon isotope ratio of selected androgens in urine specimens were determined by means of gas chromatography/combustion/isotope ratio mass spectrometry (GC-C-IRMS). Results: Significant differences have been observed between all ethnic groups. After estimation of the prevalence of the UGT2B17 deletion/deletion genotype (African:22%; Asian:81%; Caucasian:10%; Hispanic:7%), ethnicspecific thresholds were developed for a specificity of 99% for the T/E (African:5.6; Asian:3.8; Caucasian:5.7; Hispanic:5.8). Italian and Swiss populations recorded an enrichment in 13C of the urinary steroids with respect to the other groups, thereby supporting consumption of a relatively larger proportion of C3 plants in their diet. Noteworthy, detection criteria based on the difference in the carbon isotope ratio of androsterone and pregnanediol for each population were well below the established threshold value for positive cases. Conclusion: These profiling results demonstrate that a unique and nonspecific threshold to evidence testosterone misuse is not fit for purpose. In addition, the carbon isotopic ratio from these different diet groups highlight the importance to adapt the criteria for increasing the sensitivity in the detection of exogenous testosterone. In conclusion, it may be emphasized that combining the use of isotope ratio mass spectrometry including refined interpretation criteria for positivity and the subject-based profiling of steroids will most probably improve the efficiency of the confirmatory test.

The adaptive dynamics of niche constructing traits in spatially subdivided populations: evolving posthumous extended phenotypes.

Niche construction, by which organisms modify the environment in which they live, has been proposed to affect the evolution of many phenotypic traits. But what about the evolution of a niche constructing trait itself, whose expression changes the pattern of natural selection to which the trait is exposed in subsequent generations? This article provides an inclusive fitness analysis of selection on niche constructing phenotypes, which can affect their environment from local to global scales in arbitrarily spatially subdivided populations. The model shows that phenotypic effects of genes extending far beyond the life span of the actor can be affected by natural selection, provided they modify the fitness of those individuals living in the future that are likely to have inherited the niche construction lineage of the actor. Present benefits of behaviors are thus traded off against future indirect costs. The future costs will generally result from a complicated interplay of phenotypic effects, population demography and environmental dynamics. To illustrate these points, I derive the adaptive dynamics of a trait involved in the consumption of an abiotic resource, where resource abundance in future generations feeds back to the evolutionary dynamics of the trait.

Two specific populations of GABAergic neurons originating from the medial and the caudal ganglionic eminences aid in proper navigation of callosal axons.

The corpus callosum (CC) plays a crucial role in interhemispheric communication. It has been shown that CC formation relies on the guidepost cells located in the midline region that include glutamatergic and GABAergic neurons as well as glial cells. However, the origin of these guidepost GABAergic neurons and their precise function in callosal axon pathfinding remain to be investigated. Here, we show that two distinct GABAergic neuronal subpopulations converge toward the midline prior to the arrival of callosal axons. Using in vivo and ex vivo fate mapping we show that CC GABAergic neurons originate in the caudal and medial ganglionic eminences (CGE and MGE) but not in the lateral ganglionic eminence (LGE). Time lapse imaging on organotypic slices and in vivo analyses further revealed that CC GABAergic neurons contribute to the normal navigation of callosal axons. The use of Nkx2.1 knockout (KO) mice confirmed a role of these neurons in the maintenance of proper behavior of callosal axons while growing through the CC. Indeed, using in vitro transplantation assays, we demonstrated that both MGE- and CGE-derived GABAergic neurons exert an attractive activity on callosal axons. Furthermore, by combining a sensitive RT-PCR technique with in situ hybridization, we demonstrate that CC neurons express multiple short and long range guidance cues. This study strongly suggests that MGE- and CGE-derived interneurons may guide CC axons by multiple guidance mechanisms and signaling pathways. © 2013 Wiley Periodicals, Inc. Develop Neurobiol 73: 647-672, 2013.