Color polymorphism and sex ratio distortion in a cichlid fish as an incipient stage in sympatric speciation by sexual selection

We investigated a Lake Victoria cichlid with a complex colour polymorphism that apparently represents one original species and two incipient species, all of which are sympatric. In laboratory breeding experiments we observed sex ratio distortion in certain matings between original and incipient species. Mate choice experiments show that males of the incipient species exhibit mating preferences against the original species, and males and females of the original species exhibit strong mating preferences against the incipient species. Mating preferences might evolve by sex ratio selection to avoid matings with distorted progeny sex ratios. Phenotype frequencies in nature suggest that mating preferences translate into mating frequencies, thus restricting gene flow and exerting disruptive sexual selection between the original and incipient species. The incipient species do not differ in morphology or ecology from the original species, implying that colour polymorphism, associated with sex ratio distortion, can be an incipient stage in sympatric speciation, and that disruption of gene flow can precede ecological differentiation

Trends of mean and extreme temperature indices since 1874 at low-elevation sites in the Southern Alps

We describe the recovery of three daily meteorological records for the southern Alps (Domodossola, Riva del Garda, and Rovereto), all starting in the second half of the nineteenth century. We use these new data, along with additional records, to study regional changes in the mean temperature and extreme indices of heat waves and cold spells frequency and duration over the period 1874–2015. The records are homogenized using subdaily cloud cover observations as a constraint for the statistical model, an approach that has never been applied before in the literature. A case study based on a record of parallel observations between a traditional meteorological window and a modern screen shows that the use of cloud cover can reduce the root-mean-square error of the homogenization by up to 30% in comparison to an unaided statistical correction. We find that mean temperature in the southern Alps has increased by 1.4°C per century over the analyzed period, with larger increases in daily minimum temperatures than maximum temperatures. The number of hot days in summer has more than tripled, and a similar increase is observed in duration of heat waves. Cold days in winter have dropped at a similar rate. These trends are mainly caused by climate change over the last few decades.

Trends in mean and extreme precipitation in the Mount Kenya region from observations and reanalyses

Land and water management in semi-arid regions requires detailed information on precipitation distribution, including extremes, and changes therein. Such information is often lacking. This paper describes statistics of mean and extreme precipitation in a unique data set from the Mount Kenya region, encompassing around 50 stations with at least 30 years of data. We describe the data set, including quality control procedures and statistical break detection. Trends in mean precipitation and extreme indices calculated from these data for individual rainy seasons are compared with corresponding trends in reanalysis products. From 1979 to 2011, mean precipitation decreased at 75% of the stations during the ‘long rains’ (March to May) and increased at 70% of the stations during the ‘short rains’ (October to December). Corresponding trends are found in the number of heavy precipitation days, and maximum of consecutive 5-day precipitation. Conversely, an increase in consecutive dry days within both main rainy seasons is found. However, trends are only statistically significant in very few cases. Reanalysis data sets agree with observations with respect to interannual variability, while correlations are considerably lower for monthly deviations (ratios) from the mean annual cycle. While some products well reproduce the rainfall climatology and some the spatial trend pattern, no product reproduces both.

Evaluation of a Bayesian MCMC Random Effects Inference Methodology for Capture-Mark-Recapture Data

Monte Carlo simulation was used to evaluate properties of a simple Bayesian MCMC analysis of the random effects model for single group Cormack-Jolly-Seber capture-recapture data. The MCMC method is applied to the model via a logit link, so parameters p, S are on a logit scale, where logit(S) is assumed to have, and is generated from, a normal distribution with mean μ and variance σ2 . Marginal prior distributions on logit(p) and μ were independent normal with mean zero and standard deviation 1.75 for logit(p) and 100 for μ ; hence minimally informative. Marginal prior distribution on σ2 was placed on τ2=1/σ2 as a gamma distribution with α=β=0.001 . The study design has 432 points spread over 5 factors: occasions (t) , new releases per occasion (u), p, μ , and σ . At each design point 100 independent trials were completed (hence 43,200 trials in total), each with sample size n=10,000 from the parameter posterior distribution. At 128 of these design points comparisons are made to previously reported results from a method of moments procedure. We looked at properties of point and interval inference on μ , and σ based on the posterior mean, median, and mode and equal-tailed 95% credibility interval. Bayesian inference did very well for the parameter μ , but under the conditions used here, MCMC inference performance for σ was mixed: poor for sparse data (i.e., only 7 occasions) or σ=0 , but good when there were sufficient data and not small σ .

Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

MOREMATA: Stata module (Mata) to provide various functions

This package includes various Mata functions. kern(): various kernel functions; kint(): kernel integral functions; kdel0(): canonical bandwidth of kernel; quantile(): quantile function; median(): median; iqrange(): inter-quartile range; ecdf(): cumulative distribution function; relrank(): grade transformation; ranks(): ranks/cumulative frequencies; freq(): compute frequency counts; histogram(): produce histogram data; mgof(): multinomial goodness-of-fit tests; collapse(): summary statistics by subgroups; _collapse(): summary statistics by subgroups; gini(): Gini coefficient; sample(): draw random sample; srswr(): SRS with replacement; srswor(): SRS without replacement; upswr(): UPS with replacement; upswor(): UPS without replacement; bs(): bootstrap estimation; bs2(): bootstrap estimation; bs_report(): report bootstrap results; jk(): jackknife estimation; jk_report(): report jackknife results; subset(): obtain subsets, one at a time; composition(): obtain compositions, one by one; ncompositions(): determine number of compositions; partition(): obtain partitions, one at a time; npartitionss(): determine number of partitions; rsubset(): draw random subset; rcomposition(): draw random composition; colvar(): variance, by column; meancolvar(): mean and variance, by column; variance0(): population variance; meanvariance0(): mean and population variance; mse(): mean squared error; colmse(): mean squared error, by column; sse(): sum of squared errors; colsse(): sum of squared errors, by column; benford(): Benford distribution; cauchy(): cumulative Cauchy-Lorentz dist.; cauchyden(): Cauchy-Lorentz density; cauchytail(): reverse cumulative Cauchy-Lorentz; invcauchy(): inverse cumulative Cauchy-Lorentz; rbinomial(): generate binomial random numbers; cebinomial(): cond. expect. of binomial r.v.; root(): Brent's univariate zero finder; nrroot(): Newton-Raphson zero finder; finvert(): univariate function inverter; integrate_sr(): univariate function integration (Simpson's rule); integrate_38(): univariate function integration (Simpson's 3/8 rule); ipolate(): linear interpolation; polint(): polynomial inter-/extrapolation; plot(): Draw twoway plot; _plot(): Draw twoway plot; panels(): identify nested panel structure; _panels(): identify panel sizes; npanels(): identify number of panels; nunique(): count number of distinct values; nuniqrows(): count number of unique rows; isconstant(): whether matrix is constant; nobs(): number of observations; colrunsum(): running sum of each column; linbin(): linear binning; fastlinbin(): fast linear binning; exactbin(): exact binning; makegrid(): equally spaced grid points; cut(): categorize data vector; posof(): find element in vector; which(): positions of nonzero elements; locate(): search an ordered vector; hunt(): consecutive search; cond(): matrix conditional operator; expand(): duplicate single rows/columns; _expand(): duplicate rows/columns in place; repeat(): duplicate contents as a whole; _repeat(): duplicate contents in place; unorder2(): stable version of unorder(); jumble2(): stable version of jumble(); _jumble2(): stable version of _jumble(); pieces(): break string into pieces; npieces(): count number of pieces; _npieces(): count number of pieces; invtokens(): reverse of tokens(); realofstr(): convert string into real; strexpand(): expand string argument; matlist(): display a (real) matrix; insheet(): read spreadsheet file; infile(): read free-format file; outsheet(): write spreadsheet file; callf(): pass optional args to function; callf_setup(): setup for mm_callf().

Altered cardiovascular rhythmicity in children with white coat and ambulatory hypertension

Adults with ambulatory hypertension or white coat hypertension (WCH) display abnormal cardiovascular rhythms. We studied cardiovascular rhythms by Fourier analysis of 24-h ambulatory blood pressure (BP) measurement profiles in 129 hypertensive children, 54 children with WCH, and 146 age-, height-, and gender-matched healthy subjects. The day/night mean arterial pressure ratio was lower in hypertensive and patients with WCH compared with controls (1.13 versus 1.16 versus 1.21, respectively; p < 0.0001). Eighty-five percent of controls were dippers compared with 74% of WCH (n.s.) and 64% of patients with ambulatory hypertension (p < 0.0001). The prevalence of 24-h rhythms was similar among the groups, but prevalence of 12-h BP rhythms was increased in hypertensive (67%) and WCH (72%) compared with controls (51%, p < 0.0001). The amplitudes of the 24-, 8-, and 6-h BP rhythms were reduced in hypertensive and WCH compared with controls (p < 0.05). Hypertensive and patients with WCH displayed delayed 24-, 12-, 8-, 6-h acrophases in comparison with controls (p < 0.05). In conclusion, hypertensive children exhibit abnormal cardiovascular rhythmicity compared with controls, especially a higher prevalence of nondipping compared with normotensive children. Abnormalities in patients with WCH are intermediate between healthy children and patients with ambulatory hypertension.

Quantitative T2 mapping of knee cartilage: differentiation of healthy control cartilage and cartilage repair tissue in the knee with unloading - initial results

Purpose: To prospectively determine on T2 cartilage maps the effect of unloading during a clinical magnetic resonance (MR) examination in the postoperative follow-up of patients after matrix-associated autologous chondrocyte transplantation (MACT) of the knee joint. Materials and Methods: Ethical approval for this study was provided by the local ethics commission, and written informed consent was obtained. Thirty patients (mean age, 35.4 years +/- 10.5) with a mean postoperative follow-up period of 29.1 months +/- 24.4 were enrolled. A multiecho spin-echo T2-weighted sequence was performed at the beginning (early unloading) and end (late unloading) of the MR examination, with an interval of 45 minutes. Mean and zonal region of interest T2 measurements were obtained in control cartilage and cartilage repair tissue. Statistical analysis of variance was performed. Results: The change in T2 values of control cartilage (early unloading, 50.2 msec +/- 8.4; late unloading, 51.3 msec +/- 8.5) was less pronounced than the change in T2 values of cartilage repair tissue (early unloading, 51.8 msec +/- 11.7; late unloading, 56.1 msec +/- 14.4) (P = .024). The difference between control cartilage and cartilage repair tissue was not significant for early unloading (P = .314) but was significant for late unloading (P = .036). Zonal T2 measurements revealed a higher dependency on unloading for the superficial cartilage layer. Conclusion: Our results suggest that T2 relaxation can be used to assess early and late unloading values of articular cartilage in a clinical setting and that the time point of the quantitative T2 measurement affects the differentiation between native and abnormal articular cartilage. (c) RSNA, 2010.

Association between conventional risk factors and different ultrasound-based markers of atherosclerosis at carotid and femoral levels in a middle-aged population

Ultrasound detection of sub-clinical atherosclerosis (ATS) may help identify individuals at high cardiovascular risk. Most studies evaluated intima-media thickness (IMT) at carotid level. We compared the relationships between main cardiovascular risk factors (CVRF) and five indicators of ATS (IMT, mean and maximal plaque thickness, mean and maximal plaque area) at both carotid and femoral levels. Ultrasound was performed on 496 participants aged 45-64 years randomly selected from the general population of the Republic of Seychelles. 73.4 % participants had ≥ 1 plaque (IMT thickening ≥ 1.2 mm) at carotid level and 67.5 % at femoral level. Variance (adjusted R2) contributed by age, sex and CVRF (smoking, LDL-cholesterol, HDL-cholesterol, blood pressure, diabetes) in predicting any of the ATS markers was larger at femoral than carotid level. At both carotid and femoral levels, the association between CVRF and ATS was stronger based on plaque-based markers than IMT. Our findings show that the associations between CVRF and ATS markers were stronger at femoral than carotid level, and with plaque-based markers rather than IMT. Pending comparison of these markers using harder cardiovascular endpoints, our findings suggest that markers based on plaque morphology assessed at femoral artery level might be useful cardiovascular risk predictors.

Effects of verum acupuncture compared to placebo acupuncture on quantitative EEG and heart rate variability in healthy volunteers

OBJECTIVES: The aim of this single-blind randomized crossover study was to evaluate specific effects of manual acupuncture on central and vegetative nervous system activity measured by quantitative electroencephalography (qEEG) and heart rate variability (HRV). DESIGN: Twenty (20) healthy volunteers (mean: 25.2 +/- 3.6 years) were monitored simultaneously using a qEEG system and a 12-channel electrocardiogram recorder during verum acupuncture (VA) at acupuncture point Large Intestine 4 (Hegu) (LI4) or placebo acupuncture (PA) at a sham point. RESULTS: In the EEG conduction of the occipital area, needle stimulation in VA increased alpha1-frequency significantly, and the ratio alpha1/theta was shifted to the benefit of alpha1 over all electrodes. The HRV parameters showed a significant increase of the low frequency/high frequency (HF) ratio during the first minute of stimulation in VA, indicating an initial increase of sympathetic activation. However, an increase of HF power in the minute after stimulation followed by a decrease in heart rate suggests delayed vagal activation. De qi (a sensation that is typical of acupuncture needling) occurred in 16 subjects during VA and in 9 volunteers during PA (80% versus 45%). CONCLUSIONS: Manual stimulation on LI4 seems to lead to specific changes in alpha EEG-frequency and in HRV parameters. A linear relationship between the HRV parameters and the alpha EEG band might point to a specific modulation of cerebral function by vegetative effects during acupuncture.

Lung mechanics and gas exchange in ventilated preterm infants during treatment of hyaline membrane disease with multiple doses of artificial surfactant (Exosurf)

Eight premature infants ventilated for hyaline membrane disease and enrolled in the OSIRIS surfactant trial were studied. Lung mechanics, gas exchange [PaCO2, arterial/alveolar PO2 ratio (a/A ratio)], and ventilator settings were determined 20 minutes before and 20 minutes after the end of Exosurf instillation, and subsequently at 12-24 hour intervals. Respiratory system compliance (Crs) and resistance (Rrs) were measured by means of the single breath occlusion method. After surfactant instillation there were no significant immediate changes in PaCO2 (36 vs. 37 mmHg), a/A ratio (0.23 vs. 0.20), Crs (0.32 vs. 0.31 mL/cm H2O/kg), and Rrs (0.11 vs. 0.16 cmH2O/mL/s) (pooled data of 18 measurement pairs). During the clinical course, mean a/A ratio improved significantly each time from 0.17 (time 0) to 0.29 (time 12-13 hours), to 0.39 (time 24-36 hours) and to 0.60 (time 48-61 hours), although mean airway pressure was reduced substantially. Mean Crs increased significantly from 0.28 mL/cmH2O/kg (time 0) to 0.38 (time 12-13 hours), to 0.37 (time 24-38 hours), and to 0.52 (time 48-61 hours), whereas mean Rrs increased from 0.10 cm H2O/mL/s (time 0) to 0.11 (time 12-13 hours), to 0.13 (time 24-36 hours) and to (time 48-61 hours) with no overall significance. A highly significant correlation was found between Crs and a/A ratio (r = 0.698, P less than 0.001). We conclude that Exosurf does not induce immediate changes in oxygenation as does the instillation of (modified) natural surfactant preparations. However, after 12 and 24 hours of treatment oxygenation and Crs improve significantly.(ABSTRACT TRUNCATED AT 250 WORDS)

Diffusion-weighted MR imaging of upper abdominal organs: field strength and intervendor variability of apparent diffusion coefficients

PURPOSE To determine the variability of apparent diffusion coefficient (ADC) values in various anatomic regions in the upper abdomen measured with magnetic resonance (MR) systems from different vendors and with different field strengths. MATERIALS AND METHODS Ten healthy men (mean age, 36.6 years ± 7.7 [standard deviation]) gave written informed consent to participate in this prospective ethics committee-approved study. Diffusion-weighted (DW) MR imaging was performed in each subject with 1.5- and 3.0-T MR systems from each of three vendors at two institutions. Two readers independently measured ADC values in seven upper abdominal regions (left and right liver lobe, gallbladder, pancreas, spleen, and renal cortex and medulla). ADC values were tested for interobserver differences, as well as for differences related to field strength and vendor, with repeated-measures analysis of variance; coefficients of variation (CVs) and variance components were calculated. RESULTS Interreader agreement was excellent (intraclass coefficient, 0.876). ADC values were (77.5-88.8) ×10(-5) mm(2)/sec in the spleen and (250.6-278.5) ×10(-5) mm(2)/sec in the gallbladder. There were no significant differences between ADC values measured at 1.5 T and those measured at 3.0 T in any anatomic region (P >.10 for all). In two of seven regions at 1.5 T (left and right liver lobes, P < .023) and in four of seven regions at 3.0 T (left liver lobe, pancreas, and renal cortex and medulla, P < .008), intervendor differences were significant. CVs ranged from 7.0% to 27.1% depending on the anatomic location. CONCLUSION Despite significant intervendor differences in ADC values of various anatomic regions of the upper abdomen, ADC values of the gallbladder, pancreas, spleen, and kidney may be comparable between MR systems from different vendors and between different field strengths.