Bioinformatics: A challenge for statisticians

Bioinformatics is a subject that requires the skills of biologists, computer scientists, mathematicians and staisticians. This paper introduces the reader to one small aspect of the subject: the study of microarrays. It describes some of the complexities of the enormous amounts of data that are available and shows how simple statistical techniques can be used to highlight deficiencies in that data

Motif effects in Affymetrix GeneChips seriously affect probe intensities

An Affymetrix GeneChip consists of an array of hundreds of thousands of probes (each a sequence of 25 bases) with the probe values being used to infer the extent to which genes are expressed in the biological material under investigation. In this article, we demonstrate that these probe values are also strongly influenced by their precise base sequence. We use data from >28 000 CEL files relating to 10 different Affymetrix GeneChip platforms and involving nearly 1000 experiments. Our results confirm known effects (those due to the T7-primer and the formation of G-quadruplexes) but reveal other effects. We show that there can be huge variations from one experiment to another, and that there may also be sizeable disparities between batches within an experiment and between CEL files within a batch. © 2012 The Author(s)

Projected expansion of Trichodesmium’s geographical distribution and increase of growth potential in response to climate change

Estimates of marine N₂ fixation range from 52 to 73 Tg N yr‾¹, of which we calculate up to 84% is from Trichodesmium based on previous measurements of nifH gene abundance and our new model of Trichodesmium growth. Here we assess the likely effects of four major climate change‐related abiotic factors on the spatiotemporal distribution and growth potential of Trichodesmium for the last glacial maximum (LGM), the present (2006‐2015) and the end of this century (2100) by mapping our model of Trichodesmium growth onto inferred global surface ocean fields of pCO₂, temperature, light and Fe. We conclude that growth rate was severely limited by low pCO₂ at the LGM, that current pCO₂ levels do not significantly limit Trichodesmium growth and thus, the potential for enhanced growth from future increases of CO₂ is small. We also found that the area of the ocean where sea surface temperatures (SST) are within Trichodesmium’s thermal niche increased by 32% from the LGM to present, but further increases in SST due to continued global warming will reduce this area by 9%. However, the range reduction at the equator is likely to be offset by enhanced growth associated with expansion of regions with optimal or near optimal Fe and light availability. Between now and 2100, the ocean area of optimal SST and irradiance is projected to increase by 7%, and the ocean area of optimal SST, irradiance and iron is projected to increase by 173%. Given the major contribution of this keystone species to annual N₂ fixation and thus pelagic ecology, biogeochemistry and CO₂ sequestration, the projected increase in the geographical range for optimal growth could provide a negative feedback to increasing atmospheric CO₂ concentrations

Structure and functional motifs of GCR1, the only plant protein with a GPCR fold?

Whether GPCRs exist in plants is a fundamental biological question. Interest in deorphanizing new G protein coupled receptors (GPCRs), arises because of their importance in signaling. Within plants, this is controversial as genome analysis has identified 56 putative GPCRs, including GCR1 which is reportedly a remote homologue to class A, B and E GPCRs. Of these, GCR2, is not a GPCR; more recently it has been proposed that none are, not even GCR1. We have addressed this disparity between genome analysis and biological evidence through a structural bioinformatics study, involving fold recognition methods, from which only GCR1 emerges as a strong candidate. To further probe GCR1, we have developed a novel helix alignment method, which has been benchmarked against the the class A – class B - class F GPCR alignments. In addition, we have presented a mutually consistent set of alignments of GCR1 homologues to class A, class B and class F GPCRs, and shown that GCR1 is closer to class A and /or class B GPCRs than class A, class B or class F GPCRs are to each other. To further probe GCR1, we have aligned transmembrane helix 3 of GCR1 to each of the 6 GPCR classes. Variability comparisons provide additional evidence that GCR1 homologues have the GPCR fold. From the alignments and a GCR1 comparative model we have identified motifs that are common to GCR1, class A, B and E GPCRs. We discuss the possibilities that emerge from this controversial evidence that GCR1 has a GPCR fol

Assessing the effect of dynamics on the closed-loop protein-folding hypothesis

The closed-loop (loop-n-lock) hypothesis of protein folding suggests that loops of about 25 residues, closed through interactions between the loop ends (locks), play an important role in protein structure. Coarse-grain elastic network simulations, and examination of loop lengths in a diverse set of proteins, each supports a bias towards loops of close to 25 residues in length between residues of high stability. Previous studies have established a correlation between total contact distance (TCD), a metric of sequence distances between contacting residues (cf. contact order), and the log-folding rate of a protein. In a set of 43 proteins, we identify an improved correlation ( r 2 = 0.76), when the metric is restricted to residues contacting the locks, compared to the equivalent result when all residues are considered ( r 2 = 0.65). This provides qualified support for the hypothesis, albeit with an increased emphasis upon the importance of a much larger set of residues surrounding the locks. Evidence of a similar-sized protein core/extended nucleus (with significant overlap) was obtained from TCD calculations in which residues were successively eliminated according to their hydrophobicity and connectivity, and from molecular dynamics simulations. Our results suggest that while folding is determined by a subset of residues that can be predicted by application of the closed-loop hypothesis, the original hypothesis is too simplistic; efficient protein folding is dependent on a considerably larger subset of residues than those involved in lock formation. </jats:p

Normalized Affymetrix expression data are biased by G-quadruplex formation

Probes with runs of four or more guanines (G-stacks) in their sequences can exhibit a level of hybridization that is unrelated to the expression levels of the mRNA that they are intended to measure. This is most likely caused by the formation of G-quadruplexes, where inter-probe guanines form Hoogsteen hydrogen bonds, which probes with G-stacks are capable of forming. We demonstrate that for a specific microarray data set using the Human HG-U133A Affymetrix GeneChip and RMA normalization there is significant bias in the expression levels, the fold change and the correlations between expression levels. These effects grow more pronounced as the number of G-stack probes in a probe set increases. Approximately 14 of the probe sets are directly affected. The analysis was repeated for a number of other normalization pipelines and two, FARMS and PLIER, minimized the bias to some extent. We estimate that ∼15 of the data sets deposited in the GEO database are susceptible to the effect. The inclusion of G-stack probes in the affected data sets can bias key parameters used in the selection and clustering of genes. The elimination of these probes from any analysis in such affected data sets outweighs the increase of noise in the signal. © 2011 The Author(s)

Which circulating antioxidant vitamins are confounded by socioeconomic deprivation? The MIDSPAN family study

&lt;p&gt;&lt;b&gt;Background:&lt;/b&gt; Antioxidant vitamins are often described as having “independent” associations with risk of cancer, cardiovascular disease (CVD) and mortality. We aimed to compare to what extent a range of antioxidant vitamins and carotenoids are associated with adulthood and childhood markers of socioeconomic deprivation and to adverse lifestyle factors.&lt;/p&gt; &lt;p&gt;&lt;b&gt;Methods and Findings:&lt;/b&gt; Socioeconomic and lifestyle measures were available in 1040 men and 1298 women from the MIDSPAN Family Study (30–59 years at baseline) together with circulating levels of vitamins A, C, E, and carotenoids (α-carotene, β-carotene, lutein and lycopene). Markers of socioeconomic deprivation in adulthood were consistently as strongly associated with lower vitamin C and carotenoid levels as markers of adverse lifestyle; the inverse association with overcrowding was particularly consistent (vitamin C and carotenoids range from 19.1% [95% CI 30.3–6.0] to 38.8% [49.9–25.3] lower among those in overcrowded residencies). These associations were consistent after adjusting for month, classical CVD risk factors, body mass index, physical activity, vitamin supplements, dietary fat and fibre intake. Similar, but weaker, associations were seen for childhood markers of deprivation. The association of vitamin A or E were strikingly different; several adult adverse lifestyle factors associated with higher levels of vitamin A and E, including high alcohol intake for vitamin A (9.5% [5.7–13.5]) and waist hip ratio for vitamin E (9.5% [4.8–14.4]), with the latter associations partially explained by classical risk factors, particularly cholesterol levels.&lt;/p&gt; &lt;p&gt;&lt;b&gt;Conclusions:&lt;/b&gt; Plasma vitamin C and carotenoids have strong inverse associations with adulthood markers of social deprivation, whereas vitamin A and E appear positively related to specific adverse lifestyle factors. These findings should help researchers better contextualize blood antioxidant vitamin levels by illustrating the potential limitations associated with making causal inferences without consideration of social deprivation.&lt;/p&gt

Detection of myxoma viruses encoding a defective M135R gene from clinical cases of myxomatosis; possible implications for the role of the M135R protein as a virulence factor

<p>Abstract</p> <p>Background</p> <p>Myxoma virus is a member of the <it>Poxviridae </it>and causes disease in European rabbits. Laboratory confirmation of the clinical disease, which occurs in the autumn of most years in Denmark, has been achieved previously using antigen ELISA and electron microscopy.</p> <p>Results</p> <p>An unusually large number of clinically suspected cases of myxomatosis were observed in Denmark during 2007. Myxoma virus DNA was detected, using a new real time PCR assay which targets the M029L gene, in over 70% of the clinical samples submitted for laboratory confirmation. Unexpectedly, further analysis revealed that a high proportion of these viral DNA preparations contained a frame-shift mutation within the M135R gene that has previously been identified as a virulence factor. This frame-shift mutation results in expression of a greatly truncated product. The same frame-shift mutation has also been found recently within an avirulent strain of myxoma virus (6918). However, three other frame-shift mutations found in this strain (in the genes M009L, M036L and M148R) were not shared with the Danish viruses but a single nucleotide deletion in the M138R/M139R intergenic region was a common feature.</p> <p>Conclusions</p> <p>It appears that expression of the full-length myxoma virus M135R protein is not required for virulence in rabbits. Hence, the frame-shift mutation in the M135R gene in the nonpathogenic 6918 virus strain is not sufficient to explain the attenuation of this myxoma virus but one/some of the other frame-shift mutations alone or in conjunction with one/some of the thirty two amino acid substitutions must also contribute. The real time PCR assay for myxoma virus is a useful diagnostic tool for laboratory confirmation of suspected cases of myxomatosis.</p

Identification and Specification of the Mouse Skeletal Stem Cell

SummaryHow are skeletal tissues derived from skeletal stem cells? Here, we map bone, cartilage, and stromal development from a population of highly pure, postnatal skeletal stem cells (mouse skeletal stem cells, mSSCs) to their downstream progenitors of bone, cartilage, and stromal tissue. We then investigated the transcriptome of the stem/progenitor cells for unique gene-expression patterns that would indicate potential regulators of mSSC lineage commitment. We demonstrate that mSSC niche factors can be potent inducers of osteogenesis, and several specific combinations of recombinant mSSC niche factors can activate mSSC genetic programs in situ, even in nonskeletal tissues, resulting in de novo formation of cartilage or bone and bone marrow stroma. Inducing mSSC formation with soluble factors and subsequently regulating the mSSC niche to specify its differentiation toward bone, cartilage, or stromal cells could represent a paradigm shift in the therapeutic regeneration of skeletal tissues