THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

Application of detection probabilities to the design of amphibian monitoring programs in temporary ponds

Failure to detect a species in an area where it is present is a major source of error in biological surveys. We assessed whether it is possible to optimize single-visit biological monitoring surveys of highly dynamic freshwater ecosystems by framing them a priori within a particular period of time. Alternatively, we also searched for the optimal number of visits and when they should be conducted. We developed single-species occupancy models to estimate the monthly probability of detection of pond-breeding amphibians during a four-year monitoring program. Our results revealed that detection probability was species-specific and changed among sampling visits within a breeding season and also among breeding seasons. Thereby, the optimization of biological surveys with minimal survey effort (a single visit) is not feasible as it proves impossible to select a priori an adequate sampling period that remains robust across years. Alternatively, a two-survey combination at the beginning of the sampling season yielded optimal results and constituted an acceptable compromise between sampling efficacy and survey effort. Our study provides evidence of the variability and uncertainty that likely affects the efficacy of monitoring surveys, highlighting the need of repeated sampling in both ecological studies and conservation management.

Helical tomotherapy or intensity-modulated radiation therapy in the treatment of anal cancer: experience of Geneva and Lausanne

Background: To assess the early clinical outcomes and toxicities in patients treated with high precision radiation therapy (RT) consisting of helical tomotherapy (HT) or intensity-modulated radiation therapy (IMRT) for anal cancer. Materials and Methods: Since March 2006, 30 patients with stage I-IIIB anal squamous-cell carcinoma were treated curatively by IMRT or HT alone (n = 2) or by concomitant chemotherapy and IMRT or HT (n = 28). Median age was 59 years (range, 36−83 years) and the female/male ratio was 2.3 (21/9). Primary tumor site was anal canal, anal margin, or both in 26, 1, and 3 patients, respectively. Anal tumor, pelvic and inguinal nodes were irradiated with a median dose of 36 Gy using HT, or 5- or 7-field IMRT in 18 and 12 patients, respectively; After a planned gap of 1−2 weeks (median 1 week), a median boost dose of 23.4 Gwas delivered to the tumor and/or involved nodes using 3DRT (n = 24) or HT/IMRT (n = 6). The total delivered dose ranged between 59.4 and 64.8 Gy (median, 59.4 Gy). Concomitant chemotherapy consisted of mitomycin C alone (n = 1), mitomycin C and 5-fluorouracil (n = 17) or capecitabin (n = 10) in 28 patients. Common Terminology Criteria for Adverse Events v3.0 scale was used to score acute and late toxicities. Results: All but one patient, who developed progressive local and distant disease at the end of RT, achieved a complete response. Twelve months following RT, one patient had a recurrence at the primary tumor site, salvaged with brachytherapy. After a median follow-up of 7.5 months (range, 1−35 months), no deaths were observed. The 2-year actuarial locoregional control and probability of disease control without colostomy rates were 82% and 79%, respectively. RT was well tolerated without any unplanned treatment interruptions. Grade 1 or 2 acute adverse events consisted of skin toxicity in 8 and 22 patients, diarrhea in 18 and 3 patients, and cystitis in 9 and 2 patients; respectively. Only one patient developed grade 3 mucosal necrosis at the end of the treatment, requiring diverting colostomy. No difference in terms of acute toxicity was observed between patients treated with HT or IMRT. None of the 22 patients with a follow-up of more than 3 months developed grade 3 or more late toxicity. Conclusions: Our preliminary results suggest that HT or IMRT combined with concomitant chemotherapy for anal cancer is effective, and associated with favorable rates of toxicity compared with historical series. Further follow-up is warranted to assess late toxicity.

Expression of the gene encoding the high-Km glucose transporter 2 by the early postimplantation mouse embryo is essential for neural tube defects associated with diabetic embryopathy.

AIMS/HYPOTHESIS: Excess glucose transport to embryos during diabetic pregnancy causes congenital malformations. The early postimplantation embryo expresses the gene encoding the high-Km GLUT2 (also known as SLC2A2) glucose transporter. The hypothesis tested here is that high-Km glucose transport by GLUT2 causes malformations resulting from maternal hyperglycaemia during diabetic pregnancy. MATERIALS AND METHODS: Glut2 mRNA was assayed by RT-PCR. The Km of embryo glucose transport was determined by measuring 0.5-20 mmol/l 2-deoxy[3H]glucose transport. To test whether the GLUT2 transporter is required for neural tube defects resulting from maternal hyperglycaemia, Glut2+/- mice were crossed and transient hyperglycaemia was induced by glucose injection on day 7.5 of pregnancy. Embryos were recovered on day 10.5, and the incidence of neural tube defects in wild-type, Glut2+/- and Glut2-/- embryos was scored. RESULTS: Early postimplantation embryos expressed Glut2, and expression was unaffected by maternal diabetes. Moreover, glucose transport by these embryos showed Michaelis-Menten kinetics of 16.19 mmol/l, consistent with transport mediated by GLUT2. In pregnancies made hyperglycaemic on day 7.5, neural tube defects were significantly increased in wild-type embryos, but Glut2+/- embryos were partially protected from neural tube defects, and Glut2-/- embryos were completely protected from these defects. The frequency of occurrence of wild-type, Glut2+/- and Glut2-/- embryos suggests that the presence of Glut2 alleles confers a survival advantage in embryos before day 10.5. CONCLUSIONS/INTERPRETATIONS: High-Km glucose transport by the GLUT2 glucose transporter during organogenesis is responsible for the embryopathic effects of maternal diabetes.

Twenty-year trends in the prevalence of Down syndrome and other trisomies in Europe: impact of maternal age and prenatal screening.

This study examines trends and geographical differences in total and live birth prevalence of trisomies 21, 18 and 13 with regard to increasing maternal age and prenatal diagnosis in Europe. Twenty-one population-based EUROCAT registries covering 6.1 million births between 1990 and 2009 participated. Trisomy cases included live births, fetal deaths from 20 weeks gestational age and terminations of pregnancy for fetal anomaly. We present correction to 20 weeks gestational age (ie, correcting early terminations for the probability of fetal survival to 20 weeks) to allow for artefactual screening-related differences in total prevalence. Poisson regression was used. The proportion of births in the population to mothers aged 35+ years in the participating registries increased from 13% in 1990 to 19% in 2009. Total prevalence per 10 000 births was 22.0 (95% CI 21.7-22.4) for trisomy 21, 5.0 (95% CI 4.8-5.1) for trisomy 18 and 2.0 (95% CI 1.9-2.2) for trisomy 13; live birth prevalence was 11.2 (95% CI 10.9-11.5) for trisomy 21, 1.04 (95% CI 0.96-1.12) for trisomy 18 and 0.48 (95% CI 0.43-0.54) for trisomy 13. There was an increase in total and total corrected prevalence of all three trisomies over time, mainly explained by increasing maternal age. Live birth prevalence remained stable over time. For trisomy 21, there was a three-fold variation in live birth prevalence between countries. The rise in maternal age has led to an increase in the number of trisomy-affected pregnancies in Europe. Live birth prevalence has remained stable overall. Differences in prenatal screening and termination between countries lead to wide variation in live birth prevalence.

In situ stromal expression of the urokinase/plasmin system correlates with epithelial dysplasia in colorectal adenomas.

An increase of urokinase-type plasminogen activator (uPA) and a decrease of tissue-type PA (tPA) have been associated with the transition from normal to adenomatous colorectal mucosa. Serial sections from 25 adenomas were used to identify PA-related caseinolytic activities by in situ zymography, blocking selectively uPA or tPA. The distribution of uPA, tPA, and type 1 PA inhibitor mRNAs was investigated by nonradioactive in situ hybridization, and the receptor for uPA was detected by immunostaining. Low- and high-grade epithelial cell dysplasia was mapped histologically. Results show that 23 of 25 adenomas expressed uPA-related lytic activity located predominantly in the periphery whereas tPA-related activity was mainly in central areas of adenomas. In 15 of 25 adenomas, uPA mRNA was expressed in stromal cells clustered in foci that coincided with areas of uPA lytic activity. The probability of finding uPA mRNA-reactive cells was significantly higher in areas with high-grade epithelial dysplasia. uPA receptor was mainly stromal and expressed at the periphery. Type 1 PA inhibitor mRNA cellular expression was diffuse in the stroma, in endothelial cells, and in a subpopulation of alpha-smooth muscle cell actin-reactive cells. These results show that a stromal up-regulation of the uPA/plasmin system is associated with foci of severe dysplasia in a subset of colorectal adenomas.

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

Genome-wide meta-analysis of common variant differences between men and women.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Therapeutic Drug Monitoring of the new targeted anticancer agents imatinib, nilotinib, dasatinib, sunitinib, sorafenib and lapatinib by LC tandem mass spectrometry.

The treatment of some cancer patients has shifted from traditional, non-specific cytotoxic chemotherapy to chronic treatment with molecular targeted therapies. Imatinib mesylate, a selective inhibitor of tyrosine kinases (TKIs) is the most prominent example of this new era and has opened the way to the development of several additional TKIs, including sunitinib, nilotinib, dasatinib, sorafenib and lapatinib, in the treatment of various hematological malignancies and solid tumors. All these agents are characterized by an important inter-individual pharmacokinetic variability, are at risk for drug interactions, and are not devoid of toxicity. Additionally, they are administered for prolonged periods, anticipating the careful monitoring of their plasma exposure via Therapeutic Drug Monitoring (TDM) to be an important component of patients' follow-up. We have developed a liquid chromatography-tandem mass spectrometry method (LC-MS/MS) requiring 100 microL of plasma for the simultaneous determination of the six major TKIs currently in use. Plasma is purified by protein precipitation and the supernatant is diluted in ammonium formate 20 mM (pH 4.0) 1:2. Reverse-phase chromatographic separation of TKIs is obtained using a gradient elution of 20 mM ammonium formate pH 2.2 and acetonitrile containing 1% formic acid, followed by rinsing and re-equilibration to the initial solvent composition up to 20 min. Analyte quantification, using matrix-matched calibration samples, is performed by electro-spray ionization-triple quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. The method was validated according to FDA recommendations, including assessment of extraction yield, matrix effects variability (<9.6%), overall process efficiency (87.1-104.2%), as well as TKIs short- and long-term stability in plasma. The method is precise (inter-day CV%: 1.3-9.4%), accurate (-9.2 to +9.9%) and sensitive (lower limits of quantification comprised between 1 and 10 ng/mL). This is the first broad-range LC-MS/MS assay covering the major currently in-use TKIs. It is an improvement over previous methods in terms of convenience (a single extraction procedure for six major TKIs, reducing significantly the analytical time), sensitivity, selectivity and throughput. It may contribute to filling the current knowledge gaps in the pharmacokinetics/pharmacodynamics relationships of the latest TKIs developed after imatinib and better define their therapeutic ranges in different patient populations in order to evaluate whether a systematic TDM-guided dose adjustment of these anticancer drugs could contribute to minimize the risk of major adverse reactions and to increase the probability of efficient, long lasting, therapeutic response.

Comparison of copper-67- and iodine-125-labeled anti-CEA monoclonal antibody biodistribution in patients with colorectal tumors.

Copper-67 has comparable beta-particle emissions to that of 131I, but it displays more favorable gamma emission characteristics for application in radioimmunotherapy (RIT). This study investigates the potential of 67Cu-labeled monoclonal antibody (MAb) 35 for RIT of colorectal carcinoma. METHODS: Biokinetics of simultaneously injected 67Cu- and 125I-labeled MAb35 were studied in six patients scheduled for surgery of primary colorectal cancer. RESULTS: Whole-body clearance (T 1/2) of 67Cu, estimated from sequential anterior and posterior whole-body scans and corrected for decay of 67Cu, was 41 hr. Serum clearance of 67Cu was faster (27.41 hr) than that of 125I (38.33 hr). Mean tumor uptake of the 67Cu-labeled compound (0.0133% ID/g) exceeded that of 125I (0.0095% ID/g), and tumor-to-blood ratios were higher for 67Cu than for 125I, with averages of 6.07 and 2.41, respectively. The average 67Cu/125I ratio was 1.9 for tumor uptake, 0.7 for blood and 2.6 for tumor-to-blood ratios. Nonspecific liver uptake of 67Cu as calculated from whole-body scans was high in four patients, up to 25% of residual whole-body activity at 48 hr, but did not increase with time. We also observed some nonspecific bowel activity, as well as moderate to high uptake in benign polyps. CONCLUSION: Copper-67-labeled MAb35 is more favorable than its radioiodine-labeled counterpart for RIT of colorectal carcinoma due to higher tumor-to-blood ratios, but the problem of nonspecific liver and bowel uptake must first be overcome. The absolute accumulation of activity in tumor remains low, however, so the probability of cure with this compound alone is questionable. The use of 67Cu as one component of a multimodality adjuvant treatment seems to remain the most appropriate application for RIT.

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.

Lumbar spine texture enhances 10-year fracture probability assessment.

We found that lumbar spine texture analysis using trabecular bone score (TBS) is a risk factor for MOF and a risk factor for death in a retrospective cohort study from a large clinical registry for the province of Manitoba, Canada. INTRODUCTION: FRAX® estimates the 10-year probability of major osteoporotic fracture (MOF) using clinical risk factors and femoral neck bone mineral density (BMD). Trabecular bone score (TBS), derived from texture in the spine dual X-ray absorptiometry (DXA) image, is related to bone microarchitecture and fracture risk independently of BMD. Our objective was to determine whether TBS provides information on MOF probability beyond that provided by the FRAX variables. METHODS: We included 33,352 women aged 40-100 years (mean 63 years) with baseline DXA measurements of lumbar spine TBS and femoral neck BMD. The association between TBS, the FRAX variables, and the risk of MOF or death was examined using an extension of the Poisson regression model. RESULTS: During the mean of 4.7 years, 1,754 women died and 1,872 sustained one or more MOF. For each standard deviation reduction in TBS, there was a 36 % increase in MOF risk (HR 1.36, 95 % CI 1.30-1.42, p < 0.001) and a 32 % increase in death (HR 1.32, 95 % CI 1.26-1.39, p < 0.001). When adjusted for significant clinical risk factors and femoral neck BMD, lumbar spine TBS was still a significant predictor of MOF (HR 1.18, 95 % CI 1.12-1.23) and death (HR 1.20, 95 % CI 1.14-1.26). Models for estimating MOF probability, accounting for competing mortality, showed that low TBS (10th percentile) increased risk by 1.5-1.6-fold compared with high TBS (90th percentile) across a broad range of ages and femoral neck T-scores. CONCLUSIONS: Lumbar spine TBS is able to predict incident MOF independent of FRAX clinical risk factors and femoral neck BMD even after accounting for the increased death hazard.

Mutant p53 promotes epithelial-mesenchymal plasticity and enhances metastasis in mammary carcinomas of WAP-T mice.

To study the postulated mutant p53 (mutp53) "gain of function" effects in mammary tumor development, progression and metastasis, we crossed SV40 transgenic WAP-T mice with mutant p53 transgenic WAP-mutp53 mice. Compared to tumors in monotransgenic WAP-T mice, tumors in bitransgenic WAP-T x WAP-mutp53 mice showed higher tumor grading, enhanced vascularization, and significantly increased metastasis. Bitransgenic tumors revealed a gene signature associated with the oncogenic epithelial-mesenchymal transition pathway (EMT gene signature). In cultures of WAP-T tumor-derived G-2 cancer cells, which are comprised of subpopulations displaying "mesenchymal" and "epithelial" phenotypes, this EMT gene signature was associated with the "mesenchymal" compartment. Furthermore, ectopic expression of mutp53 in G-2 cells sufficed to induce a strong EMT phenotype. In contrast to these in vitro effects, monotransgenic and bitransgenic tumors were phenotypically similar suggesting that in vivo the tumor cell phenotype might be under control of the tumor microenvironment. In support, orthotopic transplantation of G-2 cells as well as of G-2 cells expressing ectopic mutp53 into syngeneic mice resulted in tumors with a predominantly epithelial phenotype, closely similar to that of endogenous primary tumors. We conclude that induction of an EMT gene signature by mutp53 in bitransgenic tumors primarily promotes tumor cell plasticity, that is, the probability of tumor cells to undergo EMT processes under appropriate stimuli, thereby possibly increasing their potential to disseminate and metastasize.

The demographic benefits of belligerence and bravery: defeated group repopulation or victorious group size expansion?

Intraspecific coalitional aggression between groups of individuals is a widespread trait in the animal world. It occurs in invertebrates and vertebrates, and is prevalent in humans. What are the conditions under which coalitional aggression evolves in natural populations? In this article, I develop a mathematical model delineating conditions where natural selection can favor the coevolution of belligerence and bravery between small-scale societies. Belligerence increases an actor's group probability of trying to conquer another group and bravery increase the actors's group probability of defeating an attacked group. The model takes into account two different types of demographic scenarios that may lead to the coevolution of belligerence and bravery. Under the first, the fitness benefits driving the coevolution of belligerence and bravery come through the repopulation of defeated groups by fission of victorious ones. Under the second demographic scenario, the fitness benefits come through a temporary increase in the local carrying capacity of victorious groups, after transfer of resources from defeated groups to victorious ones. The analysis of the model suggests that the selective pressures on belligerence and bravery are stronger when defeated groups can be repopulated by victorious ones. The analysis also suggests that, depending on the shape of the contest success function, costly bravery can evolve in groups of any size.

Can vouchers deliver? An evaluation of subsidies for maternal health care in Cambodia

OBJECTIVE: To evaluate the effect of vouchers for maternity care in public health-care facilities on the utilization of maternal health-care services in Cambodia. METHODS: The study involved data from the 2010 Cambodian Demographic and Health Survey, which covered births between 2005 and 2010. The effect of voucher schemes, first implemented in 2007, on the utilization of maternal health-care services was quantified using a difference-in-differences method that compared changes in utilization in districts with voucher schemes with changes in districts without them. FINDINGS: Overall, voucher schemes were associated with an increase of 10.1 percentage points (pp) in the probability of delivery in a public health-care facility; among women from the poorest 40% of households, the increase was 15.6 pp. Vouchers were responsible for about one fifth of the increase observed in institutional deliveries in districts with schemes. Universal voucher schemes had a larger effect on the probability of delivery in a public facility than schemes targeting the poorest women. Both types of schemes increased the probability of receiving postnatal care, but the increase was significant only for non-poor women. Universal, but not targeted, voucher schemes significantly increased the probability of receiving antenatal care. CONCLUSION: Voucher schemes increased deliveries in health centres and, to a lesser extent, improved antenatal and postnatal care. However, schemes that targeted poorer women did not appear to be efficient since these women were more likely than less poor women to be encouraged to give birth in a public health-care facility, even with universal voucher schemes.