Microsatellite primers for Ficus racemosa and Ficus rubiginosa

We developed 11 polymorphic microsatellite loci each for the figs Ficus (Sycomorus) racemosa and Ficus (Urostigma) rubiginosa from AG- and TG-enriched genomic libraries. These 22 loci were investigated for cross-species amplification and polymorphism in 17–21 F. racemosa and 16–24 F. rubiginosa individuals from Townsville, Australia. Observed heterozygosities range from 0.12 to 0.90 in F. racemosa and from 0.25 to 1.0 in F. rubiginosa.

100 years on: a century of genetics

In 2006 we celebrated the centenary of a remarkable year that saw the birth of genetics as a scientific discipline. This birth had its origins in horticulture and was supervised by a remarkable Cambridge academic, accompanied by a loyal group of female colleagues who worked together in underfunded conditions with little institutional support. Despite this deprivation, they established the foundations of an ongoing revolution, with huge academic and commercial consequences that we can recognize today in the shape of genomics and its application to biomedicine.

Benchmarking consensus model quality assessment for protein fold recognition

Background: Selecting the highest quality 3D model of a protein structure from a number of alternatives remains an important challenge in the field of structural bioinformatics. Many Model Quality Assessment Programs (MQAPs) have been developed which adopt various strategies in order to tackle this problem, ranging from the so called "true" MQAPs capable of producing a single energy score based on a single model, to methods which rely on structural comparisons of multiple models or additional information from meta-servers. However, it is clear that no current method can separate the highest accuracy models from the lowest consistently. In this paper, a number of the top performing MQAP methods are benchmarked in the context of the potential value that they add to protein fold recognition. Two novel methods are also described: ModSSEA, which based on the alignment of predicted secondary structure elements and ModFOLD which combines several true MQAP methods using an artificial neural network. Results: The ModSSEA method is found to be an effective model quality assessment program for ranking multiple models from many servers, however further accuracy can be gained by using the consensus approach of ModFOLD. The ModFOLD method is shown to significantly outperform the true MQAPs tested and is competitive with methods which make use of clustering or additional information from multiple servers. Several of the true MQAPs are also shown to add value to most individual fold recognition servers by improving model selection, when applied as a post filter in order to re-rank models. Conclusion: MQAPs should be benchmarked appropriately for the practical context in which they are intended to be used. Clustering based methods are the top performing MQAPs where many models are available from many servers; however, they often do not add value to individual fold recognition servers when limited models are available. Conversely, the true MQAP methods tested can often be used as effective post filters for re-ranking few models from individual fold recognition servers and further improvements can be achieved using a consensus of these methods.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Transcriptional profiling of the LPS induced NF-kappaB response in macrophages

BACKGROUND: Exposure of macrophages to bacterial products such as lipopolysaccharide (LPS) results in activation of the NF-kappaB transcription factor, which orchestrates a gene expression programme that underpins the macrophage-dependent immune response. These changes include the induction or repression of a wide range of genes that regulate inflammation, cell proliferation, migration and cell survival. This process is tightly regulated and loss of control is associated with conditions such as septic shock, inflammatory diseases and cancer. To study this response, it is important to have in vitro model systems that reflect the behaviour of cells in vivo. In addition, it is necessary to understand the natural differences that can occur between individuals. In this report, we have investigated and compared the LPS response in macrophage derived cell lines and peripheral blood mononuclear cell (PBMC) derived macrophages. RESULTS: Gene expression profiles were determined following LPS treatment of THP-1 cells for 1 and 4 hours. LPS significantly induced or repressed 72 out of 465 genes selected as being known or putative NF-kappaB target genes, which exhibited 4 temporal patterns of expression. Results for 34 of these genes, including several genes not previously identified as LPS target genes, were validated using real time PCR. A high correlation between microarray and real time PCR data was found. Significantly, the LPS induced expression profile of THP-1 cells, as determined using real time PCR, was found to be very similar to that of human PBMC derived macrophages. Interestingly, some differences were observed in the LPS response between the two donor PBMC macrophage populations. Surprisingly, we found that the LPS response in U937 cells was dramatically different to both THP-1 and PBMC derived macrophages. CONCLUSION: This study revealed a dynamic and diverse transcriptional response to LPS in macrophages, involving both the induction and repression of gene expression in a time dependent manner. Moreover, we demonstrated that the LPS induced transcriptional response in the THP-1 cell line is very similar to primary PBMC derived macrophages. Therefore, THP-1 cells represent a good model system for studying the mechanisms of LPS and NF-kappaB dependent gene expression.

Burkholderia Hep_Hag autotransporter (BuHA) proteins elicit a strong antibody response during experimental glanders but not human melioidosis

BACKGROUND: The bacterial biothreat agents Burkholderia mallei and Burkholderia pseudomallei are the cause of glanders and melioidosis, respectively. Genomic and epidemiological studies have shown that B. mallei is a recently emerged, host restricted clone of B. pseudomallei. RESULTS: Using bacteriophage-mediated immunoscreening we identified genes expressed in vivo during experimental equine glanders infection. A family of immunodominant antigens were identified that share protein domain architectures with hemagglutinins and invasins. These have been designated Burkholderia Hep_Hag autotransporter (BuHA) proteins. A total of 110/207 positive clones (53%) of a B. mallei expression library screened with sera from two infected horses belonged to this family. This contrasted with 6/189 positive clones (3%) of a B. pseudomallei expression library screened with serum from 21 patients with culture-proven melioidosis. CONCLUSION: Members of the BuHA proteins are found in other Gram-negative bacteria and have been shown to have important roles related to virulence. Compared with other bacterial species, the genomes of both B. mallei and B. pseudomallei contain a relative abundance of this family of proteins. The domain structures of these proteins suggest that they function as multimeric surface proteins that modulate interactions of the cell with the host and environment. Their effect on the cellular immune response to B. mallei and their potential as diagnostics for glanders requires further study.

Characterization and PCR multiplexing of polymorphic microsatellite loci in cashew (Anacardium occidentale L.) and their cross-species utilization

Cashew (Anacardium occidentale L.) is the most economically important tropical nut crop in the world, and yet there are no sequence tagged site (STS) markers available for its study. Here we use an automated, high-throughput system to isolate cashew microsatellites from a non-enriched genomic library blotted onto membranes at high density for screening. Sixty-five sequences contained a microsatellite array, of which 21 proved polymorphic among a closely related seed garden population of 49 genotypes. Twelve markers were suitable for multiplex analysis. Of these, 10 amplified in all three related tropical tree species tested: Anacardium microcarpum, Anacardium pumilum and Anacardium nanum.

Expression of target and reference genes in Daphnia magna exposed to ibuprofen

Background: Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e. g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application. Results: We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l(-1). Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was upregulated following 24 h exposure to 63-81 mg IB l(-1). Conclusions: As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.

High throughput profile-profile based fold recognition for the entire human proteome

BACKGROUND: In order to maintain the most comprehensive structural annotation databases we must carry out regular updates for each proteome using the latest profile-profile fold recognition methods. The ability to carry out these updates on demand is necessary to keep pace with the regular updates of sequence and structure databases. Providing the highest quality structural models requires the most intensive profile-profile fold recognition methods running with the very latest available sequence databases and fold libraries. However, running these methods on such a regular basis for every sequenced proteome requires large amounts of processing power.In this paper we describe and benchmark the JYDE (Job Yield Distribution Environment) system, which is a meta-scheduler designed to work above cluster schedulers, such as Sun Grid Engine (SGE) or Condor. We demonstrate the ability of JYDE to distribute the load of genomic-scale fold recognition across multiple independent Grid domains. We use the most recent profile-profile version of our mGenTHREADER software in order to annotate the latest version of the Human proteome against the latest sequence and structure databases in as short a time as possible. RESULTS: We show that our JYDE system is able to scale to large numbers of intensive fold recognition jobs running across several independent computer clusters. Using our JYDE system we have been able to annotate 99.9% of the protein sequences within the Human proteome in less than 24 hours, by harnessing over 500 CPUs from 3 independent Grid domains. CONCLUSION: This study clearly demonstrates the feasibility of carrying out on demand high quality structural annotations for the proteomes of major eukaryotic organisms. Specifically, we have shown that it is now possible to provide complete regular updates of profile-profile based fold recognition models for entire eukaryotic proteomes, through the use of Grid middleware such as JYDE.

A phyloclimatic study of Cyclamen

Background: The impact of global climate change on plant distribution, speciation and extinction is of current concern. Examining species climatic preferences via bioclimatic niche modelling is a key tool to study this impact. There is an established link between bioclimatic niche models and phylogenetic diversification. A next step is to examine future distribution predictions from a phylogenetic perspective. We present such a study using Cyclamen (Myrsinaceae), a group which demonstrates morphological and phenological adaptations to its seasonal Mediterranean-type climate. How will the predicted climate change affect future distribution of this popular genus of garden plants? Results: We demonstrate phylogenetic structure for some climatic characteristics, and show that most Cyclamen have distinct climatic niches, with the exception of several wide-ranging, geographically expansive, species. We reconstruct climate preferences for hypothetical ancestral Cyclamen. The ancestral Cyclamen lineage has a preference for the seasonal Mediterranean climate characteristic of dry summers and wet winters. Future bioclimatic niches, based on BIOCLIM and Maxent models, are examined with reference to a future climate scenario for the 2050s. Over the next 50 years we predict a northward shift in the area of climatic suitability, with many areas of current distribution becoming climatically unsuitable. The area of climatic suitability for every Cyclamen species is predicted to decrease. For many species, there may be no areas with a suitable climate regardless of dispersal ability, these species are considered to be at high risk of extinction. This risk is examined from a phylogenetic perspective. Conclusion: Examining bioclimatic niches from a phylogenetic perspective permits novel interpretations of these models. In particular, reconstruction of ancestral niches can provide testable hypothesis about the historical development of lineages. In the future we can expect a northwards shift in climatic suitability for the genus Cyclamen. If this proves to be the case then dispersal is the best chance of survival, which seems highly unlikely for ant-dispersed Cyclamen. Human-assisted establishment of Cyclamen species well outside their native ranges offers hope and could provide the only means of dispersal to potentially suitable future environments. Even without human intervention the phylogenetic perspective demonstrates that major lineages could survive climate change even if many species are lost.

Resistance is costly: trade-offs between immunity, fecundity and survival in the pea aphid

Parasitoids are among the most important natural enemies of insects in many environments. Acyrthosiphon pisum, the pea aphid, is a common pest of the leguminous crops in temperate regions. Pea aphids are frequently attacked by a range of endoparasitic wasps, including the common aphidiine, Aphidius ervi. Immunity to parasitoid attack is thought to involve secondary symbiotic bacteria, the presence of which is associated with the death of the parasitoid egg. It has been suggested that there is a fecundity cost of resistance, as individuals carrying the secondary symbionts associated with parasitoid resistance have fewer offspring. Supporting this hypothesis, we find a positive relationship between fecundity and susceptibility to parasitoid attack. There is also a negative relationship between fecundity and off-plant survival time (which positively correlates with resistance to parasitoid attack). Taken together, these results suggest that the aphids can either invest in defence (parasitoid resistance, increased off-plant survival time) or reproduction, and speculate that this may be mediated by changes in the aphids' endosymbiont fauna. Furthermore, there is a positive relationship between aphid size and resistance, suggesting that successful resistance to parasitoid attack may involve physical, as well as physiological, defences.

API-CALC 1.0: a computer program for calculating the average probability of identity allowing for substructure, inbreeding and the presence of close relatives

Individual identification via DNA profiling is important in molecular ecology, particularly in the case of noninvasive sampling. A key quantity in determining the number of loci required is the probability of identity (PIave), the probability of observing two copies of any profile in the population. Previously this has been calculated assuming no inbreeding or population structure. Here we introduce formulae that account for these factors, whilst also accounting for relatedness structure in the population. These formulae are implemented in API-CALC 1.0, which calculates PIave for either a specified value, or a range of values, for F-IS and F-ST.

The Bayesian revolution in genetics

Bayesian statistics allow scientists to easily incorporate prior knowledge into their data analysis. Nonetheless, the sheer amount of computational power that is required for Bayesian statistical analyses has previously limited their use in genetics. These computational constraints have now largely been overcome and the underlying advantages of Bayesian approaches are putting them at the forefront of genetic data analysis in an increasing number of areas.

Amino terminal interaction in the prion protein identified using fusion to green fluorescent protein

In contrast to the well-characterized carboxyl domain, the amino terminal half of the mature cellular prion protein has no defined structure. Here, following fusion of mouse prion protein fragments to green fluorescence protein as a reporter of protein stability, we report extreme variability in fluorescence level that is dependent on the prion fragment expressed. In particular, exposure of the extreme amino terminus in the context of a truncated prion protein molecule led to rapid degradation, whereas the loss of only six amino terminal residues rescued high level fluorescence. Study of the precise endpoints and residue identity associated with high fluorescence suggested a domain within the amino terminal half of the molecule defined by a long-range intramolecular interaction between 23KKRPKP28 and 143DWED146 and dependent upon the anti-parallel beta-sheet ending at residue 169 and normally associated with the structurally defined carboxyl terminal domain. This previously unreported interaction may be significant for understanding prion bioactivity and for structural studies aimed at the complete prion structure.