Counterfactuals and causal inference: Methods and principles for social research

"Most quantitative empirical analyses are motivated by the desire to estimate the causal effect of an independent variable on a dependent variable. Although the randomized experiment is the most powerful design for this task, in most social science research done outside of psychology, experimental designs are infeasible. (Winship & Morgan, 1999, p. 659)." This quote from earlier work by Winship and Morgan, which was instrumental in setting the groundwork for their book, captures the essence of our review of Morgan and Winship's book: It is about causality in nonexperimental settings.

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Genetics of cadmium hyperaccumulation variance in Thlaspi caerulescens natural populations from Switzerland

SUMMARY Heavy metal presence in the environment is a serious concern since some of them can be toxic to plants, animals and humans once accumulated along the food chain. Cadmium (Cd) is one of the most toxic heavy metal. It is naturally present in soils at various levels and its concentration can be increased by human activities. Several plants however have naturally developed strategies allowing them to grow on heavy metal enriched soils. One of them consists in the accumulation and sequestration of heavy metals in the above-ground biomass. Some plants present in addition an extreme strategy by which they accumulate a limited number of heavy metals in their shoots in amounts 100 times superior to those expected for a non-accumulating plant in the same conditions. Understanding the genetic basis of the hyperaccumulation trait - particularly for Cd - remains an important challenge which may lead to biotechnological applications in the soil phytoremediation. In this thesis, Thlaspi caerulescens J. & C. Presl (Brassicaceae) was used as a model plant to study the Cd hyperaccumulation trait, owing to its physiological and genetic characteristics. Twenty-four wild populations were sampled in different regions of Switzerland. They were characterized for environmental and soil parameters as well as intrinsic characteristics of plants (i.e. metal concentrations in shoots). They were as well genetically characterized by AFLPs, plastid DNA polymorphism and genes markers (CAPS and microsatellites) mainly developed in this thesis. Some of the investigated genes were putatively linked to the Cd hyperaccumulation trait. Since the study of the Cd hyperaccumulation in the field is important as it allows the identification of patterns of selection, the present work offered a methodology to define the Cd hyperaccumulation capacity of populations from different habitats permitting thus their comparison in the field. We showed that Cd, Zn, Fe and Cu accumulations were linked and that populations with higher Cd hyperaccumulation capacity had higher shoot and reproductive fitness. Using our genetic data, statistical methods (Beaumont & Nichols's procedure, partial Mantel tests) were applied to identify genomic signatures of natural selection related to the Cd hyperaccumulation capacity. A significant genetic difference between populations related to their Cd hyperaccumulation capacity was revealed based on somè specific markers (AFLP and candidate genes). Polymorphism at the gene encoding IRTl (Iron-transporter also participating to the transport of Zn) was suggested as explaining part of the variation in Cd hyperaccumulation capacity of populations supporting previous physiological investigations. RÉSUMÉ La présence de métaux lourds dans l'environnement est un phénomène préoccupant. En effet, certains métaux lourds - comme le cadmium (Cd) -sont toxiques pour les plantes, les animaux et enfin, accumulés le long de la chaîne alimentaire, pour les hommes. Le Cd est naturellement présent dans le sol et sa concentration peut être accrue par différentes activités humaines. Certaines plantes ont cependant développé des stratégies leur permettant de pousser sur des sols contaminés en métaux lourds. Parmi elles, certaines accumulent et séquestrent les métaux lourds dans leurs parties aériennes. D`autres présentent une stratégie encore plus extrême. Elles accumulent un nombre limité de métaux lourds en quantités 100 fois supérieures à celles attendues pour des espèces non-accumulatrices sous de mêmes conditions. La compréhension des bases génétiques de l'hyperaccumulation -particulièrement celle du Cd - représente un défi important avec des applications concrètes en biotechnologies, tout particulièrement dans le but appliqué de la phytoremediation des sols contaminés. Dans cette thèse, Thlaspi caerulescens J. & C. Presl (Brassicaceae) a été utilisé comme modèle pour l'étude de l'hyperaccumulation du Cd de par ses caractéristiques physiologiques et génétiques. Vingt-quatre populations naturelles ont été échantillonnées en Suisse et pour chacune d'elles les paramètres environnementaux, pédologique et les caractéristiques intrinsèques aux plantes (concentrations en métaux lourds) ont été déterminés. Les populations ont été caractérisées génétiquement par des AFLP, des marqueurs chloroplastiques et des marqueurs de gènes spécifiques, particulièrement ceux potentiellement liés à l'hyperaccumulation du Cd (CAPS et microsatellites). La plupart ont été développés au cours de cette thèse. L'étude de l'hyperaccumulation du Cd en conditions naturelles est importante car elle permet d'identifier la marque, éventuelle de sélection naturelle. Ce travail offre ainsi une méthodologie pour définir et comparer la capacité des populations à hyperaccumuler le Cd dans différents habitats. Nous avons montré que les accumulations du Cd, Zn, Fe et Cu sont liées et que les populations ayant une grande capacité d'hyperaccumuler le Cd ont également une meilleure fitness végétative et reproductive. Des méthodes statistiques (l'approche de Beaumont & Nichols, tests de Martel partiels) ont été utilisées sur les données génétiques pour identifier la signature génomique de la sélection naturelle liée à la capacité d'hyperaccumuler le Cd. Une différenciation génétique des populations liée à leur capacité d'hyperaccumuler le Cd a été mise en évidence sur certains marqueurs spécifiques. En accord avec les études physiologiques connues, le polymorphisme au gène codant IRT1 (un transporteur de Fe impliqué dans le transport du Zn) pourrait expliquer une partie de la variance de la capacité des populations à hyperaccumuler le Cd.

Abnormalities of fuel utilization as predisposing to the development of obesity in humans.

A number of recent investigations in man have demonstrated that a low ratio of fat to carbohydrate oxidation (i.e., a high respiratory quotient or RQ) was associated with actual and/or subsequent body weight gain in obese non-diabetic Pima Indians, in American men of various ages and in post-obese European women investigated shortly after the cessation of a hypocaloric diet. It is well known that numerous exogenous and endogenous factors influence the RQ at rest such as: the level of feeding (positive vs. negative energy balance), the composition of food eaten (high vs. low carbohydrate), the size of the glycogen stores, the amount of adipose tissue as well as genetic factors. It should be stressed that some nutritional situations can co-exist during which a low ratio of fat to carbohydrate is observed (i.e., a high RQ) despite weight loss. Furthermore, in most studies mentioned above, the low fat to carbohydrate oxidation ratio explains less than 10% of the variance in weight gain, suggesting that numerous additional factors also play a substantial role in the onset of weight gain. It is concluded that: 1) a low fat to carbohydrate oxidation ratio or an abnormal fat oxidation is difficult to define quantitatively since it is largely influenced by the energy level and the composition of the diet.(ABSTRACT TRUNCATED AT 250 WORDS)

Breeding system and genetic variance in the monogamous, semi-social shrew, Crocidura russula

The population-genetic consequences of monogamy and male philopatry (a rare breeding system in mammals) were investigated using microsatellite markers in the semisocial and anthropophilic shrew Crocidura russula. A hierarchical sampling design over a 16-km geographical transect revealed a large genetic diversity (h = 0.813) with significant differentiation among subpopulations (F-ST = 5-6%), which suggests an exchange of 4.4 migrants per generation. Demic effective-size estimates were very high, due both to this limited gene inflow and to the inner structure of subpopulations. These were made of 13-20 smaller units (breeding groups), comprising an estimate of four breeding pairs each. Members of the same breeding groups displayed significant coancestries (F-LS = 9-10%), which was essentially due to strong male kinship: syntopic males were on average related at the half-sib level. Female dispersal among breeding groups was not complete (similar to 39%), and insufficient to prevent inbreeding. From our results, the breeding strategy of C. russula seems less efficient than classical mammalian systems (polygyny and male dispersal) in disentangling coancestry from inbreeding, but more so in retaining genetic variance.

Sample size and statistical power in the hierarchical analysis of variance: applications in morphometry of the nervous system.

Analysis of variance is commonly used in morphometry in order to ascertain differences in parameters between several populations. Failure to detect significant differences between populations (type II error) may be due to suboptimal sampling and lead to erroneous conclusions; the concept of statistical power allows one to avoid such failures by means of an adequate sampling. Several examples are given in the morphometry of the nervous system, showing the use of the power of a hierarchical analysis of variance test for the choice of appropriate sample and subsample sizes. In the first case chosen, neuronal densities in the human visual cortex, we find the number of observations to be of little effect. For dendritic spine densities in the visual cortex of mice and humans, the effect is somewhat larger. A substantial effect is shown in our last example, dendritic segmental lengths in monkey lateral geniculate nucleus. It is in the nature of the hierarchical model that sample size is always more important than subsample size. The relative weight to be attributed to subsample size thus depends on the relative magnitude of the between observations variance compared to the between individuals variance.

Effects of common origin and common rearing environment on variance in ectoparasite load and phenotype of nestling Alpine swifts

Knowledge of the quantitative genetics of resistance to parasitism is key to appraise host evolutionary responses to parasite selection. Here, we studied effects of common origin (i.e. genetic and pre-hatching parental effects) and common rearing environment (i.e. post-hatching parental effects and other environment effects) on variance in ectoparasite load in nestling Alpine swifts (Apus melba). This colonial bird is intensely parasitized by blood sucking louse-flies that impair nestling development and survival. By cross-fostering half of the hatchlings between pairs of nests, we show strong significant effect of common rearing environment on variance (90.7% in 2002 and 90.9% in 2003) in the number of louse-flies per nestling and no significant effect of common origin on variance in the number of louse-flies per nestling. In contrast, significant effects of common origin were found for all the nestling morphological traits (i.e. body mass, wing length, tail length, fork length and sternum length) under investigation. Hence, our study suggests that genetic and pre-hatching parental effects play little role in the distribution of parasites among nestling Alpine swifts, and thus that nestlings have only limited scope for evolutionary responses against parasites. Our results highlight the need to take into consideration environmental factors, including the evolution of post-hatching parental effects such as nest sanitation, in our understanding of host-parasite relationships.

Natural selection on fecundity variance in subdivided populations: kin selection meets bet hedging.

In a series of seminal articles in 1974, 1975, and 1977, J. H. Gillespie challenged the notion that the "fittest" individuals are those that produce on average the highest number of offspring. He showed that in small populations, the variance in fecundity can determine fitness as much as mean fecundity. One likely reason why Gillespie's concept of within-generation bet hedging has been largely ignored is the general consensus that natural populations are of large size. As a consequence, essentially no work has investigated the role of the fecundity variance on the evolutionary stable state of life-history strategies. While typically large, natural populations also tend to be subdivided in local demes connected by migration. Here, we integrate Gillespie's measure of selection for within-generation bet hedging into the inclusive fitness and game theoretic measure of selection for structured populations. The resulting framework demonstrates that selection against high variance in offspring number is a potent force in large, but structured populations. More generally, the results highlight that variance in offspring number will directly affect various life-history strategies, especially those involving kin interaction. The selective pressures on three key traits are directly investigated here, namely within-generation bet hedging, helping behaviors, and the evolutionary stable dispersal rate. The evolutionary dynamics of all three traits are markedly affected by variance in offspring number, although to a different extent and under different demographic conditions.

In postmenopausal women with osteoporosis, denosumab significantly improved trabecular bone score (TBS), an index of trabecular microarchitecture

Background/Purpose: The trabecular bone score (TBS), a novel graylevel texture index determined from lumbar spine DXA scans, correlates with 3D parameters of trabecular bone microarchitecture known to predict fracture. TBS may enhance the identification of patients at increased risk for vertebral fracture independently of bone mineral density (BMD) (Boutroy JBMR 2010; Hans JBMR 2011). Denosumab treatment for 36 months decreased bone turnover, increased BMD, and reduced new vertebral fractures in postmenopausal women with osteoporosis (Cummings NEJM 2009). We explored the effect of denosumab on TBS over 36 months and evaluated the association between TBS and lumbar spine BMD in women who had DXA scans obtained from eligible scanners for TBS evaluation in FREEDOM. Methods: FREEDOM was a 3-year, randomized, double-blind trial that enrolled postmenopausal women with a lumbar spine or total hip DXA T-score __2.5, but not __4.0 at both sites. Women received placebo or 60 mg denosumab every 6 months. A subset of women in FREEDOM participated in a DXA substudy where lumbar spine DXA scans were obtained at baseline and months 1, 6, 12, 24, and 36. We retrospectively applied, in a blinded-to-treatment manner, a novel software program (TBS iNsightR v1.9, Med-Imaps, Pessac, France) to the standard lumbar spine DXA scans obtained in these women to determine their TBS indices at baseline and months 12, 24, and 36. From previous studies, a TBS _1.35 is considered as normal microarchitecture, a TBS between 1.35 and _1.20 as partially deteriorated, and 1.20 reflects degraded microarchitecture. Results: There were 285 women (128 placebo, 157 denosumab) with a TBS value at baseline and _1 post-baseline visit. Their mean age was 73, their mean lumbar spine BMD T-score was _2.79, and their mean lumbar spine TBS was 1.20. In addition to the robust gains in DXA lumbar spine BMD observed with denosumab (9.8% at month 36), there were consistent, progressive, and significant increases in TBS compared with placebo and baseline (Table & Figure). BMD explained a very small fraction of the variance in TBS at baseline (r2_0.07). In addition, the variance in the TBS change was largely unrelated to BMD change, whether expressed in absolute or percentage changes, regardless of treatment, throughout the study (all r2_0.06); indicating that TBS provides distinct information, independently of BMD. Conclusion: In postmenopausal women with osteoporosis, denosumab significantly improved TBS, an index of lumbar spine trabecular microarchitecture, independently of BMD.

Stepwise evaluation of unexplained syncope in a large ambulatory population.

BACKGROUND: Up to 60% of syncopal episodes remain unexplained. We report the results of a standardized, stepwise evaluation of patients referred to an ambulatory clinic for unexplained syncope. METHODS AND RESULTS: We studied 939 consecutive patients referred for unexplained syncope, who underwent a standardized evaluation, including history, physical examination, electrocardiogram, head-up tilt testing (HUTT), carotid sinus massage (CSM) and hyperventilation testing (HYV). Echocardiogram and stress test were performed when underlying heart disease was initially suspected. Electrophysiological study (EPS) and implantable loop recorder (ILR) were used only in patients with underlying structural heart disease or major unexplained syncope. We identified a cause of syncope in 66% of patients, including 27% vasovagal, 14% psychogenic, 6% arrhythmias, and 6% hypotension. Noninvasive testing identified 92% and invasive testing an additional 8% of the causes. HUTT yielded 38%, CSM 28%, HYV 49%, EPS 22%, and ILR 56% of diagnoses. On average, patients with arrhythmic causes were older, had a lower functional capacity, longer P-wave duration, and presented with fewer prodromes than patients with vasovagal or psychogenic syncope. CONCLUSIONS: A standardized stepwise evaluation emphasizing noninvasive tests yielded 2/3 of causes in patients referred to an ambulatory clinic for unexplained syncope. Neurally mediated and psychogenic mechanisms were behind >50% of episodes, while cardiac arrhythmias were uncommon. Sudden syncope, particularly in older patients with functional limitations or a prolonged P-wave, suggests an arrhythmic cause.

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.

Patterns of oxygen consumption during establishment of cephalocaudal polarity in the early chick embryo.

Oxygen uptake was studied during the establishment of cephalocaudal polarity in the very early chick embryo, i.e., 10 hr before (stage VI) and at laying (stage X). Oxygen fluxes in minute regions of the intact blastoderms were measured in vitro by scanning microspectrophotometry in the presence or absence of glucose. The oxygen consumption of the whole blastoderm remained constant (6 nmol O2 X hr-1) throughout the period studied, although the number of cells increased more than twofold. The regional oxygen fluxes varied from 0.41 to 1.13 nmol O2 X hr-1 X mm-2 at stage VI and from 0.42 to 0.70 nmol O2 X hr-1 X mm-2 at stage X. At stage VI, the oxygen flux in the center of the blastoderm was significantly higher than that in its periphery. This pattern remained evident when the values were corrected for cell number or for cytoplasmic volume. At stage X, there was a tendency for the oxygen fluxes to decrease from the posterior to the anterior regions of the area pellucida. Thus the pattern of oxidative metabolism in the late uterine embryos seems to change from radial to bilateral. This change of symmetry probably reflects the process of formation of the embryonic axis. In addition, the fact that the oxygen uptake was similar in the presence or absence of glucose suggests that early chick embryos metabolize essentially intracellular stores.

Effects of infused amino acids on glucose production and utilization in healthy human subjects.

Amino acids have been reported to increase endogenous glucose production in normal human subjects during hyperinsulinemia: however, controversy exists as to whether insulin-mediated glucose disposal is inhibited under these conditions. The effect of an amino acid infusion on glucose oxidation rate has so far not been determined. Substrate oxidation rates, endogenous glucose production, and [13C]glucose synthesis from [13C]bicarbonate were measured in six normal human subjects during sequential infusions of exogenous glucose and exogenous glucose with (n = 5) or without (n = 5) exogenous amino acids. Amino acids increased endogenous glucose production by 84% and [13C]glucose synthesis by 235%. Glucose oxidation estimated from indirect calorimetry decreased slightly after amino acids, but glucose oxidation estimated from [13C]glucose-13CO2 data was increased by 14%. It is concluded that gluconeogenesis is the major pathway of amino acid degradation. During amino acid administration, indirect calorimetry underestimates the true rate of glucose oxidation, whereas glucose oxidation calculated from the 13C enrichment of expired CO2 during [U-13C]glucose infusion does not. A slight stimulation of glucose oxidation during amino acid infusion, concomitant with an increased plasma insulin concentration, indicates that amino acids do not inhibit glucose oxidation.