A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

Refinement of defection strategies stabilizes cooperation

Cooperation among unrelated individuals is an enduring evolutionary riddle and a number of possible solutions have been suggested. Most of these suggestions attempt to refine cooperative strategies, while little attention is given to the fact that novel defection strategies can also evolve in the population. Especially in the presence of punishment to the defectors and public knowledge of strategies employed by the players, a defecting strategy that avoids getting punished by selectively cooperating only with the punishers can get a selective benefit over non-conditional defectors. Furthermore, if punishment ensures cooperation from such discriminating defectors, defectors who punish other defectors can evolve as well. We show that such discriminating and punishing defectors can evolve in the population by natural selection in a Prisoner’s Dilemma game scenario, even if discrimination is a costly act. These refined defection strategies destabilize unconditional defectors. They themselves are, however, unstable in the population. Discriminating defectors give selective benefit to the punishers in the presence of non-punishers by cooperating with them and defecting with others. However, since these players also defect with other discriminators they suffer fitness loss in the pure population. Among the punishers, punishing cooperators always benefit in contrast to the punishing defectors, as the latter not only defect with other punishing defectors but also punish them and get punished. As a consequence of both these scenarios, punishing cooperators get stabilized in the population. We thus show ironically that refined defection strategies stabilize cooperation. Furthermore, cooperation stabilized by such defectors can work under a wide range of initial conditions and is robust to mistakes.

Left-Handed DNA in Synthetic and Topologically Constrained Form V-DNA

We have investigated structural transitions in Poly(dG-dC) and Poly(dG-Me5dC) in order to understand the exact role of cations in stabilizing left-handed helical structures in specific sequences andthe biological role, if any, of these structures. From a novel temperature dependent transition it has been shown that a minor fluctuation in Na+ concentration at ambient temperature can bring about Β to Ζ transition. Forthe first time, wehave observed a novel double transition in poly(dG-Me5dC) as the Na+ concentration is gradually increased. This suggests that a minor fluctuation in Na+ concentration in conjunction with methylation may transform small stretches of CG sequences from one conformational state to another. These stretches could probably serve as sites for regulation. Supercoiled formV DNA reconstituted from pBR322 and pßG plasmids have been studied as model systems, in order to understand the nature and role of left-handed helical conformation in natural sequences. A large portion of DNA in form V, obtained by reannealing the two complementary singlestranded circles is forced to adopt left-handed double helical structure due to topological constraints (Lk = 0). Binding studies with Z-DNA specific antibody and spectroscopic studies confirm the presence of left-handed Z-structure in the pßG and pßR322 form V DNA. Cobalt hexamine chloride, which induces Z-form in Poly(dG-dC) stabilizes the Z-conformation in form V DNA even in the non-alternating purine-pyrimidine sequences. A reverse effect is observed with ethidium bromide. Interestingly, both topoisomerase I and II (from wheat germ) act effectively on form V DNA to give rise to a species having an electrophoretic mobility on agarose gel similar to that of open circular (form II) DNA. Whether this molecule is formed as a result of the left-handed helical segments of form V DNA undergoing a transition to the right-handed B-form during the topoisomerase action remains to be solved.

Automated topometric graph generation from floor plan analysis

The world is rich with information such as signage and maps to assist humans to navigate. We present a method to extract topological spatial information from a generic bitmap floor plan and build a topometric graph that can be used by a mobile robot for tasks such as path planning and guided exploration. The algorithm first detects and extracts text in an image of the floor plan. Using the locations of the extracted text, flood fill is used to find the rooms and hallways. Doors are found by matching SURF features and these form the connections between rooms, which are the edges of the topological graph. Our system is able to automatically detect doors and differentiate between hallways and rooms, which is important for effective navigation. We show that our method can extract a topometric graph from a floor plan and is robust against ambiguous cases most commonly seen in floor plans including elevators and stairwells.

Exploring methods of preparing functional cartilage-bone xenografts for joint repair

This thesis explores the feasibility of donor-receiver concept for joint replacement where cartilage-bone tissues can be taken from either human or other mammals and prepared scientifically for repairing focal joint defects in knees, hips and shoulders. The manufactured construct is immunologically inert and is capable of acting as a scaffold for engineering new cartilage-bone laminates when placed in the joint. Innovative manufacturing procedures and assessment techniques were developed for appraising this tissue-based scaffold. This research has demonstrated that tissue replacement technology can be applied in situations where blood vessels are absent such as in articular cartilage.

Topology-Guided Tessellation Of Quadratic Elements

Topology-based methods have been successfully used for the analysis and visualization of piecewise-linear functions defined on triangle meshes. This paper describes a mechanism for extending these methods to piecewise-quadratic functions defined on triangulations of surfaces. Each triangular patch is tessellated into monotone regions, so that existing algorithms for computing topological representations of piecewise-linear functions may be applied directly to the piecewise-quadratic function. In particular, the tessellation is used for computing the Reeb graph, a topological data structure that provides a succinct representation of level sets of the function.

Conformational Dynamics in Alkyl Chains of an Anchored Bilayer: A Molecular Dynamics Study

Molecular dynamics (MD) simulations are reported for an anchored bilayer formed by the intercalation of cetyl trimethyl ammonium (CTA) and CH3(CH2)15N+(CH3) ions in a layered solid, CdPS3. The intercalated CTA ions are organized with the cationic headgroups tethered to the inorganic sheet and the hydrocarbon tails arranged as bilayers. Simulations were performed at three temperatures, 65, 180, and 298 K, using an isothermal−isobaric ensemble that was subsequently switched once macroscopic parameters had converged to a canonical isothermal−isochoric ensemble. The simulations are able to reproduce the experimental features of this system, including the formation of the bilayer and layer-to-layer separation distance. An analysis of the conformation of the chains showed that at all three temperatures a fraction of the alkyl chains retained a planar all-trans conformation, and that gauche bonds occurred as part of a “kink” (gauche+−trans−gauche−) sequence and not as isolated gauche bonds. Trans−gauche isomerization rates for the alkyl chains in the anchored bilayer are slower than those in lipid bilayers at the same temperature and show a progressive increase as the torsion numbers approach the tail. A two-dimensional periodic Voronoi tessellation analysis was performed to obtain the single-molecular area of an alkyl chain in the bilayer. The single-molecular area relaxation times are an order of magnitude longer than the trans−gauche isomerization times. The results indicate that the trans−gauche isomerization is associated with the creation and annihilation of a kink defect sequence. The results of the present MD simulation explain the apparent conflicting estimates of the gauche disorder in this system as obtained from infrared and 13C nuclear magnetic resonance measurements.

Development of a novel experimental model to investigate the influence of mechanics on bone healing

This dissertation proposed a novel experimental model combining a defect configuration with an active instrumented fixation device to investigate the influence of mechanics on bone healing. The proposed defect configuration aimed to minimise physiological loading within an experimental fracture gap and the instrumented fixator was used for the application of controlled displacements and in vivo stiffness monitoring of the healing process. This thesis has provided a novel approach to advance current knowledge and understanding of mechanobiology, which has been limited in previous experimental models.

Assessment of quality defects in macadamia kernels using NIR spectroscopy

Spectral data were collected of intact and ground kernels using 3 instruments (using Si-PbS, Si, and InGaAs detectors), operating over different areas of the spectrum (between 400 and 2500 nm) and employing transmittance, interactance, and reflectance sample presentation strategies. Kernels were assessed on the basis of oil and water content, and with respect to the defect categories of insect damage, rancidity, discoloration, mould growth, germination, and decomposition. Predictive model performance statistics for oil content models were acceptable on all instruments (R2 > 0.98; RMSECV < 2.5%, which is similar to reference analysis error), although that for the instrument employing reflectance optics was inferior to models developed for the instruments employing transmission optics. The spectral positions for calibration coefficients were consistent with absorbance due to the third overtones of CH2 stretching. Calibration models for moisture content in ground samples were acceptable on all instruments (R2 > 0.97; RMSECV < 0.2%), whereas calibration models for intact kernels were relatively poor. Calibration coefficients were more highly weighted around 1360, 740 and 840 nm, consistent with absorbance due to overtones of O-H stretching and combination. Intact kernels with brown centres or rancidity could be discriminated from each other and from sound kernels using principal component analysis. Part kernels affected by insect damage, discoloration, mould growth, germination, and decomposition could be discriminated from sound kernels. However, discrimination among these defect categories was not distinct and could not be validated on an independent set. It is concluded that there is good potential for a low cost Si photodiode array instrument to be employed to identify some quality defects of intact macadamia kernels and to quantify oil and moisture content of kernels in the process laboratory and for oil content in-line. Further work is required to examine the robustness of predictive models across different populations, including growing districts, cultivars and times of harvest.

Chemical bond approach to determining conductivity band gaps in amorphous chalcogenides and pnictides

The minimum energy required for the formation of conjugate pair of charged defects is found to be approximately equal to the experimental activation energy for d.c. conductivity in a number of amorphous chalcoganides and pnictides. This observation implies that the defect pair formation energy represents an intrinsic gap for transport in amorphous chalcogenides.

EPR study on the role of Mn in enhancing PTC of BaTiO3

Low concentration of Mn (< 0.05 atom%) added to lanthanide-doped ceramics for enhancing the PTC effect did not show any EPR signal due to Mn in the tetragonal phase. Above Tc (400 K) it showed the six-line signal arising from Mn2+. This is explained on the basis of Mn existing as Mn3+ ion with short relaxation time at room temperature. Oxidation state changes to Mn2+ above Tc; thus Mn3+ acts as an electron trap. This augments the function of activated defect centres (VBa /ag VBa) in diminishing the charge carrier concentration across the phase transformation.

Graphene-nanocrystalline metal sulphide composites produced by a one-pot reaction starting from graphite oxide

Graphene-nanocrystalline metal sulphide composites were prepared by a one-pot reaction. A dispersion of graphite oxide layers in an aqueous solution of metal ions (Cd2+/Zn2+) was reacted with H2S gas, which acts as a sulphide source as well as a reducing agent, resulting in the formation of metal sulphide nanoparticles and simultaneous reduction of graphite oxide sheets to graphene sheets. The surface defect related emissions shown by free metal sulphide particles are quenched in the composites due to the interaction of the surface of the nanoparticles with graphene sheets.

Optical properties of reduced LiNbO,

Optical characteristics of H2-reduced LiNbO3 have been studied. The optical transmission in the spectral range 0.3 to 5 mu m is progressively narrowed as the sample is reduced. An absorption band centred at about 2.48 eV is assigned to small polarons. This interpretation is consistent with the published defect structure, transport and spectroscopic reports on reduced LiNbO3. A semi-quantitative analysis of the 2.48 eV absorption band is carried out by employing the Reik and Heese theory (1967) for optical properties of small polarons.

DNase I Site Mapping and Micrococcal Nuclease Digestion of Pachytene Chromatin Reveal Novel Structural Features

A comparison of the DNase I digestion products of the 32P-5’-end-labeled pachytene nucleosome core particles (containing histones H2A, TH2A, X2, H2B, THPB, H3a, nd H4) and liver nucleosome core particles (containing somatic histones H2A, H2B, H3, and H4) revealed that the cleavage sites that are 30, 40, and 110 nucleotidesa way from the 5’-enda re significantly more accessiblei n the pachytene core particles than in the liver core particles. These cleavage sites correspond to the region wherein H2B interacts with the nucleosome core DNA. These results, therefore, suggest that the histone-DNA interactiona t these sites in the pachytene core particles is weaker, possibly because of the presence of the histone variant THBB interacting at similar topological positions in the nucleosome core as that of its somatic counterpart H2B. Such a loosened structumrea y also be maintainede ven in the native pachytene chromatin since micrococcal nuclease digestion of pachytene nuclei resulted in a higher ratio of subnucleosomes (SN4 + SN?) to mononucleosomes than that observed liinv er chromatin

XPES studies of oxides of second- and third-row transition metals including rare earths

X-ray photoelectron spectroscopy has been employed to investigate oxides of second- and third-row transition metals, including those of rare earths. Systematics in the spin—orbit splittings and binding energies of core levels of the metals are described. In most of the cases studied, the dependence of the spin—orbit splittings on the atomic number Z is given by the relation ΔE = a(Z - Z0)4, where a is the quantum defect parameter and Z0 is the effective screening. Core-level binding energies are found to increase with the oxidation state of the metal. Most of the core-level binding energies are related to the atomic number Z by the expression E = x(Z - Z0)2, giving rise to linear plots of ln E versus ln Z. Specific features of individual oxides, with respect to satellites, multiplet structure, configuration mixing, and other properties are also discussed. The spectra of PrO2, Pr6O11, TbO2 and Tb4O7 are reported for the first time.