A quantitative test of Hamilton's rule for the evolution of altruism.

The evolution of altruism is a fundamental and enduring puzzle in biology. In a seminal paper Hamilton showed that altruism can be selected for when rb - c > 0, where c is the fitness cost to the altruist, b is the fitness benefit to the beneficiary, and r is their genetic relatedness. While many studies have provided qualitative support for Hamilton's rule, quantitative tests have not yet been possible due to the difficulty of quantifying the costs and benefits of helping acts. Here we use a simulated system of foraging robots to experimentally manipulate the costs and benefits of helping and determine the conditions under which altruism evolves. By conducting experimental evolution over hundreds of generations of selection in populations with different c/b ratios, we show that Hamilton's rule always accurately predicts the minimum relatedness necessary for altruism to evolve. This high accuracy is remarkable given the presence of pleiotropic and epistatic effects as well as mutations with strong effects on behavior and fitness (effects not directly taken into account in Hamilton's original 1964 rule). In addition to providing the first quantitative test of Hamilton's rule in a system with a complex mapping between genotype and phenotype, these experiments demonstrate the wide applicability of kin selection theory.

Using high-resolution quantitative mapping of R1 as an index of cortical myelination

A fundamental tenet of neuroscience is that cortical functional differentiation is related to the cross-areal differences in cyto-, receptor-, and myeloarchitectonics that are observed in ex-vivo preparations. An ongoing challenge is to create noninvasive magnetic resonance (MR) imaging techniques that offer sufficient resolution, tissue contrast, accuracy and precision to allow for characterization of cortical architecture over an entire living human brain. One exciting development is the advent of fast, high-resolution quantitative mapping of basic MR parameters that reflect cortical myeloarchitecture. Here, we outline some of the theoretical and technical advances underlying this technique, particularly in terms of measuring and correcting for transmit and receive radio frequency field inhomogeneities. We also discuss new directions in analytic techniques, including higher resolution reconstructions of the cortical surface. We then discuss two recent applications of this technique. The first compares individual and group myelin maps to functional retinotopic maps in the same individuals, demonstrating a close relationship between functionally and myeloarchitectonically defined areal boundaries (as well as revealing an interesting disparity in a highly studied visual area). The second combines tonotopic and myeloarchitectonic mapping to localize primary auditory areas in individual healthy adults, using a similar strategy as combined electrophysiological and post-mortem myeloarchitectonic studies in non-human primates.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Detection of live and antibiotic-killed bacteria by quantitative real-time PCR of specific fragments of rRNA.

Assessing bacterial viability by molecular markers might help accelerate the measurement of antibiotic-induced killing. This study investigated whether rRNA could be suitable for this purpose. Cultures of penicillin-susceptible and penicillin-tolerant (Tol1 mutant) Streptococcus gordonii were exposed to mechanistically different penicillin and levofloxacin. Bacterial survival was assessed by viable counts and compared to quantitative real-time PCR amplification of either the 16S rRNA genes or the 16S rRNA, following reverse transcription. Penicillin-susceptible S. gordonii lost > or =4 log(10) CFU/ml of viability over 48 h of penicillin treatment. In comparison, the Tol1 mutant lost < or =1 log(10) CFU/ml. Amplification of a 427-bp fragment of 16S rRNA genes yielded amplicons that increased proportionally to viable counts during bacterial growth but did not decrease during drug-induced killing. In contrast, the same 427-bp fragment amplified from 16S rRNA paralleled both bacterial growth and drug-induced killing. It also differentiated between penicillin-induced killing of the parent and the Tol1 mutant (> or =4 log(10) CFU/ml and < or =1 log(10) CFU/ml, respectively) and detected killing by mechanistically unrelated levofloxacin. Since large fragments of polynucleotides might be degraded faster than smaller fragments, the experiments were repeated by amplifying a 119-bp region internal to the original 427-bp fragment. The amount of 119-bp amplicons increased proportionally to viability during growth but remained stable during drug treatment. Thus, 16S rRNA was a marker of antibiotic-induced killing, but the size of the amplified fragment was critical for differentiation between live and dead bacteria.

Osteoporotic fracture risk in elderly women: estimation with quantitative heel US and clinical risk factors

PURPOSE: To derive a prediction rule by using prospectively obtained clinical and bone ultrasonographic (US) data to identify elderly women at risk for osteoporotic fractures. MATERIALS AND METHODS: The study was approved by the Swiss Ethics Committee. A prediction rule was computed by using data from a 3-year prospective multicenter study to assess the predictive value of heel-bone quantitative US in 6174 Swiss women aged 70-85 years. A quantitative US device to calculate the stiffness index at the heel was used. Baseline characteristics, known risk factors for osteoporosis and fall, and the quantitative US stiffness index were used to elaborate a predictive rule for osteoporotic fracture. Predictive values were determined by using a univariate Cox model and were adjusted with multivariate analysis. RESULTS: There were five risk factors for the incidence of osteoporotic fracture: older age (>75 years) (P < .001), low heel quantitative US stiffness index (<78%) (P < .001), history of fracture (P = .001), recent fall (P = .001), and a failed chair test (P = .029). The score points assigned to these risk factors were as follows: age, 2 (3 if age > 80 years); low quantitative US stiffness index, 5 (7.5 if stiffness index < 60%); history of fracture, 1; recent fall, 1.5; and failed chair test, 1. The cutoff value to obtain a high sensitivity (90%) was 4.5. With this cutoff, 1464 women were at lower risk (score, <4.5) and 4710 were at higher risk (score, >or=4.5) for fracture. Among the higher-risk women, 6.1% had an osteoporotic fracture, versus 1.8% of women at lower risk. Among the women who had a hip fracture, 90% were in the higher-risk group. CONCLUSION: A prediction rule obtained by using quantitative US stiffness index and four clinical risk factors helped discriminate, with high sensitivity, women at higher versus those at lower risk for osteoporotic fracture.

Bilateral Sequential NAION: A Comparison of Visual Outcomes in Fellow Eyes Using Quantitative Analysis of Goldmann Visual Fields

Purpose: Previous studies of the visual outcome in bilateral non-arteritic anterior ischemic optic neuropathy (NAION) have yielded conflicting results, specifically regarding congruity between fellow eyes. Prior studies have used measures of acuity and computerized perimetry but none has compared Goldmann visual field outcomes between fellow eyes. In order to better define the concordance of visual loss in this condition, we reviewed our cases of bilateral sequential NAION, including measures of visual acuity, pupillary function and both pattern and severity of visual field loss.Methods: We performed a retrospective chart review of 102 patients with a diagnosis of bilateral sequential NAION. Of the 102 patients, 86 were included in the study for analysis of final visual outcome between the affected eyes. Visual function was assessed using visual acuity, Goldmann visual fields, color vision and RAPD. A quantitative total visual field score and score per quadrant was analyzed for each eye using the numerical Goldmann visual field scoring method previously described by Esterman and colleagues. Based upon these scores, we calculated the total deviation and pattern deviation between fellow eyes and between eyes of different patients. Statistical significance was determined using nonparametric tests.Results: A statistically significant correlation was found between fellow eyes for multiple parameters, including logMAR visual acuity (P = 0.0101), global visual field (P = 0.0001), superior visual field (P = 0.0001), and inferior visual field (P = 0.0001). In addition, the mean deviation of both total (P = 0.0000000007) and pattern (P = 0.000000004) deviation analyses was significantly less between fellow eyes ("intra"-eyes) than between eyes of different patients ("inter"-eyes).Conclusions: Visual function between fellow eyes showed a fair to moderate correlation that was statistically significant. The pattern of vision loss was also more similar in fellow eyes than between eyes of different patients. These results may help allow better prediction of visual outcome for the second eye in patients with NAION. These findings may also be useful for evaluating efficacy of therapeutic interventions.

Viscoelasticity in Achilles Tendonopathy: Quantitative Assessment by Using Real-time Shear-Wave Elastography.

Purpose To investigate the differences in viscoelastic properties between normal and pathologic Achilles tendons ( AT Achilles tendon s) by using real-time shear-wave elastography ( SWE shear-wave elastography ). Materials and Methods The institutional review board approved this study, and written informed consent was obtained from 25 symptomatic patients and 80 volunteers. One hundred eighty ultrasonographic (US) and SWE shear-wave elastography studies of AT Achilles tendon s without tendonopathy and 30 studies of the middle portion of the AT Achilles tendon in patients with tendonopathy were assessed prospectively. Each study included data sets acquired at B-mode US (tendon morphology and cross-sectional area) and SWE shear-wave elastography (axial and sagittal mean velocity and relative anisotropic coefficient) for two passively mobilized ankle positions. The presence of AT Achilles tendon tears at B-mode US and signal-void areas at SWE shear-wave elastography were noted. Results Significantly lower mean velocity was shown in tendons with tendonopathy than in normal tendons in the relaxed position at axial SWE shear-wave elastography (P < .001) and in the stretched position at sagittal (P < .001) and axial (P = .0026) SWE shear-wave elastography . Tendon softening was a sign of tendonopathy in relaxed AT Achilles tendon s when the mean velocity was less than or equal to 4.06 m · sec(-1) at axial SWE shear-wave elastography (sensitivity, 54.2%; 95% confidence interval [ CI confidence interval ]: 32.8, 74.4; specificity, 91.5%; 95% CI confidence interval : 86.3, 95.1) and less than or equal to 5.70 m · sec(-1) at sagittal SWE shear-wave elastography (sensitivity, 41.7%; 95% CI confidence interval : 22.1, 63.3; specificity, 81.8%; 95% CI confidence interval : 75.3, 87.2) and in stretched AT Achilles tendon s, when the mean velocity was less than or equal to 4.86 m · sec(-1) at axial SWE shear-wave elastography (sensitivity, 66.7%; 95% CI confidence interval : 44.7, 84.3; specificity, 75.6%; 95% CI confidence interval : 68.5, 81.7) and less than or equal to 14.58 m · sec(-1) at sagittal SWE shear-wave elastography (sensitivity, 58.3%; 95% CI confidence interval : 36.7, 77.9; specificity, 83.5%; 95% CI confidence interval : 77.2, 88.7). Anisotropic results were not significantly different between normal and pathologic AT Achilles tendon s. Six of six (100%) partial-thickness tears appeared as signal-void areas at SWE shear-wave elastography . Conclusion Whether the AT Achilles tendon was relaxed or stretched, SWE shear-wave elastography helped to confirm and quantify pathologic tendon softening in patients with tendonopathy in the midportion of the AT Achilles tendon and did not reveal modifications of viscoelastic anisotropy in the tendon. Tendon softening assessed by using SWE shear-wave elastography appeared to be highly specific, but sensitivity was relatively low. © RSNA, 2014.

Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy.

We have developed a digital holographic microscope (DHM), in a transmission mode, especially dedicated to the quantitative visualization of phase objects such as living cells. The method is based on an original numerical algorithm presented in detail elsewhere [Cuche et al., Appl. Opt. 38, 6994 (1999)]. DHM images of living cells in culture are shown for what is to our knowledge the first time. They represent the distribution of the optical path length over the cell, which has been measured with subwavelength accuracy. These DHM images are compared with those obtained by use of the widely used phase contrast and Nomarski differential interference contrast techniques.

Qualitative and quantitative analysis of thalamo-cortical axons in the somatosensory cortex of normal and barrelless mice

Résumé Les rongeurs utilisent leurs moustaches (vibrisses) pour explorer le milieu environnant. Chaque moustache est mue par un système des muscles. Les récepteurs situés à sa base transmettent les informations au système nerveux central. La transmission vers l'écorce se fait via trois neurones de relais qui se trouvent au niveau du ganglion trigéminé, du tronc cérébral et du thalamus. La représentation corticale d'une vibrisse est une concentration des axones thalamo-corticaux (ATC) autour desquelles s'organisent leurs cibles, les cellules de la couche IV. La structure peut être identifiée histologiquement en coupes tangentielles et porte le nom de « barrel » (« tonneau »). Cette correspondance vibrisse - barrel fait de ce système un model idéal pour étudier l'influence de l'activité périphérique sur l'établissement et le maintien des cartes somatotopiques. Notre laboratoire dispose d'une souche de souris qui a subi une mutation spontanée pour le gène codant l'adenylyl cyclase I (ACI). Cette enzyme membranaire catalyse la formation de l'AMPc et joue un rôle important dans le guidage axonal, la libération des neurotransmetteurs et l'intégration des signaux postsynaptiques. Nous avons démontré dans un premier temps que cette souris adulte ne développe pas de barrels. Cela est dû à un manque d'organisation des ATC et aussi des cellules de la couche IV. De plus, les résultats électrophysiologiques montrent que les informations venant des vibrisses adjacentes ne sont pas intégrées d'une manière normale. Dans ce travail de thèse, j'ai analysé la morphologie des ATC révélés individuellement avec de la biocytine. L'analyse quantitative des ATC a mis en évidence les points suivants: 1. Les axones de la souris normale (NOR) quittent le thalamus, traversent la capsule interne et la substance blanche sous-corticale et pénètrent dans le cortex somato-sensoriel primaire. A l'intérieur de l'écorce ils traversent au maximum 3 colonnes corticales adjacentes dont une contient le barrel cible. En passant à travers les couches VI et V, ces axones arborisent et convergent progressivement vers le barrel dans lequel ils forment une riche arborisation. Un petit nombre des branches « errantes », pleines de boutons synaptiques, pénètrent dans les barrels voisins. Deux axones NOR provenant de corps cellulaires très proches dans le thalamus peuvent avoir un cheminement très divergent lors de la traversée de la capsule interne et de la substance blanche sous-corticale mais, à leur entrée dans le cortex, ils sont distants d'au maximum 2 colonnes corticales de la colonne qui contient le barrel cible et ils convergent progressivement vers ce barrel. 2. Les axones de la souris mutante (BRL) ont le même trajet sous-cortical que les axones NOR, mais leur entrée dans le cortex somato-sensoriel primaire est aléatoire. A l'interface entre la substance blanche sous-corticale et le cortex, l'axone principal se divise rapidement en troncs axonaux qui traversent les couches VI et V d'une manière divergente pour arriver dans la couche IV. Cela contraste beaucoup avec la trajectoire des NOR qui convergent graduellement vers leur barrel cible. Le nombre de branches radiales que les axones BRL utilisent pour entrer dans le cortex et dans la couche IV est double par rapport aux axones NOR. Parmi ces branches, seules quelques-unes donnent des arborisations, les autres ne sont pas développées et leur morphologie est semblable à celle des branches formées par les axones de la souris normale lors du développement. Deux axones BRL issus de corps cellulaires proches dans le thalamus peuvent avoir une trajectoire très divergente jusqu'à leur entrée dans la couche IV, mais à ce niveau ils sont réorientés pour se retrouver et faire un nombre maximal de branches et boutons synaptiques dans la même région corticale. Dans un cas extrême, un des axones observés est entré dans le cortex à la limite entre l'aire somatosensorielle primaire et secondaire et a parcouru une distance de 2 mm pour retrouver son partenaire thalamique et donner avec celui-ci un nombre maximal de branches dans la même région de la couche IV. 3. Les mesures quantitatives ont montré que les arborisations corticales des axones NOR ont une longueur moyenne de 18mm et sont formées par 200 segments qui portent 1200 boutons synaptiques. Par rapport à la souris NOR, les axones BRL ont en moyenne la même longueur, le même nombre de segments et boutons synaptiques, mais donnent deux fois plus de branches radiales. La surface tangentielle occupée par les arborisations BRL dans la couche IV est 2 fois plus grande que celle des NOR. Cela signifie que les 1000 boutons synaptiques qui caractérisent les arborisations NOR et BRL dans la couche IV sont disséminés sur une surface tangentielle double chez les derniers, et donc que la densité des boutons par unité de surface corticale est en moyenne plus faible. En effet, l'augmentation de la surface corticale tangentielle des BRL est due aux surfaces de faible et moyenne densité synaptique (0 - 8 boutons / 400pn2) qui augmentent 2 fois tandis que les surfaces de haute densité synaptiques (8 - 64 boutons / 4001.tm2) sont les mêmes. Nous émettons l'hypothèse selon laquelle, durant le développement, les ATC de la souris BRL divergent et forment un nombre exubérant de branches. Grâce à cette divergence et aux branches supranuméraires, ils trouvent l'endroit de l'écorce où se trouvent leurs voisins thalamiques et arborisent abondamment dans cette région. Cependant, le déficit en AGI ne leurs permet pas par la suite, sous influence de l'activité périphérique, de retirer les branches qui se trouvent dans les endroits inappropriés de l'écorce, avec de possibles conséquences sur la discrimination tactile.

Temperature-dependent turnovers in sex-determination mechanisms: a quantitative model.

Sex determination is often seen as a dichotomous process: individual sex is assumed to be determined either by genetic (genotypic sex determination, GSD) or by environmental factors (environmental sex determination, ESD), most often temperature (temperature sex determination, TSD). We endorse an alternative view, which sees GSD and TSD as the ends of a continuum. Both effects interact a priori, because temperature can affect gene expression at any step along the sex-determination cascade. We propose to define sex-determination systems at the population- (rather than individual) level, via the proportion of variance in phenotypic sex stemming from genetic versus environmental factors, and we formalize this concept in a quantitative-genetics framework. Sex is seen as a threshold trait underlain by a liability factor, and reaction norms allow modeling interactions between genotypic and temperature effects (seen as the necessary consequences of thermodynamic constraints on the underlying physiological processes). As this formalization shows, temperature changes (due to e.g., climatic changes or range expansions) are expected to provoke turnovers in sex-determination mechanisms, by inducing large-scale sex reversal and thereby sex-ratio selection for alternative sex-determining genes. The frequency of turnovers and prevalence of homomorphic sex chromosomes in cold-blooded vertebrates might thus directly relate to the temperature dependence in sex-determination mechanisms.

Expression quantitative trait locus mapping across water availability environments reveals contrasting associations with genomic features in Arabidopsis.

The regulation of gene expression is crucial for an organism's development and response to stress, and an understanding of the evolution of gene expression is of fundamental importance to basic and applied biology. To improve this understanding, we conducted expression quantitative trait locus (eQTL) mapping in the Tsu-1 (Tsushima, Japan) × Kas-1 (Kashmir, India) recombinant inbred line population of Arabidopsis thaliana across soil drying treatments. We then used genome resequencing data to evaluate whether genomic features (promoter polymorphism, recombination rate, gene length, and gene density) are associated with genes responding to the environment (E) or with genes with genetic variation (G) in gene expression in the form of eQTLs. We identified thousands of genes that responded to soil drying and hundreds of main-effect eQTLs. However, we identified very few statistically significant eQTLs that interacted with the soil drying treatment (GxE eQTL). Analysis of genome resequencing data revealed associations of several genomic features with G and E genes. In general, E genes had lower promoter diversity and local recombination rates. By contrast, genes with eQTLs (G) had significantly greater promoter diversity and were located in genomic regions with higher recombination. These results suggest that genomic architecture may play an important a role in the evolution of gene expression.

A benchmark test for a quantitative assessment of simple neuron models

Several methods and algorithms have recently been proposed that allow for the systematic evaluation of simple neuron models from intracellular or extracellular recordings. Models built in this way generate good quantitative predictions of the future activity of neurons under temporally structured current injection. It is, however, difficult to compare the advantages of various models and algorithms since each model is designed for a different set of data. Here, we report about one of the first attempts to establish a benchmark test that permits a systematic comparison of methods and performances in predicting the activity of rat cortical pyramidal neurons. We present early submissions to the benchmark test and discuss implications for the design of future tests and simple neurons models

Combining bone resorption markers and heel quantitative ultrasound to discriminate between fracture cases and controls.

This nested case-control analysis of a Swiss ambulatory cohort of elderly women assessed the discriminatory power of urinary markers of bone resorption and heel quantitative ultrasound for non-vertebral fractures. The tests all discriminated between cases and controls, but combining the two strategies yielded no additional relevant information. INTRODUCTION: Data are limited regarding the combination of bone resorption markers and heel quantitative bone ultrasound (QUS) in the detection of women at risk for fracture. METHODS: In a nested case-control analysis, we studied 368 women (mean age 76.2 +/- 3.2 years), 195 with low-trauma non-vertebral fractures and 173 without, matched for age, BMI, medical center, and follow-up duration, from a prospective study designed to predict fractures. Urinary total pyridinolines (PYD) and deoxypyridinolines (DPD) were measured by high performance liquid chromatography. All women underwent bone evaluations using Achilles+ and Sahara heel QUS. RESULTS: Areas under the receiver operating-characteristic curve (AUC) for discriminative models of the fracture group, with 95% confidence intervals, were 0.62 (0.56-0.68) and 0.59 (0.53-0.65) for PYD and DPD, and 0.64 (0.58-0.69) and 0.65 (0.59-0.71) for Achilles+ and Sahara QUS, respectively. The combination of resorption markers and QUS added no significant discriminatory information to either measurement alone with an AUC of 0.66 (0.60-0.71) for Achilles+ with PYD and 0.68 (0.62-0.73) for Sahara with PYD. CONCLUSIONS: Urinary bone resorption markers and QUS are equally discriminatory between non-vertebral fracture patients and controls. However, the combination of bone resorption markers and QUS is not better than either test used alone.