Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

Effects of brood size manipulation and common origin on phenotype and telomere length in nestling collared flycatchers.

BACKGROUND: Evidence is accumulating that telomere length is a good predictor of life expectancy, especially early in life, thus calling for determining the factors that affect telomere length at this stage. Here, we investigated the relative influence of early growth conditions and origin (genetics and early maternal effects) on telomere length of collared flycatchers (Ficedula albicollis) at fledging. We experimentally transferred hatchlings among brood triplets to create reduced, control (i.e. unchanged final nestling number) and enlarged broods. RESULTS: Although our treatment significantly affected body mass at fledging, we found no evidence that increased sibling competition affected nestling tarsus length and telomere length. However, mixed models showed that brood triplets explained a significant part of the variance in body mass (18%) and telomere length (19%), but not tarsus length (13%), emphasizing that unmanipulated early environmental factors influenced telomere length. These models also revealed low, but significant, heritability of telomere length (h(2) = 0.09). For comparison, the heritability of nestling body mass and tarsus length was 0.36 and 0.39, respectively, which was in the range of previously published estimates for those two traits in this species. CONCLUSION: Those findings in a wild bird population demonstrate that telomere length at the end of the growth period is weakly, but significantly, determined by genetic and/or maternal factors taking place before hatching. However, we found no evidence that the brood size manipulation experiment, and by extension the early growth conditions, influenced nestling telomere length. The weak heritability of telomere length suggests a close association with fitness in natural populations.

Fast subplasma membrane Ca2+ transients control exo-endocytosis of synaptic-like microvesicles in astrocytes

Astrocytes are the most abundant glial cell type in the brain. Although not apposite for long-range rapid electrical communication, astrocytes share with neurons the capacity of chemical signaling via Ca(2+)-dependent transmitter exocytosis. Despite this recent finding, little is known about the specific properties of regulated secretion and vesicle recycling in astrocytes. Important differences may exist with the neuronal exocytosis, starting from the fact that stimulus-secretion coupling in astrocytes is voltage independent, mediated by G-protein-coupled receptors and the release of Ca(2+) from internal stores. Elucidating the spatiotemporal properties of astrocytic exo-endocytosis is, therefore, of primary importance for understanding the mode of communication of these cells and their role in brain signaling. We here take advantage of fluorescent tools recently developed for studying recycling of glutamatergic vesicles at synapses (Voglmaier et al., 2006; Balaji and Ryan, 2007); we combine epifluorescence and total internal reflection fluorescence imaging to investigate with unprecedented temporal and spatial resolution, the stimulus-secretion coupling underlying exo-endocytosis of glutamatergic synaptic-like microvesicles (SLMVs) in astrocytes. Our main findings indicate that (1) exo-endocytosis in astrocytes proceeds with a time course on the millisecond time scale (tau(exocytosis) = 0.24 +/- 0.017 s; tau(endocytosis) = 0.26 +/- 0.03 s) and (2) exocytosis is controlled by local Ca(2+) microdomains. We identified submicrometer cytosolic compartments delimited by endoplasmic reticulum tubuli reaching beneath the plasma membrane and containing SLMVs at which fast (time-to-peak, approximately 50 ms) Ca(2+) events occurred in precise spatial-temporal correlation with exocytic fusion events. Overall, the above characteristics of transmitter exocytosis from astrocytes support a role of this process in fast synaptic modulation.

Serum uric acid and adiposity: deciphering causality using a bidirectional Mendelian randomization approach.

Although the relationship between serum uric acid (SUA) and adiposity is well established, the direction of the causality is still unclear in the presence of conflicting evidences. We used a bidirectional Mendelian randomization approach to explore the nature and direction of causality between SUA and adiposity in a population-based study of Caucasians aged 35 to 75 years. We used, as instrumental variables, rs6855911 within the SUA gene SLC2A9 in one direction, and combinations of SNPs within the adiposity genes FTO, MC4R and TMEM18 in the other direction. Adiposity markers included weight, body mass index, waist circumference and fat mass. We applied a two-stage least squares regression: a regression of SUA/adiposity markers on our instruments in the first stage and a regression of the response of interest on the fitted values from the first stage regression in the second stage. SUA explained by the SLC2A9 instrument was not associated to fat mass (regression coefficient [95% confidence interval]: 0.05 [-0.10, 0.19] for fat mass) contrasting with the ordinary least square estimate (0.37 [0.34, 0.40]). By contrast, fat mass explained by genetic variants of the FTO, MC4R and TMEM18 genes was positively and significantly associated to SUA (0.31 [0.01, 0.62]), similar to the ordinary least square estimate (0.27 [0.25, 0.29]). Results were similar for the other adiposity markers. Using a bidirectional Mendelian randomization approach in adult Caucasians, our findings suggest that elevated SUA is a consequence rather than a cause of adiposity.

Targeting burn prevention in the paediatric population : a prospective study of children's burns in the Lausanne area.

QUESTION UNDER STUDY: Domestic accidents are an important problem in paediatric medicine. This study was designed to gain a better understanding of burn mechanisms and target prevention. METHODS: Children treated for burn lesions in the Department of Paediatric Surgery between August 2004 and August 2005 were included in this prospective study. The burn mechanisms, the children's ages and the circumstances in which children were burned as well as their home environment variables were analyzed. RESULTS: The current study included eighty-nine patients, aged between 2 months and 15 years. Seventy-eight percent were less than 5 years old. More than half were boys. Hot liquid scalding was the most frequent mechanism. There does not seem to be an increased risk in the immigrant population or in low economic status families. In most cases, an adult person was present at time of injury. CONCLUSIONS: If we were to describe the highest "at risk" candidate for a burn in our region, it would be a boy aged 15 months to 5 years who is burned by a cup of hot liquid on his hand, at home, around mealtime, in the presence of one or both parents. Reduced attention in the safe domestic setting is probably responsible.

Biological embedding of early life exposures and disease risk in humans: a role for DNA methylation.

BACKGROUND: Following wider acceptance of "the thrifty phenotype" hypothesis and the convincing evidence that early life exposures can influence adult health even decades after the exposure, much interest has been placed on the mechanisms through which early life exposures become biologically embedded. METHODS: In this review, we summarize the current literature regarding biological embedding of early life experiences. To this end we conducted a literature search to identify studies investigating early life exposures in relation to DNA methylation changes. In addition, we summarize the challenges faced in investigations of epigenetic effects, stemming from the peculiarities of this emergent and complex field. A proper systematic review and meta-analyses were not feasible given the nature of the evidence. RESULTS: We identified 7 studies on early life socioeconomic circumstances, 10 studies on childhood obesity, and 6 studies on early life nutrition all relating to DNA methylation changes that met the stipulated inclusion criteria. The pool of evidence gathered, albeit small, favours a role of epigenetics and DNA methylation in biological embedding, but replication of findings, multiple comparison corrections, publication bias, and causality are concerns remaining to be addressed in future investigations. CONCLUSIONS: Based on these results, we hypothesize that epigenetics, in particular DNA methylation, is a plausible mechanism through which early life exposures are biologically embedded. This review describes the current status of the field and acts as a stepping stone for future, better designed investigations on how early life exposures might become biologically embedded through epigenetic effects. This article is protected by copyright. All rights reserved.

Evolution of life history traits in a natural population of the shrew Crocidura russula

Abstract Life history traits encompass all the decisions concerning fitness an individual is faced with during his life. The study of these traits is crucial to understand the factors shaping the biology of living organisms. Up until now, most of the information on the evolution of life history traits comes from laboratory studies. While these studies are interesting to test the effect of specific parameters, their conclusions are difficult to extrapolate to natural populations. Investigating the evolution of life history traits in natural populations is of great interest. This may be tricky because it requires information on reproduction, survival and morphology of individuals. Mark-recapture methods allow most of this information to be obtained. However, when direct observations of a species are not possible due to its ecology, indirect methods must be used to infer lifetime reproductive success. In this case, molecular markers are particularly helpful in assessing the genetic relationships between individuals and allow the construction of a pedigree. This thesis focuses on a natural population of a small insectivorous mammal, the greater white-toothed shrew, Crocidura russula. Because of its hidden lifestyle, the two complementary techniques mentioned above were combined to gather information on this population. The data were used to explore diverse aspects of evolutionary biology. We demonstrated that the high genetic variance displayed by the species was not maintained by its mating system because this shrew was less monogamous than previously thought. The large genetic diversity was most likely promoted by gene flow from the neighborhood. Dispersal was thus a central topic in this thesis. We showed that dispersal was not driven by inbreeding avoidance. In addition, we did not find any inbreeding depression in the population. Dispersal was promoted by a high number of vacant territories in the population for both sexes, meaning that territory acquisition played an important role in driving dispersal. Moreover, dispersal propensity was shown to have a genetic basis and, once achieved, to have no effect on individual fitness. Body mass was found to be a life history trait strongly influenced by sexual and viability selection in both sexes. Larger individuals had higher access to reproduction through territory acquisition and defense than lighter ones. By contrast, intermediate size individuals were favored by viability selection presumably because of ecological constraints and metabolic costs. Finally, we demonstrated that the majority of the life history traits in our shrew population has the potential to evolve because they maintained substantial amounts of additive genetic variance. Nonetheless, life history traits had no significant heritability due to their high level of nonadditive or environmental variance. Résumé Les traits d'histoire de vie comprennent toutes les décisions auxquelles un individu est confronté au cours de sa vie et qui concernent sa valeur adaptative. L'étude de ces traits est cruciale pour comprendre les facteurs qui façonnent la biologie des êtres vivants. Jusqu'à ce jour, la majorité des informations sur l'évolution des traits d'histoire de vie provient d'études réalisées en laboratoire. Alors que ces études sont intéressantes pour tester l'effet de paramètres spécifiques, leurs conclusions sont difficilement extrapolables aux populations naturelles. Il est particulièrement intéressant d'étudier l'évolution des traits d'histoire de vie dans des populations naturelles. Toutefois, ces études peuvent se révéler difficiles parce qu'elles requièrent des informations sur la reproduction, la survie et la morphologie des individus. Des méthodes de marquage-recapture permettent d'obtenir ces informations. Cependant, lorsque l'écologie de l'espèce rend les obervations directes impossibles, des méthodes indirectes doivent être utilisées pour obtenir le succès reproducteur des individus. Dans ce cas, les marqueurs moléculaires sont particulièrement utiles pour évaluer les relations génétiques entre individus et permettre la construction d'un pedigree. Cette thèse porte sur une population naturelle d'un petit mammifère insectivore, la musaraigne musette, Crocidura russula. Parce que cette espèce présente un mode de vie souterrain, les deux techniques complémentaires mentionnées ci-dessus ont été combinées pour acquérir les informations nécessaires. Les données ont été utilisées pour explorer divers aspects de biologie evolutive. Nous avons montré que la grande quantité de variance génétique trouvée chez cette espèce n'est pas maintenue par son système d'appariement. Celle-ci s'est en effet avérée être moins monogame que ce qui était admis jusqu'ici. Sa grande diversité génétique est plutôt entretenue par le flux de gènes provenant du voisinage. La dispersion a donc été un sujet phare dans cette thèse. Nous avons montré qu'elle n'est pas provoquée par un évitement de la consanguinité et nous n'avons pas trouvé de dépression de consanguité dans notre population. L'acquisition d'un territoire joue par contre un rôle important dans la dispersion. En outre, la dispersion possède une base génétique chez cette espèce. De plus, une fois qu'ils ont dispersé, les individus n'ont pas une valeur adaptative differente d'individus philopatriques. Le poids s'est avéré être un trait d'histoire de vie fortement influencé par la sélection sexuelle et de viabilité chez les deux sexes. Les gros individus ont accès à la reproduction parce qu'ils acquièrent et défendent un territoire plus facilement que les plus légers. Au contraire, les individus de taille intermédiaire sont favorisés par la sélection de viabilité, certainement à cause de contraintes écologiques et de coûts métaboliques. Finalement, nous avons montré que la majorité des traits d'histoire de vie dans notre population a le potentiel d'évoluer parce qu'elle maintient des quantités considérables de variance génétique additive. Néanmoins, l'héritabilité de ces traits d'histoire de vie n'est pas significative à cause de la grande quantité de variance non-additive ou environmentale associée à ces traits.

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

The spread of incompatibility-inducing parasites in sub-divided host populations.

BACKGROUND:Maternally transmitted symbionts have evolved a variety of ways to promote their spread through host populations. One strategy is to hamper the reproduction of uninfected females by a mechanism called cytoplasmic incompatibility (CI). CI occurs in crosses between infected males and uninfected females and leads to partial to near-complete infertility. CI-infections are under positive frequency-dependent selection and require genetic drift to overcome the range of low frequencies where they are counter-selected. Given the importance of drift, population sub-division would be expected to facilitate the spread of CI. Nevertheless, a previous model concluded that variance in infection between competing groups of breeding individuals impedes the spread of CI.RESULTS:In this paper we derive a model on the spread of CI-infections in populations composed of demes linked by restricted migration. Our model shows that population sub-division facilitates the invasion of CI. While host philopatry (low migration) favours the spread of infection, deme size has a non-monotonous effect, with CI-invasion being most likely at intermediate deme size. Individual-based simulations confirm these predictions and show that high levels of local drift speed up invasion but prevent high levels of prevalence across the entire population. Additional simulations with sex-specific migration rates further show that low migration rates of both sexes are required to facilitate the spread of CI.CONCLUSION:Our analyses show that population structure facilitates the invasion of CI-infections. Since some level of sub-division is likely to occur in most natural populations, our results help to explain the high incidence of CI-infections across species of arthropods. Furthermore, our work has important implications for the use of CI-systems in order to genetically modify natural populations of disease vectors.

Effective exposure to solar UV in building workers: influence of local and individual factors.

Excessive exposure to solar UV light is the main cause of skin cancers in humans. UV exposure depends on environmental as well as individual factors related to activity. Although outdoor occupational activities contribute significantly to the individual dose received, data on effective exposure are scarce and limited to a few occupations. A study was undertaken in order to assess effective short-term exposure among building workers and characterize the influence of individual and local factors on exposure. The effective exposure of construction workers in a mountainous area in the southern part of Switzerland was investigated through short-term dosimetry (97 dosimeters). Three altitudes, of about 500, 1500 and 2500 m were considered. Individual measurements over 20 working periods were performed using Spore film dosimeters on five body locations. The postural activity of workers was concomitantly recorded and static UV measurements were also performed. Effective exposure among building workers was high and exceeded occupational recommendations, for all individuals for at least one body location. The mean daily UV dose in plain was 11.9 SED (0.0-31.3 SED), in middle mountain 21.4 SED (6.6-46.8 SED) and in high mountain 28.6 SED (0.0-91.1 SED). Measured doses between workers and anatomical locations exhibited a high variability, stressing the role of local exposure conditions and individual factors. Short-term effective exposure ranged between 0 and 200% of ambient irradiation, indicating the occurrence of intense, subacute exposures. A predictive irradiation model was developed to investigate the role of individual factors. Posture and orientation were found to account for at least 38% of the total variance of relative individual exposure, and were also found to account more than altitude on the total variance of effective daily exposures. Targeted sensitization actions through professional information channels and specific prevention messages are recommended. Altitude outdoor workers should also benefit from preventive medical examination.

Nutritional status using the Mini Nutritional Assessment questionnaire and its relationship with bone quality in a population of institutionalized elderly women.

Malnutrition, a risk factor for osteoporotic fractures, is frequent in elderly people and, is underdiagnosed and undertreated. There are only few studies on the nutritional status of elderly people in Europe. The Mini Nutritional Assessment (MNA) is a non invasive and validated questionnaire to evaluate nutritional status in elderly people, classified in three groups: 1 degree score < 17: malnourished, 2 degrees score >17 and < 24: at risk of malnutrition, 3 degrees score >24: well-nourished, with a maximum of 30 points. Quantitative ultrasound of bone (QUS) is a method for assessing quality of bone which can be easily performed in nursing homes. Therefore, these two tests allowed to study the relationships between nutritional status and ultrasonic parameters of bone in 78 institutionalized women aged 86 +/- 6 years, living in 11 nursing homes around Lausanne (Switzerland). All were assessed by the MNA, had a measurement of the tricipital skin fold and of the grip strength. Functional status was evaluated by the scale "Activity of Daily Living" (ADL), and serum albumin level was measured when permitted. All had QUS of the calcaneus (with an Achilles, GE Lunar). The measured parameters are the Broadband Ultrasound Attenuation (BUA), attenuation of a band of ultrasonic frequencies through the medium, expressed in dB/MHz, and the Speed of Sound (SOS), speed of the ultrasounds through the medium, expressed in m/s. A third parameter, the stiffness index (SI), expressed as a percentage of the values obtained by the manufacturer in a young population and derived from BUA and SOS, was calculated automatically : SI = (0.67xBUA) + (0.28xSOS) - 420, expressed in percent compared to a young adult population (%YA). Fifteen percent of the women were undernourished and 58% were at risk of malnutrition. As expected, compared with the well-nourished minority, undernourished subjects had significant lower body mass index (BMI), tricipital skin fold (TSF), ADL score and albumin level (p < 0,01). The subjects "at risk of malnutrition" had significant lower BMI, ADL score (p < 0.01), tricipital skin fold and serum albumin (p < 0.05). Ultrasound parameters were low independently of the nutritional status. MNA score correlated significantly with tricipital skin fold (r = 0.508, p < 0.01), ADL (r = 0.538, p < 0.01) and albumin serum level (r = 0.409, p = 0.01). There was a trend for a correlation between the MNA and the ultrasound parameter BUA (r = 0.207, p = 0.07), whereas no correlation was found with SOS and SI. A multivariate analysis showed that tricipital skin fold and ADL explained 61% of the variance of the MNA. In conclusion, using simple and non invasive methods, this study showed that malnutrition and osteoporosis are frequent in institutionalized elderly persons in our country, and the ultrasound parameters are influenced by many others factors in addition to nutrition, especially at this age and in elderly residents of nursing homes.

Dynamic responses of oxygen uptake at the onset and end of moderate and heavy exercise in trained subjects.

Inconsistencies about dynamic asymmetry between the on- and off-transient responses in VO2 are found in the literature. Therefore the purpose of this study was to examine VO2 on- and off-transients during moderate- and heavy-intensity cycling exercise in trained subjects. Ten men underwent an initial incremental test for the estimation of ventilatory threshold (VT) and, on different days, two bouts of square-wave exercise at moderate (<VT) and heavy (>VT) intensities. VO2 kinetics in exercise and recovery were better described by a single exponential model (<VT), or by a double exponential with two time delays (>VT). For moderate exercise, we found a symmetry of VO2 kinetics between the on- and off-transients (i.e., fundamental component), consistent with a system manifesting linear control dynamics. For heavy exercise, a slow component superimposed on the fundamental phase was expressed in both the exercise and recovery, with similar parameter estimates. But the on-transient values of the time constant were appreciably faster than the associated off-transient, and independent of the work rate imposed (<VT and >VT). Our results do not support a dynamically linear system model of VO2 during cycling exercise in the heavy-intensity domain.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Psychometric evaluation of the French version of the questionnaire attitudes towards morphine use; a cross-sectional study in Valais, Switzerland.

BACKGROUND: In Switzerland, nurses are allowed to prescribe and administer morphine in emergency situations without a doctor. Still, nurses and other health professionals are often reluctant to prescribe and administer morphine for pain management in patients. No valid French-speaking instrument is available in Switzerland to assess the attitudes of nurses and other health professionals towards the prescription and administration of morphine. In this study, we evaluated the psychometric properties of the French version of the questionnaire "Attitudes towards morphine use". METHODS: The instrument was derived from an Italian version. Forward and back translations of the questionnaire were performed. Item analysis and construct validity were assessed between April and December 2010 in a cross sectional study including five Swiss hospitals in a sample of 588 health professionals (533 nurses, mean age 38.3 ± 10.2 years). Thirty subjects participated in test-retest reliability. RESULTS: The time to complete the instrument ranged between 12 and 15 minutes and neither floor nor ceiling effect were found. The initial 24-item instrument showed an intraclass correlation (ICC) of 0.69 (95% CI: 0.64 to 0.73, P < 0.001), and a Cronbach's α of 0.700. Factor analysis led to a six-component solution explaining 52.4% of the total variance. After excluding five items, the shortened version showed an ICC of 0.74 (95% CI, 0.70 to 0.77, P < 0.001) and a Cronbach's α of 0.741. Factor analysis led to a five-component solution explaining 54.3% of the total variance. The five components were named "risk of addiction/dependence"; "operational reasons for not using morphine"; "risk of escalation"; "other (non-dependence) risks" and "external (non-operational) reasons". In test-retest, the shortened instrument showed an ICC of 0.797 (95% CI, 0.630 to 0.911, P < 0.001) and a Cronbach's α of 0.797. CONCLUSIONS: The 19-item shortened instrument assessing attitudes towards the prescription and administration of morphine showed adequate content and construct validity.