Does replication groups scoring reduce false positive rate in SNP interaction discovery?

BACKGROUNG. Computational methods that infer single nucleotide polymorphism (SNP) interactions from phenotype data may uncover new biological mechanisms in non-Mendelian diseases. However, practical aspects of such analysis face many problems. Present experimental studies typically use SNP arrays with hundreds of thousands of SNPs but record only hundreds of samples. Candidate SNP pairs inferred by interaction analysis may include a high proportion of false positives. Recently, Gayan et al. (2008) proposed to reduce the number of false positives by combining results of interaction analysis performed on subsets of data (replication groups), rather than analyzing the entire data set directly. If performing as hypothesized, replication groups scoring could improve interaction analysis and also any type of feature ranking and selection procedure in systems biology. Because Gayan et al. do not compare their approach to the standard interaction analysis techniques, we here investigate if replication groups indeed reduce the number of reported false positive interactions. RESULTS. A set of simulated and false interaction-imputed experimental SNP data sets were used to compare the inference of SNP-SNP interactions by means of replication groups to the standard approach where the entire data set was directly used to score all candidate SNP pairs. In all our experiments, the inference of interactions from the entire data set (e.g. without using the replication groups) reported fewer false positives. CONCLUSIONS. With respect to the direct scoring approach the utility of replication groups does not reduce false positive rates, and may, depending on the data set, often perform worse

Nanoparticle ordering in sandwiched polymer brushes

The organization of nano-particles inside grafted polymer layers is governed by the interplay of polymer-induced entropic interactions and the action of externally applied fields. Earlier work had shown that strong external forces can drive the formation of colloidal structures in polymer brushes. Here we show that external fields are not essential to obtain such colloidal patterns: we report Monte Carlo and Molecular dynamics simulations that demonstrate that ordered structures can be achieved by compressing a `sandwich' of two grafted polymer layers, or by squeezing a coated nanotube, with nano-particles in between. We show that the pattern formation can be efficiently controlled by the applied pressure, while the characteristic length--scale, i.e. the typical width of the patterns, is sensitive to the length of the polymers. Based on the results of the simulations, we derive an approximate equation of state for nano-sandwiches.Comment: 18 pages, 4 figure

New insights on limes and lemons origin from nuclear and cytoplasmic markers genenotyping and targeted nuclear gene sequencing

It is believed that Citrus medica, C. maxima, C. reticulate and C. micrantha have generated all cultivated citrus species. Depending on the classification, lemons and limes are classified either into two species, C. limon and C. aurantifolia (Swingle and Reece, 1967) or into more than 30 (Tanaka, 1977). In order to study the molecular phylogeny of this citrus group, we analyzed 23 targeted sequenced nuclear genes and used three mitochondrial and five chloroplastic markers for 15 lemons and limes compared with representatives of the four basic taxa. We observed three main groups, each one derived from direct interspecific hybridizations: (1) the Mexican lime group (C. aurantifolia), including C. macrophylla, arising from hybridizations between papeda (C. micrantha) and citron (C. medica); (2) the yellow lemon group (C. limon) that are hybrids between sour orange (C. aurantium, which is believed to be a hybrid between C. maxima and C. reticulata) and citron; and (3) a rootstock lemon/lime group (Rough lemon and Rangpur lime) that are hybrids between an acid mandarin and a citron. We also identified different probable backcrosses and genotypes with more complex origin. None of the analyzed limes and lemons shared the C. medica cytoplasm, while this taxon is the common nuclear contributor of all limes and lemons. Limes and lemons appear to be a very complex citrus varietal group with the contribution of the 4 basic taxa. Neither the Swingle and Reece classification nor the Tanaka fit with the genetic evidence. (Résumé d'auteur

Rational design of molecularly imprinted polymers.

Molecular imprinting is the process whereby a polymer matrix is cross-linked in the presence of molecules with surface sites that can bind selectively to certain ligands on the polymer. The cross-linking process endows the polymer matrix with a chemical 'memory', such that the target molecules can subsequently be recognized by the matrix. We present a simple model that accounts for the key features of this molecular recognition. Using a combination of analytical calculations and Monte Carlo simulations, we show that the model can account for the binding of rigid particles to an imprinted polymer matrix with valence-limited interactions. We show how the binding multivalency and the polymer material properties affect the efficiency and selectivity of molecular imprinting. Our calculations allow us to formulate design criteria for optimal molecular imprinting.This work was supported by the Fundamental Research Funds for the Central Universities of P. R. China, ERC Advanced Grant 227758 (COLSTRUCTION), ITN grant 234810 (COMPPLOIDS) and by EPSRC Programme Grant EP/I001352/1. TC acknowledges support from the Herchel Smith fund.This is the final version of the article. It first appeared from RSC via http://dx.doi.org/10.1039/C5SM02144

Layering, freezing and re-entrant melting of hard spheres in soft confinement

Confinement can have a dramatic effect on the behavior of all sorts of particulate systems and it therefore is an important phenomenon in many different areas of physics and technology. Here, we investigate the role played by the softness of the confining potential. Using grand canonical Monte Carlo simulations, we determine the phase diagram of three-dimensional hard spheres that in one dimension are constrained to a plane by a harmonic potential. The phase behavior depends strongly on the density and on the stiffness of the harmonic confinement. Whilst we find the familiar sequence of confined hexagonal and square-symmetric packings, we do not observe any of the usual intervening ordered phases. Instead, the system phase separates under strong confinement, or forms a layered re-entrant liquid phase under weaker confinement. It is plausible that this behavior is due to the larger positional freedom in a soft confining potential and to the contribution that the confinement energy makes to the total free energy. The fact that specific structures can be induced or suppressed by simply changing the confinement conditions (e.g. in a dielectrophoretic trap) is important for applications that involve self-assembled structures of colloidal particles.Comment: 5 pages, 5 figure

The Effect of Attractive Interactions and Macromolecular Crowding on Crystallins Association.

In living systems proteins are typically found in crowded environments where their effective interactions strongly depend on the surrounding medium. Yet, their association and dissociation needs to be robustly controlled in order to enable biological function. Uncontrolled protein aggregation often causes disease. For instance, cataract is caused by the clustering of lens proteins, i.e., crystallins, resulting in enhanced light scattering and impaired vision or blindness. To investigate the molecular origins of cataract formation and to design efficient treatments, a better understanding of crystallin association in macromolecular crowded environment is needed. Here we present a theoretical study of simple coarse grained colloidal models to characterize the general features of how the association equilibrium of proteins depends on the magnitude of intermolecular attraction. By comparing the analytic results to the available experimental data on the osmotic pressure in crystallin solutions, we identify the effective parameters regimes applicable to crystallins. Moreover, the combination of two models allows us to predict that the number of binding sites on crystallin is small, i.e. one to three per protein, which is different from previous estimates. We further observe that the crowding factor is sensitive to the size asymmetry between the reactants and crowding agents, the shape of the protein clusters, and to small variations of intermolecular attraction. Our work may provide general guidelines on how to steer the protein interactions in order to control their association

Heterogeneous computing architecture for fast detection of SNP-SNP interactions

The extent of data in a typical genome-wide association study (GWAS) poses considerable computational challenges to software tools for gene-gene interaction discovery. Exhaustive evaluation of all interactions among hundreds of thousands to millions of single nucleotide polymorphisms (SNPs) may require weeks or even months of computation. Massively parallel hardware within a modern Graphic Processing Unit (GPU) and Many Integrated Core (MIC) coprocessors can shorten the run time considerably. While the utility of GPU-based implementations in bioinformatics has been well studied, MIC architecture has been introduced only recently and may provide a number of comparative advantages that have yet to be explored and tested. We have developed a heterogeneous, GPU and Intel MIC-accelerated software module for SNP-SNP interaction discovery to replace the previously single-threaded computational core in the interactive web-based data exploration program SNPsyn. We report on differences between these two modern massively parallel architectures and their software environments. Their utility resulted in an order of magnitude shorter execution times when compared to the single-threaded CPU implementation. GPU implementation on a single Nvidia Tesla K20 runs twice as fast as that for the MIC architecture-based Xeon Phi P5110 coprocessor, but also requires considerably more programming effort. General purpose GPUs are a mature platform with large amounts of computing power capable of tackling inherently parallel problems, but can prove demanding for the programmer. On the other hand the new MIC architecture, albeit lacking in performance reduces the programming effort and makes it up with a more general architecture suitable for a wider range of problems

Multilocus haplotyping by parallel sequencing to decipher the interspecific mosaic genome structure of cultivated citrus

The most important economic Citrus species originated from natural interspecific hybridization between four ancestral taxa (C. reticulata, C. maxima, C. medica and C. micrantha) with limited further interspecific recombination due to apomixis and vegetative propagation. Such reticulate evolution coupled with vegetative propagation results in genomes that are mosaics of large chromosome fragments of the basic taxa, in frequent interspecific heterozygosity. Breeding of these species is hampered by their complex heterozygous genomic structures. Haplotyping of multiple gene fragments along the genome should be a powerful approach to resolve the evolutionary history of the gene pools, to reveal the admixture genomic structure of current species and to develop innovative breeding schemes. We have analysed the efficiency of parallel sequencing with 454 methodology to decipher the hybrid structure of modern citrus species and cultivars along chromosome 2. Four hundred fifty four amplicon libraries were established with the fluidigm array system for 48 genotypes and 16 gene fragments of chromosome 2. Haplotypes were established from the reads of each accession and phylogenetic analyses were performed from the haplotypic data of each gene fragment. The length of 454 reads and the level of differentiation between the ancestral taxa of modern citrus allowed efficient haplotype phylogenetic assignations for 12 of the 16 gene fragments. The analysis of the mixed genomic structure of modern species and cultivars (i) revealed C. maxima introgressions in modern mandarins; (ii) was consistent with previous hypothesis regarding the origin of secondary species; and (iii) provided a new picture of the evolution of chromosome 2. Perspectives to rebuild the main secondary species from the basic taxa are discussed. (Résumé d'auteur