Reverse Engineering Time Discrete Finite Dynamical Systems: A Feasible Undertaking?

With the advent of high-throughput profiling methods, interest in reverse engineering the structure and dynamics of biochemical networks is high. Recently an algorithm for reverse engineering of biochemical networks was developed by Laubenbacher and Stigler. It is a top-down approach using time discrete dynamical systems. One of its key steps includes the choice of a term order, a technicality imposed by the use of Gröbner-bases calculations. The aim of this paper is to identify minimal requirements on data sets to be used with this algorithm and to characterize optimal data sets. We found minimal requirements on a data set based on how many terms the functions to be reverse engineered display. Furthermore, we identified optimal data sets, which we characterized using a geometric property called “general position”. Moreover, we developed a constructive method to generate optimal data sets, provided a codimensional condition is fulfilled. In addition, we present a generalization of their algorithm that does not depend on the choice of a term order. For this method we derived a formula for the probability of finding the correct model, provided the data set used is optimal. We analyzed the asymptotic behavior of the probability formula for a growing number of variables n (i.e. interacting chemicals). Unfortunately, this formula converges to zero as fast as , where and . Therefore, even if an optimal data set is used and the restrictions in using term orders are overcome, the reverse engineering problem remains unfeasible, unless prodigious amounts of data are available. Such large data sets are experimentally impossible to generate with today's technologies

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θsil=~0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima's D=0.316 vs D=−0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (πtot=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4Ne generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

<p>Abstract</p> <p>Background</p> <p>High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera.</p> <p>Results</p> <p>We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of <it>Eucalyptus </it>from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for <it>E. grandis</it>. A systematic assessment of <it>in silico </it>SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous <it>in silico </it>constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species.</p> <p>SNP reliability was high across nine <it>Eucalyptus </it>species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased.</p> <p>Conclusions</p> <p>This study indicates that the GGGT performs well both within and across species of <it>Eucalyptus </it>notwithstanding its nucleotide diversity ≥2%. The development of a much larger array of informative SNPs across multiple <it>Eucalyptus </it>species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in <it>Eucalyptus</it>.</p

In vitro phosphorylation as tool for modification of silk and keratin fibrous materials

An overview is given of the recent work on in vitro enzymatic phosphorylation of silk fibroin and human hair keratin. Opposing to many chemical "conventional" approaches, enzymatic phosphorylation is in fact a mild reaction and the treatment falls within "green chemistry" approach. Silk and keratin are not phosphorylated in vivo, but in vitro. This enzyme-driven modification is a major technological breakthrough. Harsh chemical chemicals are avoided, and mild conditions make enzymatic phosphorylation a real "green chemistry" approach. The current communication presents a novel approach stating that enzyme phosphorylation may be used as a tool to modify the surface charge of biocompatible materials such as keratin and silk

Second site escape of a T20-dependent HIV-1 variant by a single amino acid change in the CD4 binding region of the envelope glycoprotein

BACKGROUND: We previously described the selection of a T20-dependent human immunodeficiency virus type-1 (HIV-1) variant in a patient on T20 therapy. The fusion inhibitor T20 targets the viral envelope (Env) protein by blocking a conformational switch that is critical for viral entry into the host cell. T20-dependent viral entry is the result of 2 mutations in Env (GIA-SKY), creating a protein that undergoes a premature conformational switch, and the presence of T20 prevents this premature switch and rescues viral entry. In the present study, we performed 6 independent evolution experiments with the T20-dependent HIV-1 variant in the absence of T20, with the aim to identify second site compensatory changes, which may provide new mechanistic insights into Env function and the T20-dependence mechanism. RESULTS: Escape variants with improved replication capacity appeared within 42 days in 5 evolution cultures. Strikingly, 3 cultures revealed the same single amino acid change in the CD4 binding region of Env (glycine at position 431 substituted for arginine: G431R). This mutation was sufficient to abolish the T20-dependence phenotype and restore viral replication in the absence of T20. The GIA-SKY-G431R escape variant produces an Env protein that exhibits reduced syncytia formation and reduced cell-cell fusion activity. The escape variant was more sensitive to an antibody acting on an early gp41 intermediate, suggesting that the G431R mutation helps preserve a pre-fusion Env conformation, similar to T20 action. The escape variant was also less sensitive to soluble CD4, suggesting a reduced CD4 receptor affinity. CONCLUSION: The forced evolution experiments indicate that the premature conformational switch of the T20-dependent HIV-1 Env variant (GIA-SKY) can be corrected by a second site mutation in Env (GIA-SKY-G431R) that affects the interaction with the CD4 receptor

Misregulation of Scm3p/HJURP Causes Chromosome Instability in Saccharomyces cerevisiae and Human Cells

The kinetochore (centromeric DNA and associated proteins) is a key determinant for high fidelity chromosome transmission. Evolutionarily conserved Scm3p is an essential component of centromeric chromatin and is required for assembly and function of kinetochores in humans, fission yeast, and budding yeast. Overexpression of HJURP, the mammalian homolog of budding yeast Scm3p, has been observed in lung and breast cancers and is associated with poor prognosis; however, the physiological relevance of these observations is not well understood. We overexpressed SCM3 and HJURP in Saccharomyces cerevisiae and HJURP in human cells and defined domains within Scm3p that mediate its chromosome loss phenotype. Our results showed that the overexpression of SCM3 (GALSCM3) or HJURP (GALHJURP) caused chromosome loss in a wild-type yeast strain, and overexpression of HJURP led to mitotic defects in human cells. GALSCM3 resulted in reduced viability in kinetochore mutants, premature separation of sister chromatids, and reduction in Cse4p and histone H4 at centromeres. Overexpression of CSE4 or histone H4 suppressed chromosome loss and restored levels of Cse4p at centromeres in GALSCM3 strains. Using mutant alleles of scm3, we identified a domain in the N-terminus of Scm3p that mediates its interaction with CEN DNA and determined that the chromosome loss phenotype of GALSCM3 is due to centromeric association of Scm3p devoid of Cse4p/H4. Furthermore, we determined that similar to other systems the centromeric association of Scm3p is cell cycle regulated. Our results show that altered stoichiometry of Scm3p/HJURP, Cse4p, and histone H4 lead to defects in chromosome segregation. We conclude that stringent regulation of HJURP and SCM3 expression are critical for genome stability

Impact of Human Management on the Genetic Variation of Wild Pepper, Capsicum annuum var. glabriusculum

Management of wild peppers in Mexico has occurred for a long time without clear phenotypic signs of domestication. However, pre-domestication management could have implications for the population's genetic richness. To test this hypothesis we analysed 27 wild (W), let standing (LS) and cultivated (C) populations, plus 7 samples from local markets (LM), with nine polymorphic microsatellite markers. Two hundred and fifty two alleles were identified, averaging 28 per locus. Allele number was higher in W, and 15 and 40% less in LS and C populations, respectively. Genetic variation had a significant population structure. In W populations, structure was associated with ecological and geographic areas according to isolation by distance. When LM and C populations where included in the analysis, differentiation was no longer apparent. Most LM were related to distant populations from Sierra Madre Oriental, which represents their probable origin. Historical demography shows a recent decline in all W populations. Thus, pre-domestication human management is associated with a significant reduction of genetic diversity and with a loss of differentiation suggesting movement among regions by man. Measures to conserve wild and managed populations should be implemented to maintain the source and the architecture of genetic variation in this important crop relative

Protein-Binding Microarray Analysis of Tumor Suppressor AP2α Target Gene Specificity

Cheap and massively parallel methods to assess the DNA-binding specificity of transcription factors are actively sought, given their prominent regulatory role in cellular processes and diseases. Here we evaluated the use of protein-binding microarrays (PBM) to probe the association of the tumor suppressor AP2α with 6000 human genomic DNA regulatory sequences. We show that the PBM provides accurate relative binding affinities when compared to quantitative surface plasmon resonance assays. A PBM-based study of human healthy and breast tumor tissue extracts allowed the identification of previously unknown AP2α target genes and it revealed genes whose direct or indirect interactions with AP2α are affected in the diseased tissues. AP2α binding and regulation was confirmed experimentally in human carcinoma cells for novel target genes involved in tumor progression and resistance to chemotherapeutics, providing a molecular interpretation of AP2α role in cancer chemoresistance. Overall, we conclude that this approach provides quantitative and accurate assays of the specificity and activity of tumor suppressor and oncogenic proteins in clinical samples, interfacing genomic and proteomic assays