Rhinovirus Genome Variation during Chronic Upper and Lower Respiratory Tract Infections.

Routine screening of lung transplant recipients and hospital patients for respiratory virus infections allowed to identify human rhinovirus (HRV) in the upper and lower respiratory tracts, including immunocompromised hosts chronically infected with the same strain over weeks or months. Phylogenetic analysis of 144 HRV-positive samples showed no apparent correlation between a given viral genotype or species and their ability to invade the lower respiratory tract or lead to protracted infection. By contrast, protracted infections were found almost exclusively in immunocompromised patients, thus suggesting that host factors rather than the virus genotype modulate disease outcome, in particular the immune response. Complete genome sequencing of five chronic cases to study rhinovirus genome adaptation showed that the calculated mutation frequency was in the range observed during acute human infections. Analysis of mutation hot spot regions between specimens collected at different times or in different body sites revealed that non-synonymous changes were mostly concentrated in the viral capsid genes VP1, VP2 and VP3, independent of the HRV type. In an immunosuppressed lung transplant recipient infected with the same HRV strain for more than two years, both classical and ultra-deep sequencing of samples collected at different time points in the upper and lower respiratory tracts showed that these virus populations were phylogenetically indistinguishable over the course of infection, except for the last month. Specific signatures were found in the last two lower respiratory tract populations, including changes in the 5'UTR polypyrimidine tract and the VP2 immunogenic site 2. These results highlight for the first time the ability of a given rhinovirus to evolve in the course of a natural infection in immunocompromised patients and complement data obtained from previous experimental inoculation studies in immunocompetent volunteers.

Genetic variation at CHRNA5-CHRNA3-CHRNB4 interacts with smoking status to influence body mass index.

BACKGROUND: Cigarette smoking is associated with lower body mass index (BMI), and a commonly cited reason for unwillingness to quit smoking is a concern about weight gain. Common variation in the CHRNA5-CHRNA3-CHRNB4 gene region (chromosome 15q25) is robustly associated with smoking quantity in smokers, but its association with BMI is unknown. We hypothesized that genotype would accurately reflect smoking exposure and that, if smoking were causally related to weight, it would be associated with BMI in smokers, but not in never smokers. METHODS: We stratified nine European study samples by smoking status and, in each stratum, analysed the association between genotype of the 15q25 SNP, rs1051730, and BMI. We meta-analysed the results (n = 24 198) and then tested for a genotype × smoking status interaction. RESULTS: There was no evidence of association between BMI and genotype in the never smokers {difference per T-allele: 0.05 kg/m(2) [95% confidence interval (95% CI): -0.05 to 0.18]; P = 0.25}. However, in ever smokers, each additional smoking-related T-allele was associated with a 0.23 kg/m(2) (95% CI: 0.13-0.31) lower BMI (P = 8 × 10(-6)). The effect size was larger in current [0.33 kg/m(2) lower BMI per T-allele (95% CI: 0.18-0.48); P = 6 × 10(-5)], than in former smokers [0.16 kg/m(2) (95% CI: 0.03-0.29); P = 0.01]. There was strong evidence of genotype × smoking interaction (P = 0.0001). CONCLUSIONS: Smoking status modifies the association between the 15q25 variant and BMI, which strengthens evidence that smoking exposure is causally associated with reduced BMI. Smoking cessation initiatives might be more successful if they include support to maintain a healthy BMI.

Variation de la mortalité périnatale selon l'heure de naissance en Suisse. [Variations in perinatal mortality according to the hour of birth in Switzerland].

In Switzerland like in most developed countries, the number of births is strongly related to the hour of the day and the day of the week: this pattern is very probably related to the current practice in obstetrics. Less expected is the fact that the perinatal mortality shows a striking circadian rhythm according to the hour of birth. The paper presents this pattern and comments some related issues.

Genetic variation in LIN28B is associated with the timing of puberty.

The timing of puberty is highly variable. We carried out a genome-wide association study for age at menarche in 4,714 women and report an association in LIN28B on chromosome 6 (rs314276, minor allele frequency (MAF) = 0.33, P = 1.5 x 10(-8)). In independent replication studies in 16,373 women, each major allele was associated with 0.12 years earlier menarche (95% CI = 0.08-0.16; P = 2.8 x 10(-10); combined P = 3.6 x 10(-16)). This allele was also associated with earlier breast development in girls (P = 0.001; N = 4,271); earlier voice breaking (P = 0.006, N = 1,026) and more advanced pubic hair development in boys (P = 0.01; N = 4,588); a faster tempo of height growth in girls (P = 0.00008; N = 4,271) and boys (P = 0.03; N = 4,588); and shorter adult height in women (P = 3.6 x 10(-7); N = 17,274) and men (P = 0.006; N = 9,840) in keeping with earlier growth cessation. These studies identify variation in LIN28B, a potent and specific regulator of microRNA processing, as the first genetic determinant regulating the timing of human pubertal growth and development.

The political economy of childcare in OECD countries : Explaining cross-national variation in spending and coverage rates

If childcare policy has become topical in most OECD countries over the last ten years or so, actual developments display huge cross-national variations. Countries like Sweden and Denmark spend around 2 per cent of GDP on this service, and provide affordable childcare places to most children below school age. At the other extreme, in Southern Europe, only around 10 per cent of this age group has access to formal daycare. Against this background, this article aims to account for cross-national variations in childcare services. It distinguishes two dependent variables: the coverage rate and the proportion of GDP spent subsidising childcare services. Using a mix of cross-sectional and pooled times-series methods, it tests a series of hypotheses concerning the determinants of the development of this policy. Its main conclusion for the coverage rate is that key factors are public spending and wage dispersion (both positive). For spending, key factors are the proportion of women in parliaments (positive) and spending on age-related policies (negative).

Extent of resection in glioblastoma - variation between molecular risk groups in rtog-0525

Whether maximal surgical resection of glioblastoma improves patient survival has been controversial, as it is difficult to perform an unbiased assessment of extent of resection (EOR) independent of other patient-specific prognostic factors. Recently, glioblastoma has been sub-classified into 4 distinct molecular risk groups (RGs), which have been validated as prognostic biomarkers in the randomized clinical trial of temozolomide dosing in glioblastoma: the Radiation Therapy Oncology Group 0525 (RTOG-0525) trial. We sought to perform exploratory analyses examining gross total resection (GTR) versus sub-total resection (STR) within these RGs in RTOG-0525 patients. Across all randomized patients, n ¼ 354 had STR and n ¼ 450 had GTR as determined by neurosurgeon operative report. GTR was not significantly associated with survival across the overall study group. A total of 725 patients had sufficient tissue for determination of molecular RG. There were no significant differences in percentage of GTR between each of the 4 RGs (P ¼ 0.64). In exploratory subgroup analyses, GTR was associated with improved survival only for patients with tumors from RG4. Hazard ratios (95% confidence intervals) were 0.52 (0.08-2.07) for RG1 (n ¼ 28, 68% GTR), 1.74 (0.75-4.05) for RG2 (n ¼ 39, 56% GTR), 1.09 (0.84-1.42) for RG3 (n ¼ 284, 56% GTR), and 1.26 (1.01-1.56) for RG4 (n ¼ 374, 55% GTR). In univariate analysis within RG4, GTR was associated with a median survival of 14.6 months vs 12.7 months for STR (P ¼ 0.0352. In a Cox model adjusting for age, KPS, and neurologic function (NF), surgery remained an independent factor within RG4: GTR (P ¼ 0.0331), age (P ¼ 0.0014), KPS (P ¼ .3289), and NF (P ¼ 0.3804). There are important cautions in the interpretation of these data, including lack of MRI confirmation of EOR, and inclusion of a range of STR (from biopsy to near-total resection). However, these exploratory results raise the possibility that upfront characterization of tumor molecular profile may allow for personalized therapeutic strategies to improve outcomes for patients with glioblastoma.

International variation in management of screen-detected ductal carcinoma in situ of the breast.

BACKGROUND: Ductal carcinoma in situ (DCIS) incidence has grown with the implementation of screening and its detection varies across International Cancer Screening Network (ICSN) countries. The aim of this survey is to describe the management of screen-detected DCIS in ICSN countries and to evaluate the potential for treatment related morbidity. METHODS: We sought screen-detected DCIS data from the ICSN countries identified during 2004-2008. We adopted standardised data collection forms and analysis and explored DCIS diagnosis and treatment processes ranging from pre-operative diagnosis to type of surgery and radiotherapy. RESULTS: Twelve countries contributed data from a total of 15 screening programmes, all from Europe except the United States of America and Japan. Among women aged 50-69years, 7,176,050 screening tests and 5324 screen-detected DCIS were reported. From 21% to 93% of DCIS had a pre-operative diagnosis (PO); 67-90% of DCIS received breast conservation surgery (BCS), and in 41-100% of the cases this was followed by radiotherapy; 6.4-59% received sentinel lymph node biopsy (SLNB) only and 0.8-49% axillary dissection (ALND) with 0.6% (range by programmes 0-8.1%) being node positive. Among BCS patients 35% received SLNB only and 4.8% received ALND. Starting in 2006, PO and SLNB use increased while ALND remained stable. SLNB and ALND were associated with larger size and higher grade DCIS lesions. CONCLUSIONS: Variation in DCIS management among screened women is wide and includes lymph node surgery beyond what is currently recommended. This indicates the presence of varying levels of overtreatment and the potential for its reduction.

Geographic variation in the frequency of isolation and fluconazole and voriconazole susceptibilities of Candida glabrata: an assessment from the ARTEMIS DISK Global Antifungal Surveillance Program.

Geographic differences in frequency and azole resistance among Candida glabrata may impact empiric antifungal therapy choice. We examined geographic variation in isolation and azole susceptibility of C. glabrata. We examined 23 305 clinical isolates of C. glabrata during ARTEMIS DISK global surveillance. Susceptibility testing to fluconazole and voriconazole was assessed by disk diffusion, and the results were grouped by geographic location: North America (NA) (2470 isolates), Latin America (LA) (2039), Europe (EU) (12 439), Africa and the Middle East (AME) (728), and Asia-Pacific (AP) (5629). Overall, C. glabrata accounted for 11.6% of 201 653 isolates of Candida and varied as a proportion of all Candida isolated from 7.4% in LA to 21.1% in NA. Decreased susceptibility (S) to fluconazole was observed in all geographic regions and ranged from 62.8% in AME to 76.7% in LA. Variation in fluconazole susceptibility was observed within each region: AP (range, 50-100% S), AME (48-86.9%), EU (44.8-88%), LA (43-92%), and NA (74.5-91.6%). Voriconazole was more active than fluconazole (range, 82.3-84.2% S) with similar regional variation. Among 22 sentinel sites participating in ARTEMIS from 2001 through 2007 (84 140 total isolates, 8163 C. glabrata), the frequency of C. glabrata isolation increased in 14 sites and the frequency of fluconazole resistance (R) increased in 11 sites over the 7-year period of study. The sites with the highest cumulative rates of fluconazole R were in Poland (22% R), the Czech Republic (27% R), Venezuela (27% R), and Greece (33% R). C. glabrata was most often isolated from blood, normally sterile body fluids and urine. There is substantial geographic and institutional variation in both frequency of isolation and azole resistance among C. glabrata. Prompt species identification and fluconazole susceptibility testing are necessary to optimize therapy for invasive candidiasis.

Interhuman transmission as a potential key parameter for geographical variation in the prevalence of Pneumocystis jirovecii dihydropteroate synthase mutations.

BACKGROUND: Pneumocystis jirovecii dihydropteroate synthase (DHPS) mutations are associated with failure of prophylaxis with sulfa drugs. This retrospective study sought to better understand the geographical variation in the prevalence of these mutations. METHODS: DHPS polymorphisms in 394 clinical specimens from immunosuppressed patients who received a diagnosis of P. jirovecii pneumonia and who were hospitalized in 3 European cities were examined using polymerase chain reaction (PCR) single-strand conformation polymorphism. Demographic and clinical characteristics were obtained from patients' medical charts. RESULTS: Of the 394 patients, 79 (20%) were infected with a P. jirovecii strain harboring one or both of the previously reported DHPS mutations. The prevalence of DHPS mutations was significantly higher in Lyon than in Switzerland (33.0% vs 7.5%; P < .001). The proportion of patients with no evidence of sulfa exposure who harbored a mutant P. jirovecii DHPS genotype was significantly higher in Lyon than in Switzerland (29.7% vs 3.0%; P < .001). During the study period in Lyon, in contrast to the Swiss hospitals, measures to prevent dissemination of P. jirovecii from patients with P. jirovecii pneumonia were generally not implemented, and most patients received suboptimal prophylaxis, the failure of which was strictly associated with mutated P. jirovecii. Thus, nosocomial interhuman transmission of mutated strains directly or indirectly from other individuals in whom selection of mutants occurred may explain the high proportion of mutations without sulfa exposure in Lyon. CONCLUSIONS: Interhuman transmission of P. jirovecii, rather than selection pressure by sulfa prophylaxis, may play a predominant role in the geographical variation in the prevalence in the P. jirovecii DHPS mutations.

Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Age-specific fitness components and their temporal variation in the barn owl.

Theory predicts that temporal variability plays an important role in the evolution of life histories, but empirical studies evaluating this prediction are rare. In constant environments, fitness can be measured by the population growth rate lambda, and the sensitivity of lambda to changes in fitness components estimates selection on these traits. In variable environments, fitness is measured by the stochastic growth rate lambda(S), and stochastic sensitivities estimate selection pressure. Here we examine age-specific schedules for reproduction and survival in a barn owl population (Tyto alba). We estimated how temporal variability affected fitness and selection, accounting for sampling variance. Despite large sample sizes of old individuals, we found no strong evidence for senescence. The most variable fitness components were associated with reproduction. Survival was less variable. Stochastic simulations showed that the observed variation decreased fitness by about 30%, but the sensitivities of lambda and lambda(S) to changes in all fitness components were almost equal, suggesting that temporal variation had negligible effects on selection. We obtained these results despite high observed variability in the fitness components and relatively short generation time of the study organism, a situation in which temporal variability should be particularly important for natural selection and early senescence is expected.

Seasonal variation of overall and cardiovascular mortality: a study in 19 countries from different geographic locations.

BACKGROUND: Cardiovascular diseases (CVD) mortality has been shown to follow a seasonal pattern. Several studies suggested several possible determinants of this pattern, including misclassification of causes of deaths. We aimed at assessing seasonality in overall, CVD, cancer and non-CVD/non-cancer mortality using data from 19 countries from different latitudes. METHODS AND FINDINGS: Monthly mortality data were compiled from 19 countries, amounting to over 54 million deaths. We calculated ratios of the observed to the expected numbers of deaths in the absence of a seasonal pattern. Seasonal variation (peak to nadir difference) for overall and cause-specific (CVD, cancer or non-CVD/non-cancer) mortality was analyzed using the cosinor function model. Mortality from overall, CVD and non-CVD/non-cancer showed a consistent seasonal pattern. In both hemispheres, the number of deaths was higher than expected in winter. In countries close to the Equator the seasonal pattern was considerably lower for mortality from any cause. For CVD mortality, the peak to nadir differences ranged from 0.185 to 0.466 in the Northern Hemisphere, from 0.087 to 0.108 near the Equator, and from 0.219 to 0.409 in the Southern Hemisphere. For cancer mortality, the seasonal variation was nonexistent in most countries. CONCLUSIONS: In countries with seasonal variation, mortality from overall, CVD and non-CVD/non-cancer show a seasonal pattern with mortality being higher in winter than in summer. Conversely, cancer mortality shows no substantial seasonality.

Geographic variation in sex-chromosome differentiation in the common frog (Rana temporaria).

In sharp contrast with birds and mammals, sex-determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex-determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500-km transect across Sweden. Genetic differentiation between sexes is strongest (FST  = 0.152) in a northern-boreal population, where male-specific alleles and heterozygote excesses (FIS  = -0.418 in males, +0.025 in females) testify to a male-heterogametic system and lack of X-Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (FST  = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male-specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X-Y recombination, co-option of an alternative sex-chromosome pair, or a mixed sex-determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the 'sexual races' described in common frogs in the 1930s.

Segmental copy number variation shapes tissue transcriptomes.

Copy number variation (CNV) is a key source of genetic diversity, but a comprehensive understanding of its phenotypic effect is only beginning to emerge. We have generated a CNV map in wild mice and classical inbred strains. Genome-wide expression data from six major organs show not only that expression of genes within CNVs tend to correlate with copy number changes, but also that CNVs influence the expression of genes in their vicinity, an effect that extends up to half a megabase. Genes within CNVs show lower expression and more specific spatial expression patterns than genes mapping elsewhere. Our analyses reveal differential constraint on copy number changes of genes expressed in different tissues. Dosage alterations of brain-expressed genes are less frequent than those of other genes and are buffered by tighter transcriptional regulation. Our study provides initial evidence that CNVs shape tissue transcriptomes on a global scale.