An Introductory Guide to Aligning Networks Using SANA, the Simulated Annealing Network Aligner.

Sequence alignment has had an enormous impact on our understanding of biology, evolution, and disease. The alignment of biological networks holds similar promise. Biological networks generally model interactions between biomolecules such as proteins, genes, metabolites, or mRNAs. There is strong evidence that the network topology-the "structure" of the network-is correlated with the functions performed, so that network topology can be used to help predict or understand function. However, unlike sequence comparison and alignment-which is an essentially solved problem-network comparison and alignment is an NP-complete problem for which heuristic algorithms must be used.Here we introduce SANA, the Simulated Annealing Network Aligner. SANA is one of many algorithms proposed for the arena of biological network alignment. In the context of global network alignment, SANA stands out for its speed, memory efficiency, ease-of-use, and flexibility in the arena of producing alignments between two or more networks. SANA produces better alignments in minutes on a laptop than most other algorithms can produce in hours or days of CPU time on large server-class machines. We walk the user through how to use SANA for several types of biomolecular networks

High quality copy number and genotype data from FFPE samples using Molecular Inversion Probe (MIP) microarrays

BACKGROUND:A major challenge facing DNA copy number (CN) studies of tumors is that most banked samples with extensive clinical follow-up information are Formalin-Fixed Paraffin Embedded (FFPE). DNA from FFPE samples generally underperforms or suffers high failure rates compared to fresh frozen samples because of DNA degradation and cross-linking during FFPE fixation and processing. As FFPE protocols may vary widely between labs and samples may be stored for decades at room temperature, an ideal FFPE CN technology should work on diverse sample sets. Molecular Inversion Probe (MIP) technology has been applied successfully to obtain high quality CN and genotype data from cell line and frozen tumor DNA. Since the MIP probes require only a small (~40 bp) target binding site, we reasoned they may be well suited to assess degraded FFPE DNA. We assessed CN with a MIP panel of 50,000 markers in 93 FFPE tumor samples from 7 diverse collections. For 38 FFPE samples from three collections we were also able to asses CN in matched fresh frozen tumor tissue.RESULTS:Using an input of 37 ng genomic DNA, we generated high quality CN data with MIP technology in 88% of FFPE samples from seven diverse collections. When matched fresh frozen tissue was available, the performance of FFPE DNA was comparable to that of DNA obtained from matched frozen tumor (genotype concordance averaged 99.9%), with only a modest loss in performance in FFPE.CONCLUSION:MIP technology can be used to generate high quality CN and genotype data in FFPE as well as fresh frozen samples.This item is part of the UA Faculty Publications collection. For more information this item or other items in the UA Campus Repository, contact the University of Arizona Libraries at [email protected]

Prediction of epigenetically regulated genes in breast cancer cell lines

Methylation of CpG islands within the DNA promoter regions is one mechanism that leads to aberrant gene expression in cancer. In particular, the abnormal methylation of CpG islands may silence associated genes. Therefore, using high-throughput microarrays to measure CpG island methylation will lead to better understanding of tumor pathobiology and progression, while revealing potentially new biomarkers. We have examined a recently developed high-throughput technology for measuring genome-wide methylation patterns called mTACL. Here, we propose a computational pipeline for integrating gene expression and CpG island methylation profles to identify epigenetically regulated genes for a panel of 45 breast cancer cell lines, which is widely used in the Integrative Cancer Biology Program (ICBP). The pipeline (i) reduces the dimensionality of the methylation data, (ii) associates the reduced methylation data with gene expression data, and (iii) ranks methylation-expression associations according to their epigenetic regulation. Dimensionality reduction is performed in two steps: (i) methylation sites are grouped across the genome to identify regions of interest, and (ii) methylation profles are clustered within each region. Associations between the clustered methylation and the gene expression data sets generate candidate matches within a fxed neighborhood around each gene. Finally, the methylation-expression associations are ranked through a logistic regression, and their significance is quantified through permutation analysis. Our two-step dimensionality reduction compressed 90% of the original data, reducing 137,688 methylation sites to 14,505 clusters. Methylation-expression associations produced 18,312 correspondences, which were used to further analyze epigenetic regulation. Logistic regression was used to identify 58 genes from these correspondences that showed a statistically signifcant negative correlation between methylation profles and gene expression in the panel of breast cancer cell lines. Subnetwork enrichment of these genes has identifed 35 common regulators with 6 or more predicted markers. In addition to identifying epigenetically regulated genes, we show evidence of differentially expressed methylation patterns between the basal and luminal subtypes. Our results indicate that the proposed computational protocol is a viable platform for identifying epigenetically regulated genes. Our protocol has generated a list of predictors including COL1A2, TOP2A, TFF1, and VAV3, genes whose key roles in epigenetic regulation is documented in the literature. Subnetwork enrichment of these predicted markers further suggests that epigenetic regulation of individual genes occurs in a coordinated fashion and through common regulators

Heterodimers of photoreceptor-specific nuclear receptor (PNR/NR2E3) and peroxisome proliferator-activated receptor (PPARγ) are disrupted by retinal disease-associated mutations

Photoreceptor-specific nuclear receptor (PNR/NR2E3) and Tailless homolog (TLX/NR2E1) are human orthologs of the NR2E group, a subgroup of phylogenetically related members of the Nuclear Receptor (NR) superfamily of transcription factors. We assessed the ability of these NRs to form heterodimers with other members of the human NRs representing all major subgroups. The TLX ligand binding domain (LBD) did not appear to form homodimers or interact directly with any other NR tested. The PNR LBD was able to form homodimers, but also exhibited robust interactions with the LBDs of PPARγ/NR1C3 and TRβ/NR1A2. The binding of PNR to PPARγ was specific for this paralog, as no interaction was observed with the LBDs of PPARαNR1C1 or PPARδNR1C2. In support of these findings, PPARγ and PNR were found to be co-expressed in human retinal tissue extracts and could be co-immunoprecipitated as a native complex. Selected sequence variants in the PNR LBD associated with human retinopathies, or a mutation in the dimerization region of PPARγ LBD associated with familial partial lipodystrophy type 3, were found to disrupt PNR/PPARγ complex formation. Wild type PNR, but not a PNR309G mutant, was able to repress PPARγ-mediated transcription in reporter assays. In summary our results reveal novel heterodimer interactions in the NR superfamily, suggesting previously unknown functional interactions of PNR with PPARγ and TRβ that have potential importance in retinal development and disease

Exploring Uncoupling Proteins and Antioxidant Mechanisms under Acute Cold Exposure in Brains of Fish

Exposure to fluctuating temperatures accelerates the mitochondrial respiration and increases the formation of mitochondrial reactive oxygen species (ROS) in ectothermic vertebrates including fish. To date, little is known on potential oxidative damage and on protective antioxidative defense mechanisms in the brain of fish under cold shock. In this study, the concentration of cellular protein carbonyls in brain was significantly increased by 38% within 1 h after cold exposure (from 28°C to 18°C) of zebrafish (Danio rerio). In addition, the specific activity of superoxide dismutase (SOD) and the mRNA level of catalase (CAT) were increased after cold exposure by about 60% (6 h) and by 60%–90% (1 and 24 h), respectively, while the specific glutathione content as well as the ratio of glutathione disulfide to glutathione remained constant and at a very low level. In addition, cold exposure increased the protein level of hypoxia-inducible factor (HIF) by about 50% and the mRNA level of the glucose transporter zglut3 in brain by 50%–100%. To test for an involvement of uncoupling proteins (UCPs) in the cold adaptation of zebrafish, five UCP members were annotated and identified (zucp1-5). With the exception of zucp1, the mRNA levels of the other four zucps were significantly increased after cold exposure. In addition, the mRNA levels of four of the fish homologs (zppar) of the peroxisome proliferator-activated receptor (PPAR) were increased after cold exposure. These data suggest that PPARs and UCPs are involved in the alterations observed in zebrafish brain after exposure to 18°C. The observed stimulation of the PPAR-UCP axis may help to prevent oxidative damage and to maintain metabolic balance and cellular homeostasis in the brains of ectothermic zebrafish upon cold exposure

Robust simplifications of multiscale biochemical networks

<p>Abstract</p> <p>Background</p> <p>Cellular processes such as metabolism, decision making in development and differentiation, signalling, etc., can be modeled as large networks of biochemical reactions. In order to understand the functioning of these systems, there is a strong need for general model reduction techniques allowing to simplify models without loosing their main properties. In systems biology we also need to compare models or to couple them as parts of larger models. In these situations reduction to a common level of complexity is needed.</p> <p>Results</p> <p>We propose a systematic treatment of model reduction of multiscale biochemical networks. First, we consider linear kinetic models, which appear as "pseudo-monomolecular" subsystems of multiscale nonlinear reaction networks. For such linear models, we propose a reduction algorithm which is based on a generalized theory of the limiting step that we have developed in <abbrgrp><abbr bid="B1">1</abbr></abbrgrp>. Second, for non-linear systems we develop an algorithm based on dominant solutions of quasi-stationarity equations. For oscillating systems, quasi-stationarity and averaging are combined to eliminate time scales much faster and much slower than the period of the oscillations. In all cases, we obtain robust simplifications and also identify the critical parameters of the model. The methods are demonstrated for simple examples and for a more complex model of NF-<it>κ</it>B pathway.</p> <p>Conclusion</p> <p>Our approach allows critical parameter identification and produces hierarchies of models. Hierarchical modeling is important in "middle-out" approaches when there is need to zoom in and out several levels of complexity. Critical parameter identification is an important issue in systems biology with potential applications to biological control and therapeutics. Our approach also deals naturally with the presence of multiple time scales, which is a general property of systems biology models.</p

Morphological and Pathological Evolution of the Brain Microcirculation in Aging and Alzheimer’s Disease

Key pathological hallmarks of Alzheimer’s disease (AD), including amyloid plaques, cerebral amyloid angiopathy (CAA) and neurofibrillary tangles do not completely account for cognitive impairment, therefore other factors such as cardiovascular and cerebrovascular pathologies, may contribute to AD. In order to elucidate the microvascular changes that contribute to aging and disease, direct neuropathological staining and immunohistochemistry, were used to quantify the structural integrity of the microvasculature and its innervation in three oldest-old cohorts: 1) nonagenarians with AD and a high amyloid plaque load; 2) nonagenarians with no dementia and a high amyloid plaque load; 3) nonagenarians without dementia or amyloid plaques. In addition, a non-demented (ND) group (average age 71 years) with no amyloid plaques was included for comparison. While gray matter thickness and overall brain mass were reduced in AD compared to ND control groups, overall capillary density was not different. However, degenerated string capillaries were elevated in AD, potentially suggesting greater microvascular “dysfunction” compared to ND groups. Intriguingly, apolipoprotein ε4 carriers had significantly higher string vessel counts relative to non-ε4 carriers. Taken together, these data suggest a concomitant loss of functional capillaries and brain volume in AD subjects. We also demonstrated a trend of decreasing vesicular acetylcholine transporter staining, a marker of cortical cholinergic afferents that contribute to arteriolar vasoregulation, in AD compared to ND control groups, suggesting impaired control of vasodilation in AD subjects. In addition, tyrosine hydroxylase, a marker of noradrenergic vascular innervation, was reduced which may also contribute to a loss of control of vasoconstriction. The data highlight the importance of the brain microcirculation in the pathogenesis and evolution of AD

Search for neutral MSSM Higgs bosons decaying to a pair of tau leptons in pp collisions

Peer reviewe

An introductory guide to aligning networks using SANA, the Simulated Annealing Network Aligner

Sequence alignment has had an enormous impact on our understanding of biology, evolution, and disease. The alignment of biological {\em networks} holds similar promise. Biological networks generally model interactions between biomolecules such as proteins, genes, metabolites, or mRNAs. There is strong evidence that the network topology -- the "structure" of the network -- is correlated with the functions performed, so that network topology can be used to help predict or understand function. However, unlike sequence comparison and alignment -- which is an essentially solved problem -- network comparison and alignment is an NP-complete problem for which heuristic algorithms must be used. Here we introduce SANA, the {\it Simulated Annealing Network Aligner}. SANA is one of many algorithms proposed for the arena of biological network alignment. In the context of global network alignment, SANA stands out for its speed, memory efficiency, ease-of-use, and flexibility in the arena of producing alignments between 2 or more networks. SANA produces better alignments in minutes on a laptop than most other algorithms can produce in hours or days of CPU time on large server-class machines. We walk the user through how to use SANA for several types of biomolecular networks. Availability: https://github.com/waynebhayes/SAN

Searches for electroweak production of charginos, neutralinos, and sleptons decaying to leptons and W, Z, and Higgs bosons in pp collisions at 8 TeV

Peer reviewe