Inferring the effects of demography and selection on Drosophila melanogaster populations from a chromosome-wide scan of DNA variation

Identifying regions of the Drosophila melanogaster genome that have been recent targets of positive Darwinian selection will provide evidence for adaptations that have helped this species to colonize temperate habitats. We have begun a search for such genomic regions by analyzing multiple loci (about 250) dispersed across the X chromosome in a putatively ancestral population from East Africa and a derived European population. For both populations we found evidence for past changes in population size. We estimated that a major bottleneck associated with the colonization of Europe occurred about 3,500-16,000 years ago. We also found that while this bottleneck can account for most of the reduction in variation observed in the European sample, there is a deficit of polymorphism in some genomic regions that cannot be explained by demography alone.     

ParIS Genome Rearrangement server

SUMMARY: ParIS Genome Rearrangement is a web server for a Bayesian analysis of unichromosomal genome pairs. The underlying model allows inversions, transpositions and inverted transpositions. The server generates a Markov chain using a Partial Importance Sampler technique, and samples trajectories of mutations from this chain. The user can specify several marginalizations to the posterior: the posterior distribution of number of mutations needed to transform one genome into another, length distribution of mutations, number of mutations that have occurred at a given site. Both text and graphical outputs are available. We provide a limited server, a downloadable unlimited server that can be installed locally on any linux/Unix operating system, and a database of mitochondrial gene orders.     

Monitoring presence and advance of the alien invasive western corn rootworm beetle in eastern Slovenia with highly sensitive Metcalf traps

The American Chrysomelid beetle Diabrotica virgifera virgifera LeConte (D.v.v.), also called the western corn rootworm, spread from the location of its original introduction into Europe, Belgrade airport (BACA 1993), in all directions. Within a decade it occupied almost all countries of South-eastern and Central Europe. However, it reached Slovenia as late as 2003: Only 19 specimen were found in maize fields of the eastern and also western provinces. Already in the summer of 2004, their number had risen to 386 which were mainly found in the eastern provinces near the borders to Croatia, Austria and Hungary. For their monitoring, a simple trap is being described which can be acquired in high numbers at a very low price and can guaranty a most sensitive detection of beetles. The Metcalf cup trap which in Slovenia so far has been unknown performed the task of monitoring quite well at five locations. It will also facilitate the future search for new and increasing infestations. Comparing 2004 with the year 2003, D.v.v. expanded its range in eastern Slovenia by about 15 km. These new infestations will include territories in which fields of the regionally Important oil seed pumpkin Cucurbita pepo are located. With important traffic connections between South-eastern and Central Europe, Slovenia will occupy a bridgehead function in the preventive protection of maize from D.v.v. spreading into more northerly European regions including western Austria and southern Germany.     

A promising biotechnical approach to pest management of Diabrotica virgifera virgifera in Illinois maize fields under kairomonal shielding with the new MSD technique

Environmentally compatible and sustainable plant protection requires novel approaches to pest management characterized by minimal emphasis on toxicants. Classical toxicants traditionally dominated economic entomology for half a century. But worldwide problems with environmental pollution and with increasing resistance levels in all major pesticide classes and in many key insect species including Diabrotica virgifera virgifera (D.v.v.) strongly advocate a rethinking and a change in management paradigms used. Soft, minimally invasive, biological, biotechnical and cultural approaches should replace hard pesticides which are in favor up to now. Fortunately, pheromones, kairomones, plant attractants, better traps, new plant varieties and cultural methods like crop rotation, in short more sophisticated methods are now available as pressure for finding and exploring novel strategies increases. Facing this situation, a new biotechnical approach of population reduction of D.v.v., called "MSD" technique, is introduced. MSD is characterized as an approach combining mass trapping, shielding and deflecting of adult insects along an invisible odor barrier of synthetic kairomone which diminishes the flux of insects across a high capacity trap line baited with kairomone, thus reducing both the population fluctuation and number and its reproductive success within the shielded area. In the case of D.v.v. in Zea mays fields, effects realized by the MSD technique have been measured simultaneously by a number of independent criteria during the summers of 2003 and 2004 at 2 different locations in Illinois maize fields of up to one half hectare size. Results observed are statistically significant and cannot be explained by mass trapping alone. There is also an additional shielding and deflection, in short "diversion" effect whose basic sensory and behavioral mechanisms call for future exploration.     

Ascent exhalations of Antarctic fur seals: a behavioural adaptation for breath-hold diving?

Novel observations collected from video, acoustic and conductivity sensors showed that Antarctic fur seals consistently exhale during the last 50-85% of ascent from all dives (10-160 m, n > 8000 dives from 50 seals). The depth of initial bubble emission was best predicted by maximum dive depth, suggesting an underlying physical mechanism. Bubble sound intensity recorded from one seal followed predictions of a simple model based on venting expanding lung air with decreasing pressure. Comparison of air release between dives, together with lack of variation in intensity of thrusting movement during initial descent regardless of ultimate dive depth, suggested that inhaled diving lung volume was constant for all dives. The thrusting intensity in the final phase of ascent was greater for dives in which ascent exhalation began at a greater depth, suggesting an energetic cost to this behaviour, probably as a result of loss of buoyancy from reduced lung volume. These results suggest that fur seals descend with full lung air stores, and thus face the physiological consequences of pressure at depth. We suggest that these regular and predictable ascent exhalations could function to reduce the potential for a precipitous drop in blood oxygen that would result in shallow-water blackout.     

A secreted form of P-cadherin is expressed in malignant melanoma

Cadherins are Ca-dependent homophilic cell-cell adhesion molecules which are responsible for correct location of cells and tissue integrity. They are critical for development and maintenance of epithelial architecture. Aberrantly expressed cadherins are known to be involved in malignant transformation of different types of tissues. In this study, we show the expression of a short truncated 50 kDa form of the N-terminal part of P-cadherin in seven melanoma cell lines compared to melanocytes and keratinocytes. In vitro protein analysis on cell culture supernatant as well as immunohistochemistry of primary and metastatic melanoma tissue revealed the expression of this short form of P-cadherin. Furthermore, analysis showed that this short 50 kDa form of P-cadherin is secreted by melanoma cells in contrast to the membrane bound form in melanocytes. Analysis on mRNA level detected only exon 1 to 10 of P-cadherin resulting in the 50 kDa form missing the transmembrane and cytoplasmatic region. Genomic sequence analysis did not show any mutations in melanoma cells neither in the exons nor in the exon-intron boundaries. Furthermore, there was no loss of exons 11-16 on the genomic level. Functionally, the secreted form of P-cadherin could play a role as regulator of the homophilic interaction between P-cadherin molecules by antagonizing their biological role acting as a dominant negative form to interrupt cell-cell attachment.     

BGX: a fully Bayesian integrated approach to the analysis of Affymetrix GeneChip data

We present Bayesian hierarchical models for the analysis of Affymetrix GeneChip data. The approach we take differs from other available approaches in two fundamental aspects. Firstly, we aim to integrate all processing steps of the raw data in a common statistically coherent framework, allowing all components and thus associated errors to be considered simultaneously. Secondly, inference is based on the full posterior distribution of gene expression indices and derived quantities, such as fold changes or ranks, rather than on single point estimates. Measures of uncertainty on these quantities are thus available. The models presented represent the first building block for integrated Bayesian Analysis of Affymetrix GeneChip data: the models take into account additive as well as multiplicative error, gene expression levels are estimated using perfect match and a fraction of mismatch probes and are modeled on the log scale. Background correction is incorporated by modeling true signal and cross-hybridization explicitly, and a need for further normalization is considerably reduced by allowing for array-specific distributions of nonspecific hybridization. When replicate arrays are available for a condition, posterior distributions of condition-specific gene expression indices are estimated directly, by a simultaneous consideration of replicate probe sets, avoiding averaging over estimates obtained from individual replicate arrays. The performance of the Bayesian model is compared to that of standard available point estimate methods on subsets of the well known GeneLogic and Affymetrix spike-in data. The Bayesian model is found to perform well and the integrated procedure presented appears to hold considerable promise for further development.     

Modeling the effect of an associated single-nucleotide polymorphism in linkage studies

For linkage analysis in affected sibling pairs, we propose a regression model to incorporate information from a disease-associated single-nucleotide polymorphism located under the linkage peak. This model can be used to study if the associated single-nucleotide polymorphism marker partly explains the original linkage peak. Two sources of information are used for performing this task, namely the genotypes of the parents and the genotypes of the siblings. We applied the methods to three significantly disease-associated single-nucleotide polymorphisms and five microsatellite markers at the end of chromosome 3 of replicate 1 of Aipotu population. Two out of five of the microsatellite markers showed a LOD score higher than 3. The question to be answered was whether one of the single-nucleotide polymorphisms partly explains these high LOD scores. We did not have the answers when we analyzed the data.     

Bayesian coestimation of phylogeny and sequence alignment

Background  

Two central problems in computational biology are the determination of the alignment and phylogeny of a set of biological sequences. The traditional approach to this problem is to first build a multiple alignment of these sequences, followed by a phylogenetic reconstruction step based on this multiple alignment. However, alignment and phylogenetic inference are fundamentally interdependent, and ignoring this fact leads to biased and overconfident estimations. Whether the main interest be in sequence alignment or phylogeny, a major goal of computational biology is the co-estimation of both.