Universal DNA methylation age across mammalian tissues

Aging, often considered a result of random cellular damage, can be accurately estimated using DNA methylation profiles, the foundation of pan-tissue epigenetic clocks. Here, we demonstrate the development of universal pan-mammalian clocks, using 11,754 methylation arrays from our Mammalian Methylation Consortium, which encompass 59 tissue types across 185 mammalian species. These predictive models estimate mammalian tissue age with high accuracy (r > 0.96). Age deviations correlate with human mortality risk, mouse somatotropic axis mutations and caloric restriction. We identified specific cytosines with methylation levels that change with age across numerous species. These sites, highly enriched in polycomb repressive complex 2-binding locations, are near genes implicated in mammalian development, cancer, obesity and longevity. Our findings offer new evidence suggesting that aging is evolutionarily conserved and intertwined with developmental processes across all mammals.

Between-population differences in the genetic and maternal components of body mass in roe deer

Background

Understanding the genetic and environmental mechanisms governing variation in morphology or phenology in wild populations is currently an important challenge. While there is a general consensus that selection is stronger under stressful conditions, it remains unclear whether the evolutionary potential of traits should increase or decrease with increasingly stressful conditions. Here, we investigate how contrasting environmental conditions during growth may affect the maternal and genetic components of body mass in roe deer, the most abundant and widespread wild ungulate in Western Europe. Body mass is a key life history trait that strongly influences both survival and reproductive performance in large herbivores. We used pedigrees and animal models to determine the variance components of juvenile and adult winter body mass in two populations experiencing contrasting early-life conditions.

Results

Our analyses showed that roe deer at Chizé, where habitat was poor and unpredictable, exhibited very low genetic variance in juvenile body mass. Instead, variance in mass was mainly driven by among-cohort differences in early-life conditions and maternal environment. In contrast, roe deer at Bogesund, where resource availability during the critical period of fawn rearing was higher, displayed a substantial level of genetic variance in body mass. We discuss the potential role of past demography and viability selection on fawn body mass on the erosion of genetic variance in the poor habitat.

Conclusions

Our study highlights the importance of accounting for both spatial (i.e. between-population variation) and temporal (i.e. cohort variation) heterogeneity in environmental conditions, especially in early life, to understand the potential for adaptive responses of wild populations to selection.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1154-9) contains supplementary material, which is available to authorized users.

The link between behavioural type and natal dispersal propensity reveals a dispersal syndrome in a large herbivore

When individuals disperse, they modify the physical and social composition of their reproductive environment, potentially impacting their fitness. The choice an individual makes between dispersal and philopatry is thus critical, hence a better understanding of the mechanisms involved in the decision to leave the natal area is crucial. We explored how combinations of behavioural (exploration, mobility, activity and stress response) and morphological (body mass) traits measured prior to dispersal were linked to the subsequent dispersal decision in 77 roe deer Capreolus capreolus fawns. Using an unusually detailed multi-trait approach, we identified two independent behavioural continuums related to dispersal. First, a continuum of energetic expenditure contrasted individuals of low mobility, low variability in head activity and low body temperature with those that displayed opposite traits. Second, a continuum of neophobia contrasted individuals that explored more prior to dispersal and were more tolerant of capture with those that displayed opposite traits. While accounting for possible confounding effects of condition-dependence (body mass), we showed that future dispersers were less neophobic and had higher energetic budgets than future philopatric individuals, providing strong support for a dispersal syndrome in this species.

The allometry between secondary sexual traits and body size is nonlinear among cervids

Allometric relationships between sexually selected traits and body size have been extensively studied in recent decades. While sexually selected traits generally display positive allometry, a few recent reports have suggested that allometric relationships are not always linear. In male cervids, having both long antlers and large size provides benefits in terms of increased mating success. However, such attributes are costly to grow and maintain, and these costs might constrain antler length from increasing at the same rate as body mass in larger species if the quantity of energy that males can extract from their environment is limiting. We tested for possible nonlinearity in the relationship between antler size and body mass (on a log-log scale) among 31 cervids and found clear deviation from linearity in the allometry of antler length. Antler length increased linearly until a male body mass threshold at approximately 110 kg. Beyond this threshold, antler length did not change with increasing mass. We discuss this evidence of nonlinear allometry in the light of life-history theory and stress the importance of testing for nonlinearity when studying allometric relationships.

Modeling reproductive trajectories of roe deer females: fixed or dynamic heterogeneity?

The relative role of dynamic and fixed heterogeneity in shaping the individual heterogeneity observed in most life-history traits remains difficult to quantify. In a recent work, Tuljapurkar et al. (2009) suggested modeling individual heterogeneity in lifetime reproductive success by a null model building reproductive trajectories from a first-order Markov chain. According to this model, among-individual differences in reproductive trajectories would be generated by the stochastic transitions among reproductive states (such as breeder and non-breeder) due to dynamic heterogeneity. In this work, we analyze the individual variation in three reproductive metrics (reproductive status, fecundity, and reproductive success) in two populations of roe deer intensively monitored using Tuljapurkar et al. (2009)’s dynamic model. Moreover, we challenge the Tuljapurkar model previously used as a biological null model to test whether the observed distribution of reproductive success over the lifetime was generated by a stochastic process by modifying two steps of the previous model to build a full stochastic model. We show that a distribution generated by the full dynamic model proposed by Tuljapurkar et al. (2009) can be consistently interpreted as only generated from a stochastic biological process provided that the probabilities of transition among reproductive states used are independent of the current reproductive state and that the positive covariation that usually occurs between survival and reproduction among individuals is removed. Only the reproductive status of roe deer females could be restricted to a stochastic process described by the full stochastic model, probably because most females (>90%) were breeders in a given year. The fecundity of roe deer females could not be adequately described by the full dynamic and full stochastic model, and the observed distribution of female reproductive success differed from the one generated by a full dynamic model in which each individual reproductive trajectory was independent of the individual lifespan (second step of the full dynamic model changed). While there was clear evidence that dynamic heterogeneity occurred and accounted for a large part of the observed among-individual variation in reproductive trajectories of roe deer females, a stochastic process did not provide a suitable model for all reproductive metrics. Consequently, models including additional fixed and dynamic traits need to be built in order to separate the relative role of fixed and dynamic heterogeneities in generating reproductive trajectories.

Assessing the impact of climate variation on survival in vertebrate populations

The impact of the ongoing rapid climate change on natural systems is a major issue for human societies. An important challenge for ecologists is to identify the climatic factors that drive temporal variation in demographic parameters, and, ultimately, the dynamics of natural populations. The analysis of long-term monitoring data at the individual scale is often the only available approach to estimate reliably demographic parameters of vertebrate populations. We review statistical procedures used in these analyses to study links between climatic factors and survival variation in vertebrate populations. We evaluated the efficiency of various statistical procedures from an analysis of survival in a population of white stork, Ciconia ciconia, a simulation study and a critical review of 78 papers published in the ecological literature. We identified six potential methodological problems: (i) the use of statistical models that are not well-suited to the analysis of long-term monitoring data collected at the individual scale; (ii) low ratios of number of statistical units to number of candidate climatic covariates; (iii) collinearity among candidate climatic covariates; (iv) the use of statistics, to assess statistical support for climatic covariates effects, that deal poorly with unexplained variation in survival; (v) spurious detection of effects due to the co-occurrence of trends in survival and the climatic covariate time series; and (vi) assessment of the magnitude of climatic effects on survival using measures that cannot be compared across case studies. The critical review of the ecological literature revealed that five of these six methodological problems were often poorly tackled. As a consequence we concluded that many of these studies generated hypotheses but only few provided solid evidence for impacts of climatic factors on survival or reliable measures of the magnitude of such impacts. We provide practical advice to solve efficiently most of the methodological problems identified. The only frequent issue that still lacks a straightforward solution was the low ratio of the number of statistical units to the number of candidate climatic covariates. In the perspective of increasing this ratio and therefore of producing more robust analyses of the links between climate and demography, we suggest leads to improve the procedures for designing field protocols and selecting a set of candidate climatic covariates. Finally, we present recent statistical methods with potential interest for assessing the impact of climatic factors on demographic parameters.

Lifetime reproductive success and composition of the home range in a large herbivore

The relationship between individual performance and nonrandom use of habitat is fundamental to ecology; however, empirical tests of this relationship remain limited, especially for higher orders of selection like that of the home range. We quantified the association between lifetime reproductive success (LRS) and variables describing lifetime home ranges during the period of maternal care (spring to autumn) for 77 female roe deer (Capreolus capreolus) at Trois-Fontaines, Champagne-Ardenne, France (1976-2000). We maintained population growth rate (adjusted to account for removals of non-focal animals) near rmax, which enabled us to define the fitness-habitat relationship in the absence of density effects. Using a negative binomial model, we showed that a roe deer's incorporation into its home range of habitat components important to food, cover, and edge (meadows, thickets, and increased density of road allowances) was significantly related to LRS. Further, LRS decreased with increasing age of naturally reclaimed meadows at the time of a deer's birth, which may have reflected a cohort effect related to, but not entirely explained by, a decline in quality of meadows through time. Predictive capacity of the selected model, estimated as the median correlation (rs) between predicted and observed LRS among deer of cross-validation samples, was 0.55. The strength of this relationship suggests that processes like selection of the site of a home range during dispersal may play a more important role in determining fitness of individuals than previously thought. Individual fitness of highly sedentary income breeders with high reproductive output such as roe deer should be more dependent on home range quality during the period of maternal care compared to capital breeders with low reproductive output. Identification of the most important habitat attributes to survival and reproduction at low density (low levels of intraspecific competition) may prove useful for defining habitat value ("intrinsic habitat value").

The effects of a severe drought on mouflon lamb survival

The mouflon population of Caroux-Espinouse, southern France, inhabits a highly seasonal area with dry summers. We monitored summer lamb survival during a severe drought in 2003, from early June to late August. The survival of 35 radio-tagged lambs over nine two-week periods was strongly affected by the timing of rainfall. Survival depended on the amount of rainfall recorded at a given 14 day period and in the previous 14-21 day period. Survival was not influenced by the exceptionally high mean daily temperature recorded during some periods. Male lamb survival (0.68) tended to be less than female survival (0.81), although not significantly, possibly because of a low sample size. The high lamb mortality (25.7%) recorded during a four-month period is much higher than previous estimates of first-year mortality (less than 10%). We recommend accounting for climatic variation in summer when studying the population dynamics of ungulates.

Observation of the decay B-->J/psietaK and search for X(3872)-->J/psieta

We report the observation of the B meson decay B+/- -->J/psietaK+/- and evidence for the decay B0-->J/psietaK0S, using 90 x 10(6) BB; events collected at the Upsilon(4S) resonance with the BABAR detector at the SLAC PEP-II e+e- asymmetric-energy storage ring. We obtain branching fractions of B(B+/- -->J/psietaK+/-) = [10.8 +/- 2.3(stat) +/- 2.4(syst)] x 10(-5) and B(B0-->J/psietaK0S) = [8.4 +/- 2.6(stat) +/- 2.7(syst)] x 10(-5). We search for the new narrow mass state, the X(3872), recently reported by the Belle Collaboration, in the decay B+/- -->X(3872)K+/-,X(3872)-->J/psieta and determine an upper limit of B[B +/- -->X(3872)K+/- -->J/psietaK+/-] < 7.7 x 10(-6) at 90% confidence level.

Measurement of the direct CP asymmetry in b-->s gamma Decays

We describe a measurement of the direct CP asymmetry between inclusive b-->s gamma and b-->s gamma decays. This asymmetry is expected to be less than 0.01 in the standard model, but could be enhanced up to about 0.10 by new physics contributions. We use a sample of 89 x 10(6) BB pairs recorded with the BABAR detector at SLAC PEP-II, from which we reconstruct a set of 12 exclusive b-->s gamma final states containing one charged or neutral kaon and one to three pions. We measure an asymmetry of A(CP)(b-->s gamma)=0.025+/-0.050(stat)+/-0.015(syst), corresponding to an allowed range of -0.06<A(CP)(b-->s gamma)<+0.11 at 90% confidence level.

Search for the decay B0-->gammagamma

We present a limit on the branching fraction for the decay B0-->gammagamma using data collected at the Upsilon(4S) resonance with the BABAR detector at the PEP-II asymmetric energy e+e- collider. Based on the observation of one event in the signal region, out of a sample of 21.3x10(6) e+e--->Upsilon(4S)-->BB decays, we establish an upper limit on the branching fraction of B(B0-->gammagamma)<1.7x10(-6) at the 90% confidence level. This result substantially improves upon existing limits.

Measurement of the B--> J/psiK*(892) decay amplitudes

We present a measurement of the decay amplitudes in B-->J/psiK*(892) channels using 20.7 fb(-1) of data collected at the Upsilon(4S) resonance with the BABAR detector at PEP-II. We measure a P-wave fraction R(perpendicular) = (16.0 +/- 3.2 +/- 1.4)% and a longitudinal polarization fraction (59.7 +/- 2.8 +/- 2.4)%. The measurement of a relative phase that is neither 0 nor pi, phi = 2.50 +/- 0.20 +/-0.08 radians, favors a departure from the factorization hypothesis. Although the decay B-->/psiK(pi) proceeds mainly via K*(892), there is also evidence for K2*(1430) and K(pi) S-wave contributions.

Measurements of the branching fractions of exclusive charmless B meson decays with eta(') or omega mesons

We present the results of searches for B decays to charmless two-body final states containing eta(') or omega mesons, based on 20.7 fb(-1) of data collected with the BABAR detector. We find the branching fractions Beta(B(+)-->eta(')K(+)) = (70+/-8+/-5) x 10(-6), Beta(B(0)-->eta(')K(0)) = (42(+13)(-11) +/- 4) x 10(-6), and Beta(B(+)-->omega pi(+)) = (6.6(+2.1)(-1.8) +/- 0.7) x 10(-6), where the first error quoted is statistical and the second is systematic. We give measurements of four additional modes for which the 90% confidence level upper limits are Beta(B(+)-->eta(')pi(+)) < 12 x 10(-6), Beta(B(+)-->omega K(+)) < 4 x 10(-6), Beta(B(0)-->omega K(0)) < 13 x 10(-6), and Beta(B(0)-->omega pi(0)) < 3 x 10(-6).

Measurement of the B(0) and B(+) meson lifetimes with fully reconstructed hadronic final states

The B(0) and B(+) meson lifetimes have been measured in e(+)e(-) annihilation data collected in 1999 and 2000 with the BABAR detector at center-of-mass energies near the Upsilon(4S) resonance. Events are selected in which one B meson is fully reconstructed in a hadronic final state while the second B meson is reconstructed inclusively. A combined fit to the B(0) and the B(+) decay time difference distributions yields tau(B(0)) = 1.546+/-0.032(stat)+/-0.022(syst) ps, tau(B(+)) = 1.673+/-0.032(stat)+/-0.023(syst) ps, and tau(B(+))/tau(B(0)) = 1.082+/-0.026(stat)+/-0.012(syst).

Measurement of J/psi production in continuum e(+)e(-) annihilations near square root of s = 10.6 GeV

The production of J/psi mesons in continuum e(+)e(-) annihilations has been studied with the BABAR detector at energies near the Upsilon(4S) resonance. The mesons are distinguished from J/psi production in B decays through their center-of-mass momentum and energy. We measure the cross section e(+)e(-)-->J/psi X to be 2.52+/-0.21+/-0.21 pb. We set a 90% C.L. upper limit on the branching fraction for direct Upsilon(4S)-->J/psi X decays at 4.7 x 10(-4).

Measurement of the decays B--> phiK and B--> phiK*

We have observed the decays B--> phiK and phiK(*) in a sample of over 45 million B mesons collected with the BABAR detector at the PEP-II collider. The measured branching fractions are B(B+--> phiK+) = (7.7(+1.6)(-1.4)+/-0.8)x10(-6), B(B0--> phiK0) = (8.1(+3.1)(-2.5)+/-0.8)x10(-6), B(B+--> phiK(*+)) = (9.7(+4.2)(-3.4)+/-1.7)x10(-6), and B(B0--> phiK(*0)) = (8.7(+2.5)(-2.1)+/-1.1)x10(-6). We also report the upper limit B(B+--> phipi(+))<1.4x10(-6) ( 90% C.L.).

Measurement of branching fractions and search for CP-violating charge asymmetries in charmless two-body B decays into pions and kaons

We present measurements, based on a sample of approximately 23x10(6) BB pairs, of the branching fractions and a search for CP-violating charge asymmetries in charmless hadronic decays of B mesons into two-body final states of kaons and pions. We find the branching fractions B(B0-->pi(+)pi(-)) = (4.1+/-1.0+/-0.7)x10(-6), B(B0-->K+pi(-)) = (16.7+/-1.6+/-1.3)x10(-6), B(B+-->K+pi(0)) = (10.8(+2.1)(-1.9)+/-1.0)x10(-6), B(B+-->K0pi(+)) = (18.2(+3.3)(-3.0)+/-2.0)x10(-6), B(B0-->K0pi(0)) = (8.2(+3.1)(-2.7)+/-1.2)x10(-6). We also report 90% confidence level upper limits for B meson decays to the pi(+)pi(0), K+K-, and K0K+ final states. In addition, charge asymmetries have been found to be consistent with zero, where the statistical precision is in the range of +/-0.10 to +/-0.18, depending on the decay mode.

Observation of CP violation in the B(0) meson system

We present an updated measurement of time-dependent CP-violating asymmetries in neutral B decays with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. This result uses an additional sample of Upsilon(4S) decays collected in 2001, bringing the data available to 32 x 10(6) BB macro pairs. We select events in which one neutral B meson is fully reconstructed in a final state containing charmonium and the flavor of the other neutral B meson is determined from its decay products. The amplitude of the CP-violating asymmetry, which in the standard model is proportional to sin2 beta, is derived from the decay time distributions in such events. The result sin2 beta = 0.59+/-0.14(stat)+/-0.05(syst) establishes CP violation in the B(0) meson system. We also determine absolute value of lambda = 0.93+/-0.09(stat)+/-0.03(syst), consistent with no direct CP violation.