Multiple refugia and barriers explain the phylogeography of the Valais shrew, Sorex antinorii (Mammalia: Soricomorpha), Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)-based hindcasting approach. We also used landscape genetics (including isolation-by-resistance) to infer the determinants of current intra-specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median-joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north-eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi.

The second internal transcribed spacer of nuclear ribosomal DNA as a tool for Latin American anopheline taxonomy: a critical review

Among the molecular markers commonly used for mosquito taxonomy, the internal transcribed spacer 2 (ITS2) of the ribosomal DNA is useful for distinguishing among closely-related species. Here we review 178 GenBank accession numbers matching ITS2 sequences of Latin American anophelines. Among those, we found 105 unique sequences corresponding to 35 species. Overall the ITS2 sequences distinguish anopheline species, however, information on intraspecific and geographic variations is scarce. Intraspecific variations ranged from 0.2% to 19% and our analysis indicates that misidentification and/or sequencing errors could be responsible for some of the high values of divergence. Research in Latin American malaria vector taxonomy profited from molecular data provided by single or few field capture mosquitoes. However we propose that caution should be taken and minimum requirements considered in the design of additional studies. Future studies in this field should consider that: (1) voucher specimens, assigned to the DNA sequences, need to be deposited in collections, (2) intraspecific variations should be thoroughly evaluated, (3) ITS2 and other molecular markers, considered as a group, will provide more reliable information, (4) biological data about vector populations are missing and should be prioritized, (5) the molecular markers are most powerful when coupled with traditional taxonomic tools.

Representing diffusion MRI in 5D for segmentation of white matter tracts with a level set method.

We present a method for segmenting white matter tracts from high angular resolution diffusion MR. images by representing the data in a 5 dimensional space of position and orientation. Whereas crossing fiber tracts cannot be separated in 3D position space, they clearly disentangle in 5D position-orientation space. The segmentation is done using a 5D level set method applied to hyper-surfaces evolving in 5D position-orientation space. In this paper we present a methodology for constructing the position-orientation space. We then show how to implement the standard level set method in such a non-Euclidean high dimensional space. The level set theory is basically defined for N-dimensions but there are several practical implementation details to consider, such as mean curvature. Finally, we will show results from a synthetic model and a few preliminary results on real data of a human brain acquired by high angular resolution diffusion MRI.

Biological aspects of crosses between Triatoma maculata (Erichson, 1848) and Triatoma pseudomaculata Corrêa & Espínola, 1964 (Hemiptera: Reduviidae)

Due to morphological similarities between Triatoma maculata and T. pseudomaculata, which comprise the "maculate complex", both had been regarded as the same species until 1964. Considering that the studies on triatomine hybridization permit hypotheses formulation concerning origin and divergence of species, enabling a quantitative analysis of taxonomic relationships between species, the present investigation was aimed at broadening further understanding related to the capacity of hybrid production by determining the degree of reproductive isolation between T. maculata and T. pseudomaculata. Our results have demonstrated that T. maculata and T. pseudomaculata showed no differences regarding reproduction patterns and they are able to cross, generating infertile hybrids.

Genetic variability of Aedes aegypti in the Americas using a mitochondrial gene: evidence of multiple introductions

To analyze the genetic relatedness and phylogeographic structure of Aedes aegypti, we collected samples from 36 localities throughout the Americas (Brazil, Peru, Venezuela, Guatemala, US), three from Africa (Guinea, Senegal, Uganda), and three from Asia (Singapore, Cambodia, Tahiti). Amplification and sequencing of a fragment of the mitochondrial NADH dehydrogenase subunit 4 gene identified 20 distinct haplotypes, of which 14 are exclusive to the Americas, four to African/Asian countries, one is common to the Americas and Africa, and one to the Americas and Asia. Nested clade analysis (NCA), pairwise distribution, statistical parsimony, and maximum parsimony analyses were used to infer evolutionary and historic processes, and to estimate phylogenetic relationships among haplotypes. Two clusters were found in all the analyses. Haplotypes clustered in the two clades were separated by eight mutational steps. Phylogeographic structure detected by the NCA was consistent with distant colonization within one clade and fragmentation followed by range expansion via long distance dispersal in the other. Three percent of nucleotide divergence between these two clades is suggestive of a gene pool division that may support the hypothesis of occurrence of two subspecies of Ae. aegypti in the Americas.

Ambulatory gait analysis results after total knee arthroplasty : arguments for mobile or fixed bearing?

Background: Mobile-bearing knee replacements have some theoretical advantages over fixed-bearing devices. However, very few randomized controlled clinical trials have been published to date, and studies showed little clinical and subjective advantages for the mobile-bearing using traditional systems of scoring. The choice of the ideal outcome measure to assess total joint replacement remains a complex issue. However, gait analysis provides objective and quantifying evidences of treatment evaluation. Significant methodological advances are currently made in gait analysis laboratories and ambulatory gait devices are now available. The goal of this study was to provide gait parameters as a new objective method to assess total knee arthroplasty outcome between patients with fixed- and mobile-bearing, using an ambulatory device with minimal sensor configuration. Methods: This randomized controlled double-blind study included to date 31 patients: the gait signatures of 12 patients with mobile-bearing were compared to the gait signatures of 19 patients with fixed-bearing pre-operatively and post-operatively at 6 weeks, 3 months and 6 months. Each participant was asked to perform two walking trials of 30m long at his/her preferred speed and to complete a EQ-5D questionnaire, a WOMAC and Knee Society Score (KSS). Lower limbs rotations were measured by four miniature angular rate sensors mounted respectively, on each shank and thigh. Results: Better relative differences between pre-operative and post-operative 3 months and 6 months KSS (122% vs 34% at 3 months, 138% vs 36% at 6 months) and KSS function (154% vs 8% at 3 months, 183% vs 42% at 6 months) scores were observed for the fixed-bearing compared to the mobile-bearing. The same better improvements for fixed-bearing were also found with the range of knee angles (Affected side: 31% vs -5% at 3 months, 47% vs 5% at 6 months), (Unaffected side: 16% vs 5% at 3 months, 15% vs 6% at 6 months) and peak swing speeds of shank (Affected side: 18% vs -2% at 3 months, 30% vs 4% at 6 months), (Unaffected side: 8% vs -3% at 3 months, 7% vs 4% at 6 months). Conclusions: A new method for a portable system for gait analysis has been developed with very encouraging results regarding the objective outcome of total knee arthroplasty using mobile- and fixed-bearings.

Chromosome variability in the Chagas disease vector Rhodnius pallescens (Hemiptera, Reduviidae, Rhodniini)

Rhodnius pallescens is the main vector of Trypanosoma cruzi in Panama and one of the most relevant secondary vectors in Colombia. Despite the importance of this species, there is limited knowledge about the genetic variability along its geographical distribution. In order to evaluate the degree of karyotype variability we analyzed the meiotic behavior and banding pattern of the chromosomes of 112 males of R. pallescens coming from different regions of Colombia and Panama. Using the C-banding technique we identified two chromosomal patterns or cytotypes characterized by differences in the amount, size and distribution of constitutive heterochromatic regions in the chromosome complement (2n = 20 autosomes plus XY in males). The individuals can be easily classified in each cytotype by the analysis of the chromosomes during first meiotic prophase. The frequencies of the cytotypes are variable according to the geographic origin of the populations. This chromosomal divergence together with morphological data supports the existence of three genetically different populations of R. pallescens and provides new information to understand the distribution dynamics of this species.

Does a shift in host plants trigger speciation in the Alpine leaf beetle Oreina speciosissima (Coleoptera, Chrysomelidae)?

Background: Within the Coleoptera, the largest order in the animal kingdom, the exclusively herbivorous Chrysomelidae are recognized as one of the most species rich beetle families. The evolutionary processes that have fueled radiation into the more than thirty-five thousand currently recognized leaf beetle species remain partly unresolved. The prominent role of leaf beetles in the insect world, their omnipresence across all terrestrial biomes and their economic importance as common agricultural pest organisms make this family particularly interesting for studying the mechanisms that drive diversification. Here we specifically focus on two ecotypes of the alpine leaf beetle Oreina speciosissima (Scop.), which have been shown to exhibit morphological differences in male genitalia roughly corresponding to the subspecies Oreina speciosissima sensu stricto and Oreina speciosissima troglodytes. In general the two ecotypes segregate along an elevation gradient and by host plants: Oreina speciosissima sensu stricto colonizes high forb vegetation at low altitude and Oreina speciosissima troglodytes is found in stone run vegetation at higher elevations. Both host plants and leaf beetles have a patchy geographical distribution. Through use of gene sequencing and genome fingerprinting (AFLP) we analyzed the genetic structure and habitat use of Oreina speciosissima populations from the Swiss Alps to examine whether the two ecotypes have a genetic basis. By investigating a wide range of altitudes and focusing on the structuring effect of habitat types, we aim to provide answers regarding the factors that drive adaptive radiation in this phytophagous leaf beetle.Results: While little phylogenetic resolution was observed based on the sequencing of four DNA regions, the topology and clustering resulting from AFLP genotyping grouped specimens according to their habitat, mostly defined by plant associations. A few specimens with intermediate morphologies clustered with one of the two ecotypes or formed separate clusters consistent with habitat differences. These results were discussed in an ecological speciation framework.Conclusions: The question of whether this case of ecological differentiation occurred in sympatry or allopatry remains open. Still, the observed pattern points towards ongoing divergence between the two ecotypes which is likely driven by a recent shift in host plant use.

Brain network analysis of patients with Temporal Lobe Epilepsy

Patients with Temporal Lobe Epilepsy (TLE) suffer from widespread subtle white matter abnormalities and abnormal functional connectivity extending beyond the affected lobe, as revealed by Diffusion Tensor MR Imaging, volumetric and functional MRI studies. Diffusion Spectrum Imaging (DSI) is a diffusion imaging technique with high angular resolution for improving the mapping of white matter pathways. In this study, we used DSI, connectivity matrices and topological measures to investigate how the alteration in structural connectivity influences whole brain structural networks. Eleven patients with right-sided TLE and hippocampal sclerosis and 18 controls underwent our DSI protocol at 3T. The cortical and subcortical grey matters were parcellated into 86 regions of interest and the connectivity between every region pair was estimated using global tractography and a connectivity matrix (the adjacency matrix of the structural network). We then compared the networks of patients and controls using topological measures. In patients, we found a higher characteristic path length and a lower clustering coefficient compared to controls. Local measures at node level of the clustering and efficiency showed a significant difference after a multiple comparison correction (Bonferroni). These significant nodes were located within as well outside the temporal lobe, and the localisation of most of them was consistent with regions known to be part of epileptic networks in TLE. Our results show altered connectivity patterns that are concordant with the mapping of functional epileptic networks in patients with TLE. Further studies are needed to establish the relevance of these findings for the propagation of epileptic activity, cognitive deficits in medial TLE and outcome of epilepsy surgery in individual patients.

Ecological niche modelling and differentiation between Rhodnius neglectus Lent, 1954 and Rhodnius nasutus Stål, 1859 (Hemiptera: Reduviidae: Triatominae) in Brazil

Ecological niche modelling was used to predict the potential geographical distribution of Rhodnius nasutus Stål and Rhodnius neglectus Lent, in Brazil and to investigate the niche divergence between these morphologically similar triatomine species. The distribution of R. neglectus covered mainly the cerrado of Central Brazil, but the prediction maps also revealed its occurrence in transitional areas within the caatinga, Pantanal and Amazon biomes. The potential distribution of R. nasutus covered the Northeastern Region of Brazil in the semi-arid caatinga and the Maranhão babaçu forests. Clear ecological niche differences between these species were observed. R. nasutus occurred more in warmer and drier areas than R. neglectus. In the principal component analysis PC1 was correlated with altitude and temperature (mainly temperature in the coldest and driest months) and PC2 with vegetation index and precipitation. The prediction maps support potential areas of co-occurrence for these species in the Maranhão babaçu forests and in caatinga/cerrado transitional areas, mainly in state of Piaui. Entomologists engaged in Chagas disease vector surveillance should be aware that R. neglectus and R. nasutus can occur in the same localities of Northeastern Brazil. Thus, the identification of bugs in these areas should be improved by applying morphometrical and/or molecular methods.

How confident can we be that orthologs are similar, but paralogs differ?

Homologous genes are classified into orthologs and paralogs, depending on whether they arose by speciation or duplication. It is widely assumed that orthologs share similar functions, whereas paralogs are expected to diverge more from each other. But does this assumption hold up on further examination? We present evidence that orthologs and paralogs are not so different in either their evolutionary rates or their mechanisms of divergence. We emphasize the importance of appropriately designed studies to test models of gene evolution between orthologs and between paralogs. Thus, functional change between orthologs might be as common as between paralogs, and future studies should be designed to test the impact of duplication against this alternative model.

Compositional Data Analysis with R

R from http://www.r-project.org/ is ‘GNU S’ – a language and environment for statistical computingand graphics. The environment in which many classical and modern statistical techniques havebeen implemented, but many are supplied as packages. There are 8 standard packages and many moreare available through the cran family of Internet sites http://cran.r-project.org .We started to develop a library of functions in R to support the analysis of mixtures and our goal isa MixeR package for compositional data analysis that provides support foroperations on compositions: perturbation and power multiplication, subcomposition with or withoutresiduals, centering of the data, computing Aitchison’s, Euclidean, Bhattacharyya distances,compositional Kullback-Leibler divergence etc.graphical presentation of compositions in ternary diagrams and tetrahedrons with additional features:barycenter, geometric mean of the data set, the percentiles lines, marking and coloring ofsubsets of the data set, theirs geometric means, notation of individual data in the set . . .dealing with zeros and missing values in compositional data sets with R procedures for simpleand multiplicative replacement strategy,the time series analysis of compositional data.We’ll present the current status of MixeR development and illustrate its use on selected data sets

The molecular signature of the stroma response in prostate cancer-induced osteoblastic bone metastasis highlights expansion of hematopoietic and prostate epithelial stem cell niches.

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.

Evaluation of a mixed approach combining stationary and wearable systems to monitor gait over long distance.

Thanks to decades of research, gait analysis has become an efficient tool. However, mainly due to the price of the motion capture systems, standard gait laboratories have the capability to measure only a few consecutive steps of ground walking. Recently, wearable systems were proposed to measure human motion without volume limitation. Although accurate, these systems are incompatible with most of existing calibration procedures and several years of research will be necessary for their validation. A new approach consisting of using a stationary system with a small capture volume for the calibration procedure and then to measure gait using a wearable system could be very advantageous. It could benefit from the knowledge related to stationary systems, allow long distance monitoring and provide new descriptive parameters. The aim of this study was to demonstrate the potential of this approach. Thus, a combined system was proposed to measure the 3D lower body joints angles and segmental angular velocities. It was then assessed in terms of reliability towards the calibration procedure, repeatability and concurrent validity. The dispersion of the joint angles across calibrations was comparable to those of stationary systems and good reliability was obtained for the angular velocities. The repeatability results confirmed that mean cycle kinematics of long distance walks could be used for subjects' comparison and pointed out an interest for the variability between cycles. Finally, kinematics differences were observed between participants with different ankle conditions. In conclusion, this study demonstrated the potential of a mixed approach for human movement analysis.

Regional and hemispheric determinants of language laterality: Implications for preoperative fMRI.

Language is typically a function of the left hemisphere but the right hemisphere is also essential in some healthy individuals and patients. This inter-subject variability necessitates the localization of language function, at the individual level, prior to neurosurgical intervention. Such assessments are typically made by comparing left and right hemisphere language function to determine "language lateralization" using clinical tests or fMRI. Here, we show that language function needs to be assessed at the region and hemisphere specific level, because laterality measures can be misleading. Using fMRI data from 82 healthy participants, we investigated the degree to which activation for a semantic word matching task was lateralized in 50 different brain regions and across the entire cortex. This revealed two novel findings. First, the degree to which language is lateralized across brain regions and between subjects was primarily driven by differences in right hemisphere activation rather than differences in left hemisphere activation. Second, we found that healthy subjects who have relatively high left lateralization in the angular gyrus also have relatively low left lateralization in the ventral precentral gyrus. These findings illustrate spatial heterogeneity in language lateralization that is lost when global laterality measures are considered. It is likely that the complex spatial variability we observed in healthy controls is more exaggerated in patients with brain damage. We therefore highlight the importance of investigating within hemisphere regional variations in fMRI activation, prior to neuro-surgical intervention, to determine how each hemisphere and each region contributes to language processing. Hum Brain Mapp, 2010. © 2010 Wiley-Liss, Inc.