Integrative analyses of speciation and divergence in Psammodromus hispanicus (Squamata: Lacertidae).

BackgroundGenetic, phenotypic and ecological divergence within a lineage is the result of past and ongoing evolutionary processes, which lead ultimately to diversification and speciation. Integrative analyses allow linking diversification to geological, climatic, and ecological events, and thus disentangling the relative importance of different evolutionary drivers in generating and maintaining current species richness.ResultsHere, we use phylogenetic, phenotypic, geographic, and environmental data to investigate diversification in the Spanish sand racer (Psammodromus hispanicus). Phylogenetic, molecular clock dating, and phenotypic analyses show that P. hispanicus consists of three lineages. One lineage from Western Spain diverged 8.3 (2.9-14.7) Mya from the ancestor of Psammodromus hispanicus edwardsianus and P. hispanicus hispanicus Central lineage. The latter diverged 4.8 (1.5-8.7) Mya. Molecular clock dating, together with population genetic analyses, indicate that the three lineages experienced northward range expansions from southern Iberian refugia during Pleistocene glacial periods. Ecological niche modelling shows that suitable habitat of the Western lineage and P. h. edwardsianus overlap over vast areas, but that a barrier may hinder dispersal and genetic mixing of populations of both lineages. P. h. hispanicus Central lineage inhabits an ecological niche that overlaps marginally with the other two lineages.ConclusionsOur results provide evidence for divergence in allopatry and niche conservatism between the Western lineage and the ancestor of P. h. edwardsianus and P. h. hispanicus Central lineage, whereas they suggest that niche divergence is involved in the origin of the latter two lineages. Both processes were temporally separated and may be responsible for the here documented genetic and phenotypic diversity of P. hispanicus. The temporal pattern is in line with those proposed for other animal lineages. It suggests that geographic isolation and vicariance played an important role in the early diversification of the group, and that lineage diversification was further amplified through ecological divergence.

International perspectives on general internal medicine and the case for "globalization" of a discipline.

General internal medicine (GIM) has flourished in the United States (U.S.). Unlike other subspecialties of internal medicine, however, GIM's evolution has not been global in scope, but rather appears to have occurred in isolation within countries. Here, we describe international models of GIM from Canada, Switzerland, Australia/New Zealand, Argentina, and Japan, and compare these with the U.S. model. There are notable differences in the typical clinical roles assumed by General Internists across these 7 countries, but also important overlap in clinical and academic domains. Despite this overlap, there has been a relative lack of contact among General Internists from these and other countries at a truly international GIM meeting; the time is now for increased international exchange and the "globalization" of GIM.

On stand by: host genetics of HIV control.

The impact of host genetic variation on determining the differential outcomes after HIV infection has been studied by two approaches: targeting of candidate genes and genome-wide association studies (GWASs). The overlap in genetic variants that has been identified by these two means has essentially been restricted to variants near to the human leukocyte antigen (HLA) class I genes, although variation in the CCR5 locus, which was first shown to have an effect on HIV outcomes using the candidate gene approach, does reach significance genome-wide when very large samples sizes (i.e. thousands) are used in GWAS. Overall, many of the variants identified by the candidate gene approach are likely to be spurious, as no additional variants apart from a novel variant near the HLA-C gene have been consistently identified by GWAS. Variants with low frequency and/or low impact on HIV outcomes are likely to exist in the genome and there could be many of them, but these are not identifiable, given current GWAS sample sizes. Several loci centrally involved in the immune response, including the immunoglobulin genes, T-cell receptor loci, or leukocyte receptor complex, are either poorly covered on the GWAS chips or difficult to interpret due to their repetitive nature and/or the presence of insertion/deletion polymorphisms in the region. These loci warrant further interrogation, but genetic characterization of these regions across a range of individuals will first be required. Finally, synergistic interactions between loci may affect outcome after infection, as suggested by associations of specific, functionally relevant HLA and killer cell immunoglobulin-like receptor variants with HIV disease outcomes, and these require further consideration as well.

Hippocampal atrophy predicts conversion to dementia after STN-DBS in Parkinson's disease.

BACKGROUND: Hippocampal atrophy (HA) is a known predictor of dementia in Alzheimer's disease. HA has been found in advanced Parkinson's disease (PD), but no predicting value has been demonstrated yet. The identification of such a predictor in candidates for subthalamic deep brain stimulation (STN-DBS) would be of value. Our objective was to compare preoperative hippocampal volumes (HV) between PD patients who subsequently converted to dementia (PDD) after STN-DBS and those who did not (PDnD). METHODS: From a cohort of 70 consecutive STN-DBS treated PD patients, 14 converted to dementia over 25.6+/-20.2 months (PDD). They were compared to 14 matched controls (PDnD) who did not convert to dementia after 43.9+/-11.7 months. On the preoperative 3D MPRAGE MRI images, HV and total brain volumes (TBV) were measured by a blinded investigator using manual and automatic segmentation respectively. RESULTS: PDD had smaller preoperative HV than PDnD (1.95+/-0.29 ml; 2.28+/-0.33 ml; p<0.01). This difference reinforced after normalization for TBV (3.28+/-0.48, 3.93+/-0.60; p<0.01). Every 0.1 ml decrease of HV increased the likelihood to develop dementia by 24.6%. A large overlap was found between PD and PDnD HVs, precluding the identification of a cut-off score. CONCLUSIONS: As in Alzheimer's disease, HA may be a predictor of the conversion to dementia in PD. This preoperative predictor suggests that the development of dementia after STN-DBS is related to the disease progression, rather then the procedure. Further studies are needed to define a cut-off score for HA, in order to affine its predictive value for an individual patient.

Taxonomic status of Hylomys parvus and Hylomys suillus (Insectivora: Erinaceidae): biochemical and morphological analyses

The genus Hylomys was thought to be represented by a single widespread species. Biochemical and morphometric analyses of several Southeast Asian populations reveal that Sumatra is inhabited by two distinct species, the dwarf gymnure (H. parvus) and the lesser gymnure (H. suillus). The absence of interbreeding between these two groups along with their relatively ancient common origins are documented by several diagnostic loci and a large Nei's genetic distance (D = 0.353 +/- 0.035). The dwarf gymnure has been reported only from the slopes of the Mt. Kerinci volcano in Sumatra, where the species lives at higher elevations than its potential competitor, the lesser gymnure. Other populations of Hylomys from Java, Borneo, and Malaysia are more closely related to the Sumatran sample of H. suillus, but they exhibit strong interpopulational genetic differentiation (D = 0.165 +/- 0.040) that may be accounted for by their isolated montane habitat. In addition, a principal-components analysis based on 16 measurements of the skull clearly separates adult specimens of both species. There is little overlap in the measurements between H. suillus (which is larger) and H. parvus. On Sumatra where both species may be sympatric, the notched space between premaxillary tips, soft texture of the fur, and more delicate skull and dentition are diagnostic of H. parvus.

Genome-wide profiling of blood pressure in adults and children.

Hypertension is an important determinant of cardiovascular morbidity and mortality and has a substantial heritability, which is likely of polygenic origin. The aim of this study was to assess to what extent multiple common genetic variants contribute to blood pressure regulation in both adults and children and to assess overlap in variants between different age groups, using genome-wide profiling. Single nucleotide polymorphism sets were defined based on a meta-analysis of genome-wide association studies on systolic blood pressure and diastolic blood pressure performed by the Cohort for Heart and Aging Research in Genome Epidemiology (n=29 136), using different P value thresholds for selecting single nucleotide polymorphisms. Subsequently, genetic risk scores for systolic blood pressure and diastolic blood pressure were calculated in an independent adult population (n=2072) and a child population (n=1034). The explained variance of the genetic risk scores was evaluated using linear regression models, including sex, age, and body mass index. Genetic risk scores, including also many nongenome-wide significant single nucleotide polymorphisms, explained more of the variance than scores based only on very significant single nucleotide polymorphisms in adults and children. Genetic risk scores significantly explained ≤1.2% (P=9.6*10(-8)) of the variance in adult systolic blood pressure and 0.8% (P=0.004) in children. For diastolic blood pressure, the variance explained was similar in adults and children (1.7% [P=8.9*10(-10)] and 1.4% [P=3.3*10(-5)], respectively). These findings suggest the presence of many genetic loci with small effects on blood pressure regulation both in adults and children, indicating also a (partly) common polygenic regulation of blood pressure throughout different periods of life.

Habitat selection in zones of parapatric contact between the common shrew Sorex araneaus and Millet's shrew S. coronatus

(1) The common shrew Sorex araneus and Millet's shrew S. coronatusare sibling species.They are morphologically and genetically very similar but do not hybridize.Their parapatric distribution throughout south-western Europe, with a few narrow zones of distributional overlap, suggests that they are in competitive parapatry. (2) Two of these contact zones were studied; there was evidence of coexistence over periods of 2 years as well as habitat segregation. In both zones, the species segregated on litter thickness and humidity variables. (3) A simple analysis of spatial distribution showed that habitats visible in the field corresponded to the habitats selected by the species. Habitat selection was found throughout the annual life-cycle of the shrews. (4) In one contact zone, a removal experiment was performed to test whether habitat segregation is induced by interspecific interactions. The experiment showed that the species select habitats differentially when both are present and abandon habitat selection when their competitor removed. (5)These results confirm the role of resource partitioning in promoting narrow ranges of distributional overlap between such parapatric species and qualitatively support the prediction of habitat selection theory that, in a two-species system, coexistence may be achieved by differential habitat selection to avoid competition. The results also support the view that the common shrew and Millet's shrew are in competitive parapatry.

An operator formulation of the multiscale finite-volume method with correction function

The multiscale finite-volume (MSFV) method has been derived to efficiently solve large problems with spatially varying coefficients. The fine-scale problem is subdivided into local problems that can be solved separately and are coupled by a global problem. This algorithm, in consequence, shares some characteristics with two-level domain decomposition (DD) methods. However, the MSFV algorithm is different in that it incorporates a flux reconstruction step, which delivers a fine-scale mass conservative flux field without the need for iterating. This is achieved by the use of two overlapping coarse grids. The recently introduced correction function allows for a consistent handling of source terms, which makes the MSFV method a flexible algorithm that is applicable to a wide spectrum of problems. It is demonstrated that the MSFV operator, used to compute an approximate pressure solution, can be equivalently constructed by writing the Schur complement with a tangential approximation of a single-cell overlapping grid and incorporation of appropriate coarse-scale mass-balance equations.

Complete nucleotide sequence of Birmingham IncP alpha plasmids. Compilation and comparative analysis.

The IncP alpha promiscuous plasmid (R18, R68, RK2, RP1 and RP4) comprises 60,099 bp of nucleotide sequence, encoding at least 74 genes. About 40 kb of the genome, designated the IncP core and including all essential replication and transfer functions, can be aligned with equivalent sequences in the IncP beta plasmid R751. The compiled IncP alpha sequence revealed several previously unidentified reading frames that are potential genes. IncP alpha plasmids carry genetic information very efficiently: the coding sequences of the genes are closely packed but rarely overlap, and occupy almost 86% of the genome's nucleotide sequence. All of the 74 genes should be expressed, although there is as yet experimental evidence for expression of only 60 of them. Six examples of tandem-in-frame initiation sites specifying two gene products each are known. Two overlapping gene arrangements occupy different reading frames of the same region. Intergenic regions include most of the 25 promoters; transcripts are usually polycistronic. Translation of most of the open reading frames seems to be initiated independently, each from its own ribosomal binding and initiation site, although, a few cases of coupled translation have been reported. The most frequently used initiation codon is AUG but translation for a few open reading frames begins at GUG or UUG. The most common stop-codon is UGA followed by UAA and then UAG. Regulatory circuits are complex and largely dependent on two components of the central control operon. KorA and KorB are transcriptional repressors controlling at least seven operons. KorA and KorB act synergistically in several cases by recognizing and binding to conserved nucleotide sequences. Twelve KorB binding sites were found around the IncP alpha sequence and these are conserved in R751 (IncP beta) with respect to both sequence and location. Replication of IncP alpha plasmids requires oriV and the plasmid-encoded initiator protein TrfA in combination with the host-encoded replication machinery. Conjugative plasmid transfer depends on two separate regions occupying about half of the genome. The primary segregational stability system designated Par/Mrs consists of a putative site-specific recombinase, a possible partitioning apparatus and a post-segregational lethality mechanism, all encoded in two divergent operons. Proteins related to the products of F sop and P1 par partitioning genes are separately encoded in the central control operon.

Residence time, expansion toward the equator in the invaded range and native range size matter to climatic niche shifts in non-native species

Aim Identifying climatic niche shifts and their drivers is important to accurately predict the risk of biological invasions. The niches of non-native plants and birds have recently been assessed in large-scale multi-species studies, but such large-scale tests are lacking for non-native reptiles and amphibians (herpetofauna). Furthermore, little is known about the factors contributing to niche shifts when they occur. Based on the occurrence of 71 reptile and amphibian species, we compared native and non-native realized niches in 101 invaded ranges at a worldwide scale and identified the factors that affect niche shifts. Location The world except the Antarctic. Methods We assessed climatic niche dynamics in a gridded environmental space allowing the quantification of niche overlap and expansion into climatic conditions not colonized by the species in their native range. We analyzed the factors affecting niche shifts using a model averaging approach based on generalized linear mixed-effects models. Results Approximately 57% of the invaded ranges (51% for amphibians and 61% for reptiles) showed niche shifts (≥10% expansion in the realized climatic niche). Island endemics, species introduced to Oceania and invaded ranges outside the native biogeographic realm showed a higher proportion of niche shifts. Niche shifts were more likely for species that had smaller native range sizes, were introduced earlier into a new range or invaded areas located at lower latitudes than the native range. Main conclusions The proportion of niche shifts for non-native herpetofauna was higher than those for Holarctic non-native plants and European non-native birds. The 'climate matching hypothesis' should be used with caution for species shifting their niche because it could underestimate the risk of their establishment.

Testing the competitive exclusion principle using various niche parameters in a native (Natrix maura) and an introduced (N. tessellata) colubrid

Despite the increase of animal and plant introductions worldwide and the strong augmentation of the reptile trade, few invasive snake populations have been studied. Dice snakes (Natrix tessellata) were introduced to the shores of Lake Geneva (Switzerland) in the early 1920s, and are now well established. This region of introduction was previously inhabited by Viperine snakes (N. maura). Ever since these two species have been under monitoring (which began in 1996) the Viperine snake population has shown drastic decline. We examine here the possibility of trophic competition by analysing diet composition, prey size and trophic niche overlap. Spatial distribution is also assessed in order to address the question of spatial competitive exclusion. We found very similar diets, and thus a high trophic niche overlap, indicating no partitioning of the trophic resource. No arguments in favour of spatial competitive exclusion were found. Our study suggests that trophic competition may occur between the two natricines and that it may give an explanation for the drastic decline of the Viperine snake in this area. Other pathways potentially playing a role in the exclusion of the Viperine snake are discussed.

Magnetoencephalography demonstrates multiple asynchronous generators during human sleep spindles.

Sleep spindles are approximately 1 s bursts of 10-16 Hz activity that occur during stage 2 sleep. Spindles are highly synchronous across the cortex and thalamus in animals, and across the scalp in humans, implying correspondingly widespread and synchronized cortical generators. However, prior studies have noted occasional dissociations of the magnetoencephalogram (MEG) from the EEG during spindles, although detailed studies of this phenomenon have been lacking. We systematically compared high-density MEG and EEG recordings during naturally occurring spindles in healthy humans. As expected, EEG was highly coherent across the scalp, with consistent topography across spindles. In contrast, the simultaneously recorded MEG was not synchronous, but varied strongly in amplitude and phase across locations and spindles. Overall, average coherence between pairs of EEG sensors was approximately 0.7, whereas MEG coherence was approximately 0.3 during spindles. Whereas 2 principle components explained approximately 50% of EEG spindle variance, >15 were required for MEG. Each PCA component for MEG typically involved several widely distributed locations, which were relatively coherent with each other. These results show that, in contrast to current models based on animal experiments, multiple asynchronous neural generators are active during normal human sleep spindles and are visible to MEG. It is possible that these multiple sources may overlap sufficiently in different EEG sensors to appear synchronous. Alternatively, EEG recordings may reflect diffusely distributed synchronous generators that are less visible to MEG. An intriguing possibility is that MEG preferentially records from the focal core thalamocortical system during spindles, and EEG from the distributed matrix system.

Phenotip - a web-based instrument to help diagnosing fetal syndromes antenatally.

Prenatal ultrasound can often reliably distinguish fetal anatomic anomalies, particularly in the hands of an experienced ultrasonographer. Given the large number of existing syndromes and the significant overlap in prenatal findings, antenatal differentiation for syndrome diagnosis is difficult. We constructed a hierarchic tree of 1140 sonographic markers and submarkers, organized per organ system. Subsequently, a database of prenatally diagnosable syndromes was built. An internet-based search engine was then designed to search the syndrome database based on a single or multiple sonographic markers. Future developments will include a database with magnetic resonance imaging findings as well as further refinements in the search engine to allow prioritization based on incidence of syndromes and markers.

On the tractability of the piecewise-linear approximation for general discrete-choice network revenue management

The choice network revenue management (RM) model incorporates customer purchase behavioras customers purchasing products with certain probabilities that are a function of the offeredassortment of products, and is the appropriate model for airline and hotel network revenuemanagement, dynamic sales of bundles, and dynamic assortment optimization. The underlyingstochastic dynamic program is intractable and even its certainty-equivalence approximation, inthe form of a linear program called Choice Deterministic Linear Program (CDLP) is difficultto solve in most cases. The separation problem for CDLP is NP-complete for MNL with justtwo segments when their consideration sets overlap; the affine approximation of the dynamicprogram is NP-complete for even a single-segment MNL. This is in contrast to the independentclass(perfect-segmentation) case where even the piecewise-linear approximation has been shownto be tractable. In this paper we investigate the piecewise-linear approximation for network RMunder a general discrete-choice model of demand. We show that the gap between the CDLP andthe piecewise-linear bounds is within a factor of at most 2. We then show that the piecewiselinearapproximation is polynomially-time solvable for a fixed consideration set size, bringing itinto the realm of tractability for small consideration sets; small consideration sets are a reasonablemodeling tradeoff in many practical applications. Our solution relies on showing that forany discrete-choice model the separation problem for the linear program of the piecewise-linearapproximation can be solved exactly by a Lagrangian relaxation. We give modeling extensionsand show by numerical experiments the improvements from using piecewise-linear approximationfunctions.

Selective resonance suppression 1H-[13C] NMR spectroscopy with asymmetric adiabatic RF pulses.

Despite obvious improvements in spectral resolution at high magnetic field, the detection of 13C labeling by 1H-[13C] NMR spectroscopy remains hampered by spectral overlap, such as in the spectral region of 1H resonances bound to C3 of glutamate (Glu) and glutamine (Gln), and C6 of N-acetylaspartate (NAA). The aim of this study was to develop, implement, and apply a novel 1H-[13C] NMR spectroscopic editing scheme, dubbed "selective Resonance suppression by Adiabatic Carbon Editing and Decoupling single-voxel STimulated Echo Acquisition Mode" (RACED-STEAM). The sequence is based on the application of two asymmetric narrow-transition-band adiabatic RF inversion pulses at the resonance frequency of the 13C coupled to the protons that need to be suppressed during the mixing time (TM) period, alternating the inversion band downfield and upfield from the 13C resonance on odd and even scans, respectively, thus suppressing the detection of 1H resonances bound to 13C within the transition band of the inversion pulse. The results demonstrate the efficient suppression of 1H resonances bound to C3 of Glu and Gln, and C4 of Glu, which allows the 1H resonances bound to C6 of NAA and C4 of Gln to be revealed. The measured time course of the resolved labeling into NAA C6 with the new scheme was consistent with the slow turnover of NAA.