Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct monte carlo methods.

Dose kernel convolution (DK) methods have been proposed to speed up absorbed dose calculations in molecular radionuclide therapy. Our aim was to evaluate the impact of tissue density heterogeneities (TDH) on dosimetry when using a DK method and to propose a simple density-correction method. METHODS: This study has been conducted on 3 clinical cases: case 1, non-Hodgkin lymphoma treated with (131)I-tositumomab; case 2, a neuroendocrine tumor treatment simulated with (177)Lu-peptides; and case 3, hepatocellular carcinoma treated with (90)Y-microspheres. Absorbed dose calculations were performed using a direct Monte Carlo approach accounting for TDH (3D-RD), and a DK approach (VoxelDose, or VD). For each individual voxel, the VD absorbed dose, D(VD), calculated assuming uniform density, was corrected for density, giving D(VDd). The average 3D-RD absorbed dose values, D(3DRD), were compared with D(VD) and D(VDd), using the relative difference Δ(VD/3DRD). At the voxel level, density-binned Δ(VD/3DRD) and Δ(VDd/3DRD) were plotted against ρ and fitted with a linear regression. RESULTS: The D(VD) calculations showed a good agreement with D(3DRD). Δ(VD/3DRD) was less than 3.5%, except for the tumor of case 1 (5.9%) and the renal cortex of case 2 (5.6%). At the voxel level, the Δ(VD/3DRD) range was 0%-14% for cases 1 and 2, and -3% to 7% for case 3. All 3 cases showed a linear relationship between voxel bin-averaged Δ(VD/3DRD) and density, ρ: case 1 (Δ = -0.56ρ + 0.62, R(2) = 0.93), case 2 (Δ = -0.91ρ + 0.96, R(2) = 0.99), and case 3 (Δ = -0.69ρ + 0.72, R(2) = 0.91). The density correction improved the agreement of the DK method with the Monte Carlo approach (Δ(VDd/3DRD) < 1.1%), but with a lesser extent for the tumor of case 1 (3.1%). At the voxel level, the Δ(VDd/3DRD) range decreased for the 3 clinical cases (case 1, -1% to 4%; case 2, -0.5% to 1.5%, and -1.5% to 2%). No more linear regression existed for cases 2 and 3, contrary to case 1 (Δ = 0.41ρ - 0.38, R(2) = 0.88) although the slope in case 1 was less pronounced. CONCLUSION: This study shows a small influence of TDH in the abdominal region for 3 representative clinical cases. A simple density-correction method was proposed and improved the comparison in the absorbed dose calculations when using our voxel S value implementation.

Interrogating eleven fast-evolving genes for signatures of recent positive selection in worldwide human populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs inchimpanzee, mouse, rat and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes were then further explored for signatures of recent positive selection using SNP data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP Diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintainance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal ofrecent selection in this region may come from a different, neighbouring gene (CD5) ratherthan from the candidate gene itself (VPS37C). For this set of positively-selected genes in thehuman lineage, we find no indication that these genes maintained their rapid evolutionarypace among human populations. Based on these data, it therefore appears that adaptation forhuman-specific and for population-specific traits may have involved different genes.

Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment and innate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints. Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection are frequent on genes that are known to be at bifurcation points, and that are identified as being in key position by a network-level analysis such as MGAT3 and GCS1. These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

Recognition of Gram-positive Intestinal Bacteria by Hybridoma- and Colostrum-derived Secretory Immunoglobulin A Is Mediated by Carbohydrates.

Humans live in symbiosis with 10(14) commensal bacteria among which >99% resides in their gastrointestinal tract. The molecular bases pertaining to the interaction between mucosal secretory IgA (SIgA) and bacteria residing in the intestine are not known. Previous studies have demonstrated that commensals are naturally coated by SIgA in the gut lumen. Thus, understanding how natural SIgA interacts with commensal bacteria can provide new clues on its multiple functions at mucosal surfaces. Using fluorescently labeled, nonspecific SIgA or secretory component (SC), we visualized by confocal microscopy the interaction with various commensal bacteria, including Lactobacillus, Bifidobacteria, Escherichia coli, and Bacteroides strains. These experiments revealed that the interaction between SIgA and commensal bacteria involves Fab- and Fc-independent structural motifs, featuring SC as a crucial partner. Removal of glycans present on free SC or bound in SIgA resulted in a drastic drop in the interaction with Gram-positive bacteria, indicating the essential role of carbohydrates in the process. In contrast, poor binding of Gram-positive bacteria by control IgG was observed. The interaction with Gram-negative bacteria was preserved whatever the molecular form of protein partner used, suggesting the involvement of different binding motifs. Purified SIgA and SC from either mouse hybridoma cells or human colostrum exhibited identical patterns of recognition for Gram-positive bacteria, emphasizing conserved plasticity between species. Thus, sugar-mediated binding of commensals by SIgA highlights the currently underappreciated role of glycans in mediating the interaction between a highly diverse microbiota and the mucosal immune system.

Interpreting Breast International Group (BIG) 1-98: a randomized, double-blind, phase III trial comparing letrozole and tamoxifen as adjuvant endocrine therapy for postmenopausal women with hormone receptor-positive, early breast cancer.

The Breast International Group (BIG) 1-98 study is a four-arm trial comparing 5 years of monotherapy with tamoxifen or with letrozole or with sequences of 2 years of one followed by 3 years of the other for postmenopausal women with endocrine-responsive early invasive breast cancer. From 1998 to 2003, BIG -98 enrolled 8,010 women. The enhanced design f the trial enabled two complementary analyses of efficacy and safety. Collection of tumor specimens further enabled treatment comparisons based on tumor biology. Reports of BIG 1-98 should be interpreted in relation to each individual patient as she weighs the costs and benefits of available treatments. Clinicaltrials.gov ID: NCT00004205.

The complete genome sequence of the gram-positive bacterium Bacillus subtilis.

Bacillus subtilis is the best-characterized member of the Gram-positive bacteria. Its genome of 4,214,810 base pairs comprises 4,100 protein-coding genes. Of these protein-coding genes, 53% are represented once, while a quarter of the genome corresponds to several gene families that have been greatly expanded by gene duplication, the largest family containing 77 putative ATP-binding transport proteins. In addition, a large proportion of the genetic capacity is devoted to the utilization of a variety of carbon sources, including many plant-derived molecules. The identification of five signal peptidase genes, as well as several genes for components of the secretion apparatus, is important given the capacity of Bacillus strains to secrete large amounts of industrially important enzymes. Many of the genes are involved in the synthesis of secondary metabolites, including antibiotics, that are more typically associated with Streptomyces species. The genome contains at least ten prophages or remnants of prophages, indicating that bacteriophage infection has played an important evolutionary role in horizontal gene transfer, in particular in the propagation of bacterial pathogenesis.

Positive time-frequency distributions based on joint marginal constraints

This correspondence studies the formulation of members ofthe Cohen-Posch class of positive time-frequency energy distributions.Minimization of cross-entropy measures with respect to different priorsand the case of no prior or maximum entropy were considered. It isconcluded that, in general, the information provided by the classicalmarginal constraints is very limited, and thus, the final distributionheavily depends on the prior distribution. To overcome this limitation,joint time and frequency marginals are derived based on a "directioninvariance" criterion on the time-frequency plane that are directly relatedto the fractional Fourier transform.

Can resistance to trastuzumab be reversed by endocrine therapy? Combination of trastuzumab and letrozole after resistance to sequential trastuzumab and aromatase inhibitor monotherapies in patients with ER-positive, HER-2 positive, advanced breast cancer: a proof-of-concept trial (SAKK 23/03)

Introduction: Trastuzumab (T) is a cornerstone in the treatment of patients with HER2-overexpressing advanced breast cancer and development of resistance to T is a major therapeutic problem. HER-2 is part of a highly interactive signaling network that may impair efficacy of endocrine therapy. A sequential treatment design was chosen in this trial to ensure complete resistance to single agent therapy before receiving both a non-steroidal aromatase inhibitor (AI) and T. Any kind of clinical activity with combined treatment of AI and T after progression of single agent treatments could indicate restoration of sensitivity as a consequence of cross-talking and networking between both pathways. Methods: Key eligibility criteria included postmenopausal patients (pts.) with advanced, measurable, HER-2 positive (assessed by FISH, ratio (≥2)), HR positive disease and progression on prior treatment with a non-steroidal AI, e.g. letrozole or anastrozole, either in an adjuvant or advanced setting. Pts. received standard dose T monotherapy either weekly or three-weekly in step 1 and upon disease progression, continued T in combination with letrozole in step 2. The primary endpoint was clinical benefit response (CBR: CR, PR or SD for at least 24 weeks (+/- 1 week) according to RECIST) in step 2. Results: Thirteen pts. were enrolled in five centers in Switzerland. In step 1, six pts. (46%) achieved CBR. Median time to progression (TTP) was 161 days (Range: 50 - 627). Based on data collected until the end of May 2010, CBR was observed in seven out of the eleven evaluable pts. (64%) in step 2, including one pt. with partial response. Four of the seven pts. within step 2 that achieved CBR also had CBR in step 1. Seven out of eleven pts. have documented tumor progression during step 2 treatment. Median TTP for all eleven pts. was 184 days (range 61 - 471). Mean time on study treatment (TTP in step 1 plus TTP in step 2) for pts. reaching step 2 was 380 days (range 174 - 864). Adverse events were generally mild. Conclusion: Results of this proof-of-principle trial suggest that complete resistance to both AI and T can be overcome in a proportion of pts. by combined treatment of AI and T, as all pts. served as their own control. Our results appear promising for a new treatment strategy which offers a chemotherapy-free and well-tolerated option for at least a subset of the pts. with HR positive, HER-2 positive breast cancer. Further trials will need to corroborate this finding.

Rifaximin treatment for the irritable bowel syndrome with a positive lactulose hydrogen breath test improves symptoms for at least 3 months.

BACKGROUND: While rifaximin was able to improve symptoms in patients with irritable bowel syndrome (IBS) in phase III trials, these results are yet to be repeated in phase IV studies. AIM: To evaluate the treatment response to rifaximin in IBS patients in a phase IV trial. METHODS: IBS patients underwent lactulose hydrogen breath testing (LHBT). LHBT-positive patients were treated with rifaximin for 14 days. Prior to treatment as well as at week 4 and 14 following the start of rifaximin treatment, patients completed a questionnaire assessing symptom severity on a Likert scale from 0 to 10. RESULTS: One hundred and six of 150 IBS patients (71%) were LHBT-positive and treated with rifaximin. As assessed at week 4 following commencement of the therapy, rifaximin provided significant improvement of the following IBS-associated symptoms: bloating (5.5±2.6 before the start of the treatment vs. 3.6±2.7 at week 4, P<0.001), flatulence (5.0±2.7 vs. 4.0±2.7, P=0.015), diarrhoea (2.9±2.4 vs. 2.0±2.4, P=0.005) and abdominal pain (4.8±2.7 vs. 3.3±2.5, P<0.001). Overall well-being also significantly improved (3.9 ± 2.4 vs. 2.7 ± 2.3, P < 0.001). Similar improvements in IBS symptoms were obtained at week 14. Eighty-six per cent of patients undergoing repetitive LHBT (55/64) tested negative at week 4. CONCLUSIONS: We found a high percentage of LHBT-positive IBS patients. IBS-associated symptoms (bloating, flatulence, diarrhoea, pain) were improved for a period of 3 months following 2 weeks of treatment with rifaximin. We conclude that rifaximin treatment alleviates symptoms in LHBT-positive IBS patients.

Predicting prognosis using molecular profiling in estrogen receptor-positive breast cancer treated with tamoxifen.

Background: Estrogen receptor positive (ER+) breast cancers (BC) are heterogeneous with regard to their clinical behavior and response to therapies. The ER is currently the best predictor of response to the anti-estrogen agent tamoxifen, yet up to 30-40% of ER+ BC will relapse despite tamoxifen treatment. New prognostic biomarkers and further biological understanding of tamoxifen resistance are required. We used gene expression profiling to develop an outcome-based predictor using a training set of 255 ER+ BC samples from women treated with adjuvant tamoxifen monotherapy. We used clusters of highly correlated genes to develop our predictor to facilitate both signature stability and biological interpretation. Independent validation was performed using 362 tamoxifen-treated ER+ BC samples obtained from multiple institutions and treated with tamoxifen only in the adjuvant and metastatic settings.Results: We developed a gene classifier consisting of 181 genes belonging to 13 biological clusters. In the independent set of adjuvantly-treated samples, it was able to define two distinct prognostic groups (HR 2.01 95% CI: 1.29-3.13; p = 0.002). Six of the 13 gene clusters represented pathways involved in cell cycle and proliferation. In 112 metastatic breast cancer patients treated with tamoxifen, one of the classifier components suggesting a cellular inflammatory mechanism was significantly predictive of response.Conclusion: We have developed a gene classifier that can predict clinical outcome in tamoxifen-treated ER+ BC patients. Whilst our study emphasizes the important role of proliferation genes in prognosis, our approach proposes other genes and pathways that may elucidate further mechanisms that influence clinical outcome and prediction of response to tamoxifen.

The potential of metastatic lymph nodes to predict positive resection margins in patients with resectable pancreatic cancer

Objective: Despite progress during recent decades, long-term outcome¦of patients with pancreatic cancer remains dismal. Since positive resection¦margins and metastatic lymph nodes are known risk factors for early tumor¦recurrence, patients at risk should be identified and could potentially benefit¦from preoperative radio-chemotherapy. This study aimed to assess whether the¦presence of lymph node metastasis could be used to predict positive resection¦margins in patients with pancreatic cancer.¦Methods: A series of 146 patients (82 male, 64 female, median age 68 years)¦underwent pancreatic head resection for various malignant diseases (pancreatic¦ductal adenocarcinoma, biliary cancer, periampullary cancer) at our institution¦from 2000 to 2011. Patients were identified from our prospective database¦that collects more than 60 single items of all patients undergoing pancreatic¦resection. Lymph node metastasis and positive resection margins were all¦confirmed by histological evaluation. Positive predictive value (PPV), negative¦predictive value (NPV), sensitivity and specificity were calculated to assess¦the predictive value of metastatic lymph nodes regarding tumor-free (R0) and¦tumor-involved (R1) resection margins.¦Results: There were 110 specimens (76%) with tumor-positive lymph nodes¦and 36 specimens with tumor-negative lymph nodes. Resection margins were¦positive in 47 specimens (32%) and negative in 99 specimens. Sensitivity of¦tumor-positive lymph nodes to detect positive resection margins was 96%, and¦the NPV was 94%. In contrast, specificity was 34% and the PPV was 41%.¦Conclusion: Patients with resectable pancreatic cancer, who have no lymph¦node metastasis, are at very low risk to have positive resection margins (2 of¦36 patients, NPV 94%). In contrast, more than one third of patients with¦metastatic lymph nodes are at increased risk for an incomplete tumor resection¦(sensitivity 96%). If lymph nodesmetastases are highly suspected at preoperative¦staging, a neoadjuvant treatment strategy should be considered to increase the¦R0 resection rate.

Coping with an HIV infection. A multicenter qualitative survey on HIV positive adolescents' perceptions of their disease, therapeutic adherence and treatment.

HIV-positive adolescents face a number of challenges in dealing with their disease and its treatment. In this qualitative study, twenty-nine HIV-positive adolescents aged 13 to 20 years (22 girls), who live in Switzerland, were asked, in a semi-structured interview (duration of 40-110 minutes), to describe their perceptions and experiences with the disease itself and with therapeutic adherence. While younger adolescents most often thought of their disease as fate, older adolescents usually knew that they had received it through vertical transmission, although the topic appeared to be particularly difficult to discuss for those living with their HIV-positive mothers. Based on their attending physician's assessment, 18 subjects were judged highly adherent, 4 fairly and 7 poorly adherent. High adherence appeared linked with adequate psychological adjustment and effective coping mechanisms, as well as with the discussion and adoption of explicit medication-taking strategies. The setting and organisation of health care teams should allow for ongoing discussions with HIV-positive adolescents that focus on their perceptions of their disease, how they cope with it and with the treatment, and how they could improve their adherence.

Continuous positive airway pressure causes lung injury in a model of sepsis.

Continuous positive airway pressure, aimed at preventing pulmonary atelectasis, has been used for decades to reduce lung injury in critically ill patients. In neonatal practice, it is increasingly used worldwide as a primary form of respiratory support due to its low cost and because it reduces the need for endotracheal intubation and conventional mechanical ventilation. We studied the anesthetized in vivo rat and determined the optimal circuit design for delivery of continuous positive airway pressure. We investigated the effects of continuous positive airway pressure following lipopolysaccharide administration in the anesthetized rat. Whereas neither continuous positive airway pressure nor lipopolysaccharide alone caused lung injury, continuous positive airway pressure applied following intravenous lipopolysaccharide resulted in increased microvascular permeability, elevated cytokine protein and mRNA production, and impaired static compliance. A dose-response relationship was demonstrated whereby higher levels of continuous positive airway pressure (up to 6 cmH(2)O) caused greater lung injury. Lung injury was attenuated by pretreatment with dexamethasone. These data demonstrate that despite optimal circuit design, continuous positive airway pressure causes significant lung injury (proportional to the airway pressure) in the setting of circulating lipopolysaccharide. Although we would currently avoid direct extrapolation of these findings to clinical practice, we believe that in the context of increasing clinical use, these data are grounds for concern and warrant further investigation.

Differentiation of blue ballpoint pen inks by positive and negative mode LDI-MS

Usually, the differentiation of inks on questioned documents is carried out by optical methods and thin layer chromatography (TLC). Therefore, spectrometric methods were also proposed in forensic literature for the analysis of dyes. Between these techniques, laser desorption/ionization mass spectrometry (LDI-MS) has demonstrated a great versatility thanks to its sensitivity to blue ballpoint ink dyes and minimal sample destruction. Previous researches concentrated mostly on the LDI-MS positive mode and have shown that this analytical tool offers higher discrimination power than high performance TLC (HPTLC) for the differentiation of blue ballpoint inks. Although LDI-MS negative mode has already been applied in numerous forensic domains like the studies of works of art, automotive paints or rollerball pens, its potential for the discrimination of ballpoint pens was never studied before. The aim of the present paper is therefore to evaluate its potential for the discrimination of blue ballpoint inks. After optimization of the method, ink entries from 33 blue ballpoint pens were analyzed directly on paper in both positive and negative modes by LDI-MS. Several cationic and anionic ink components were identified in inks; therefore, pens were classified and compared according to their formulations. Results show that additional information provided by anionic dyes and pigments significantly increases the discrimination power of positive mode. In fact, it was demonstrated that classifications obtained by the two modes were, to some extent, complementary (i.e., inks with specific cationic dyes not necessarily contained the same anionic components).