The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

Earthworm communities in organic and conventional coffee cultivation

The objective of this work was to evaluate the effect of organic and conventional coffee crops on biomass, population density and diversity of earthworms, in Lerroville, district of Londrina County, Paraná state, Brazil. Earthworm communities were sampled in three areas with organic coffee cultivation (CO1, CO2 and CO3), two with conventional coffee (CC1 and CC2), and a native forest fragment (MT). The soil of the areas CO1, CC1, and MT was classified as Nitossolo Vermelho (Rhodic Kandiudox), while CO2, CO3, and CC2 were on Latossolo Vermelho (Rhodic Hapludox). Eight samples were taken in each area on two occasions, winter and summer, using the Tropical Soil Biology and Fertility (TSBF) method in the 0-20 cm soil layer. The earthworms were handsorted and preserved in 4% formaldehyde, and were later weighed, counted and identified. The highest earthworm biomass, both in winter and summer, occurred in the CO3 area. For population density, the higher numbers of individuals were found in CO1 and CO3. The highest number of species was identified in the organic cultivation. The adoption of organic practices in coffee cultivation favored the diversity, density and biomass of earthworm communities.

UPLC-TOF-MS for plant metabolomics: a sequential approach for wound marker analysis in Arabidopsis thaliana.

The model plant Arabidopsis thaliana was studied for the search of new metabolites involved in wound signalling. Diverse LC approaches were considered in terms of efficiency and analysis time and a 7-min gradient on a UPLC-TOF-MS system with a short column was chosen for metabolite fingerprinting. This screening step was designed to allow the comparison of a high number of samples over a wide range of time points after stress induction in positive and negative ionisation modes. Thanks to data treatment, clear discrimination was obtained, providing lists of potential stress-induced ions. In a second step, the fingerprinting conditions were transferred to longer column, providing a higher peak capacity able to demonstrate the presence of isomers among the highlighted compounds.

Anti-angiogenic therapies in cancer: achievements and open questions.

The approval in 2004 of bevacizumab (Avastin), a neutralizing monoclonal antibody directed against vascular endothelial growth factor (VEGF) as the first anti-angiogenic systemic drug to treat cancer patients validated the notion introduced 33 years earlier by Dr. Judah Folkman, that inhibition of tumor angiogenesis might be a valid approach to control tumor growth. Anti-angiogenic therapy was greeted in the clinic a major step forward in cancer treatment. At the same time this success recently boosted the field to the quest for new anti-angiogenic targets and drugs. In spite of this success, however, some old questions in the field have remained unanswered and new ones have emerged. They include the identification for surrogate markers of angiogenesis and anti-angiogenesis, the understanding about how anti-angiogenic therapy and chemotherapy synergize, the characterization of the biological consequences of sustained suppression of angiogenesis on tumor biology and normal tissue homeostasis, and the mechanisms of tumor escape from anti-angiogenesis. In this review we summarize some of these outstanding questions, and highlight future challenges in clinical, translational and experimental research in anti-angiogenic therapy that need to be addressed in order to improve current treatments and to design new drugs.

Living with an imperfect cell wall: compensation of femAB inactivation in Staphylococcus aureus.

BACKGROUND: Synthesis of the Staphylococcus aureus peptidoglycan pentaglycine interpeptide bridge is catalyzed by the nonribosomal peptidyl transferases FemX, FemA and FemB. Inactivation of the femAB operon reduces the interpeptide to a monoglycine, leading to a poorly crosslinked peptidoglycan. femAB mutants show a reduced growth rate and are hypersusceptible to virtually all antibiotics, including methicillin, making FemAB a potential target to restore beta-lactam susceptibility in methicillin-resistant S. aureus (MRSA). Cis-complementation with wild type femAB only restores synthesis of the pentaglycine interpeptide and methicillin resistance, but the growth rate remains low. This study characterizes the adaptations that ensured survival of the cells after femAB inactivation. RESULTS: In addition to slow growth, the cis-complemented femAB mutant showed temperature sensitivity and a higher methicillin resistance than the wild type. Transcriptional profiling paired with reporter metabolite analysis revealed multiple changes in the global transcriptome. A number of transporters for sugars, glycerol, and glycine betaine, some of which could serve as osmoprotectants, were upregulated. Striking differences were found in the transcription of several genes involved in nitrogen metabolism and the arginine-deiminase pathway, an alternative for ATP production. In addition, microarray data indicated enhanced expression of virulence factors that correlated with premature expression of the global regulators sae, sarA, and agr. CONCLUSION: Survival under conditions preventing normal cell wall formation triggered complex adaptations that incurred a fitness cost, showing the remarkable flexibility of S. aureus to circumvent cell wall damage. Potential FemAB inhibitors would have to be used in combination with other antibiotics to prevent selection of resistant survivors.

Structure and function of RecA-DNA complexes.

While the E. coli RecA protein has been the most intensively studied enzyme of homologous recombination, the unusual RecA-DNA filament has stood alone until very recently. It now appears that this protein is part of a universal family that spans all of biology, and the filament that is formed by the protein on DNA is a universal structure. With RecA's role in recombination given new and greatly increased significance, we focus in this review on the energetics of the RecA-mediated strand exchange and the relation between the energetics and recombination spanning heterologous inserts.

Characterisation of the putative effector interaction site of the regulatory HbpR protein from Pseudomonas azelaica by site-directed mutagenesis.

Bacterial transcription activators of the XylR/DmpR subfamily exert their expression control via σ(54)-dependent RNA polymerase upon stimulation by a chemical effector, typically an aromatic compound. Where the chemical effector interacts with the transcription regulator protein to achieve activation is still largely unknown. Here we focus on the HbpR protein from Pseudomonas azelaica, which is a member of the XylR/DmpR subfamily and responds to biaromatic effectors such as 2-hydroxybiphenyl. We use protein structure modeling to predict folding of the effector recognition domain of HbpR and molecular docking to identify the region where 2-hydroxybiphenyl may interact with HbpR. A large number of site-directed HbpR mutants of residues in- and outside the predicted interaction area was created and their potential to induce reporter gene expression in Escherichia coli from the cognate P(C) promoter upon activation with 2-hydroxybiphenyl was studied. Mutant proteins were purified to study their conformation. Critical residues for effector stimulation indeed grouped near the predicted area, some of which are conserved among XylR/DmpR subfamily members in spite of displaying different effector specificities. This suggests that they are important for the process of effector activation, but not necessarily for effector specificity recognition.

Proliferative and Osteogenic Differentiation Potentials of Human Fetal Bone Cells

BACKGROUND. Human primary fetal bone cells (hFBC) are being characterized for use in bone tissue regeneration. Unlike human mesenchymal stem cells (hMSC), hFBC are partially differentiated with high expansion and regeneration potential. To date, proliferative and osteoblastic differentiation capacities of fetal bone cells remain poorly examined. The goal of this study was to define an environmental culture conditions for optimal proliferation and production of extracellular bone matrix leading to efficient bone repair. METHODS. Human primary FBC derived from our dedicated, consistent banks of bone cells comprising several fetal donors. For proliferation study, monolayer cultures of both cell types were expanded in DMEM or α- MEM media. Osteoblastic differentiation potentials of both hFBC and hMSC were evaluated through RT-PCR. Regulation of osteogenic differentiation by protein ligands Wnt3a and Wnt5a was studied by ALP enzymatic activity measurement. RESULTS. Evaluation of the proliferation rate demonstrated that hFBC proliferated more rapidly in α-MEM medium. Regarding growth factors that could stimulate cell proliferation rate, we observed that PDGF, FGF2 and Wnt3a had positive effects on proliferation of hFBC. Gene expression analysis demonstrated a higher expression of runx2 in hFBC cultured in basal conditions, which was was similar than that was observed in hMSC in osteoinductive culture conditions. Expression of sox9 was very low in hBFC and hMSC, compared to expression observed in fetal cartilage cells. Looking at osteogenic differentiation capacity, ALP activity was positively regulated byWnt5awhen hFBCwere cultured inα-MEM, but not in DMEM. Conversely, Wnt3a was shown to block the effect of osteogenic inductors on differentiation of both cell types. CONCLUSION. Data presented in this study indicate that the proliferation and differentiation of fetal and mesenchymal stem cells is optimal in α- MEM. Evidence for a pre-differentiated state of hBFC was given by extracellular matrix spontaneous mineralization as well as by higher ALP activity levels observed for these cells in baseline culture conditions, in comparison with hMSC. As we showed that, in vitro, hFBC express a higher capacity to differentiate in osteoblasts, they represent an attractive and promising prospect for fundamental research, and specifically for a new generation of skeletal tissue engineering.

Identification of cancer/testis-antigen genes by massively parallel signature sequencing.

Massively parallel signature sequencing (MPSS) generates millions of short sequence tags corresponding to transcripts from a single RNA preparation. Most MPSS tags can be unambiguously assigned to genes, thereby generating a comprehensive expression profile of the tissue of origin. From the comparison of MPSS data from 32 normal human tissues, we identified 1,056 genes that are predominantly expressed in the testis. Further evaluation by using MPSS tags from cancer cell lines and EST data from a wide variety of tumors identified 202 of these genes as candidates for encoding cancer/testis (CT) antigens. Of these genes, the expression in normal tissues was assessed by RT-PCR in a subset of 166 intron-containing genes, and those with confirmed testis-predominant expression were further evaluated for their expression in 21 cancer cell lines. Thus, 20 CT or CT-like genes were identified, with several exhibiting expression in five or more of the cancer cell lines examined. One of these genes is a member of a CT gene family that we designated as CT45. The CT45 family comprises six highly similar (>98% cDNA identity) genes that are clustered in tandem within a 125-kb region on Xq26.3. CT45 was found to be frequently expressed in both cancer cell lines and lung cancer specimens. Thus, MPSS analysis has resulted in a significant extension of our knowledge of CT antigens, leading to the discovery of a distinctive X-linked CT-antigen gene family.

Implicit methods for qualitative modeling of gene regulatory networks.

Advancements in high-throughput technologies to measure increasingly complex biological phenomena at the genomic level are rapidly changing the face of biological research from the single-gene single-protein experimental approach to studying the behavior of a gene in the context of the entire genome (and proteome). This shift in research methodologies has resulted in a new field of network biology that deals with modeling cellular behavior in terms of network structures such as signaling pathways and gene regulatory networks. In these networks, different biological entities such as genes, proteins, and metabolites interact with each other, giving rise to a dynamical system. Even though there exists a mature field of dynamical systems theory to model such network structures, some technical challenges are unique to biology such as the inability to measure precise kinetic information on gene-gene or gene-protein interactions and the need to model increasingly large networks comprising thousands of nodes. These challenges have renewed interest in developing new computational techniques for modeling complex biological systems. This chapter presents a modeling framework based on Boolean algebra and finite-state machines that are reminiscent of the approach used for digital circuit synthesis and simulation in the field of very-large-scale integration (VLSI). The proposed formalism enables a common mathematical framework to develop computational techniques for modeling different aspects of the regulatory networks such as steady-state behavior, stochasticity, and gene perturbation experiments.

Cannabidiol protects against hepatic ischemia/reperfusion injury by attenuating inflammatory signaling and response, oxidative/nitrative stress, and cell death.

Ischemia/reperfusion (I/R) is a pivotal mechanism of liver damage after liver transplantation or hepatic surgery. We have investigated the effects of cannabidiol (CBD), the nonpsychotropic constituent of marijuana, in a mouse model of hepatic I/R injury. I/R triggered time-dependent increases/changes in markers of liver injury (serum transaminases), hepatic oxidative/nitrative stress (4-hydroxy-2-nonenal, nitrotyrosine content/staining, and gp91phox and inducible nitric oxide synthase mRNA), mitochondrial dysfunction (decreased complex I activity), inflammation (tumor necrosis factor α (TNF-α), cyclooxygenase 2, macrophage inflammatory protein-1α/2, intercellular adhesion molecule 1 mRNA levels; tissue neutrophil infiltration; nuclear factor κB (NF-κB) activation), stress signaling (p38MAPK and JNK), and cell death (DNA fragmentation, PARP activity, and TUNEL). CBD significantly reduced the extent of liver inflammation, oxidative/nitrative stress, and cell death and also attenuated the bacterial endotoxin-triggered NF-κB activation and TNF-α production in isolated Kupffer cells, likewise the adhesion molecule expression in primary human liver sinusoidal endothelial cells stimulated with TNF-α and attachment of human neutrophils to the activated endothelium. These protective effects were preserved in CB(2) knockout mice and were not prevented by CB(1/2) antagonists in vitro. Thus, CBD may represent a novel, protective strategy against I/R injury by attenuating key inflammatory pathways and oxidative/nitrative tissue injury, independent of classical CB(1/2) receptors.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ∼2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

GENCODE: producing a reference annotation for ENCODE.

BACKGROUND: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manual annotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results. RESULTS: The GENCODE gene features are divided into eight different categories of which only the first two (known and novel coding sequence) are confidently predicted to be protein-coding genes. 5' rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentally verify the initial annotation. Of the 420 coding loci tested, 229 RACE products have been sequenced. They supported 5' extensions of 30 loci and new splice variants in 50 loci. In addition, 46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15 putative transcripts. We assessed the comprehensiveness of the GENCODE annotation by attempting to validate all the predicted exon boundaries outside the GENCODE annotation. Out of 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only two of them in intergenic regions. CONCLUSION: In total, 487 loci, of which 434 are coding, have been annotated as part of the GENCODE reference set available from the UCSC browser. Comparison of GENCODE annotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained within the two sets, which is a reflection of the high number of alternative splice forms with unique exons annotated. Over 50% of coding loci have been experimentally verified by 5' RACE for EGASP and the GENCODE collaboration is continuing to refine its annotation of 1% human genome with the aid of experimental validation.

The origins, evolution, and functional potential of alternative splicing in vertebrates.

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.

Body image perception and weight-related behaviour among adolescents of the Seychelles

Background: We examined one's own body image perception and its association with reported weight-related behavior among adolescents of a rapidly developing country in the African region. Methods: We conducted a school-based survey of 1432 students aged 11-17 years in the Seychelles. Weight and height were measured, and thinness, normal weight and overweight were assessed along standard criteria. A self-administered and anonymous questionnaire was administered. Perception of body image was assessed using both a closed-ended question (CEQ) and the Stunkard's pictorial silhouettes (SPS). Finally, a question assessed voluntary attempts to change weight. Results: Overall, 14.1% of the students were thin, 63.9% were normal-weight, and 22.0% were overweight or obese. There was fair agreement between actual weight status and self-perceived body image based on either CEQ or SPS. However, a substantial proportion of the overweight students did not consider themselves as overweight (SPS: 24%, CEQ: 34%) and, inversely, a substantial proportion of the normal-weight students considered themselves as too thin (SPS: 29%, CEQ: 15%). Among the overweight students, an adequate attempt to lose weight was reported more often by boys and girls who perceived themselves as overweight vs. not overweight (72-88% vs. 40-71%, p <0.05 for most comparisons). Among the normal-weight students, an inadequate attempt to gain weight was reported more often by boys and girls who perceived themselves as thin vs. not thin (27-68% vs. 11-19%, p <0.05). Girls had leaner own body ideals than boys. Conclusions: We found that substantial proportions of overweight students did not perceive themselves as overweight and/or did not want to lose weight and, inversely, that many normalweight students perceived themselves as too thin and/or wanted to gain weight: this points to forces that can drive the upwards overweight trends. Appropriate perception of one's weight was associated with adequate weight-control behavior, although not strongly, emphasizing that appropriate weight perception is only one of several factors driving adequate weight-related behavior. These findings emphasize the need to address appropriate perception of one's own weight and adequate weight-related behavior in adolescents for both individual and community weight-related interventions.