Parallel biocomputing

<p>Abstract</p> <p>Background</p> <p>With the advent of high throughput genomics and high-resolution imaging techniques, there is a growing necessity in biology and medicine for parallel computing, and with the low cost of computing, it is now cost-effective for even small labs or individuals to build their own personal computation cluster.</p> <p>Methods</p> <p>Here we briefly describe how to use commodity hardware to build a low-cost, high-performance compute cluster, and provide an in-depth example and sample code for parallel execution of R jobs using MOSIX, a mature extension of the Linux kernel for parallel computing. A similar process can be used with other cluster platform software.</p> <p>Results</p> <p>As a statistical genetics example, we use our cluster to run a simulated eQTL experiment. Because eQTL is computationally intensive, and is conceptually easy to parallelize, like many statistics/genetics applications, parallel execution with MOSIX gives a linear speedup in analysis time with little additional effort.</p> <p>Conclusions</p> <p>We have used MOSIX to run a wide variety of software programs in parallel with good results. The limitations and benefits of using MOSIX are discussed and compared to other platforms.</p

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

Co-regulatory expression quantitative trait loci mapping: method and application to endometrial cancer

<p>Abstract</p> <p>Background</p> <p>Expression quantitative trait loci (eQTL) studies have helped identify the genetic determinants of gene expression. Understanding the potential interacting mechanisms underlying such findings, however, is challenging.</p> <p>Methods</p> <p>We describe a method to identify the <it>trans-</it>acting drivers of multiple gene co-expression, which reflects the action of regulatory molecules. This method-termed <it>co-regulatory expression quantitative trait locus </it>(creQTL) <it>mapping</it>-allows for evaluation of a more focused set of phenotypes within a clear biological context than conventional eQTL mapping.</p> <p>Results</p> <p>Applying this method to a study of endometrial cancer revealed regulatory mechanisms supported by the literature: a creQTL between a locus upstream of STARD13/DLC2 and a group of seven IFNβ-induced genes. This suggests that the Rho-GTPase encoded by STARD13 regulates IFNβ-induced genes and the DNA damage response.</p> <p>Conclusions</p> <p>Because of the importance of IFNβ in cancer, our results suggest that creQTL may provide a finer picture of gene regulation and may reveal additional molecular targets for intervention. An open source R implementation of the method is available at <url>http://sites.google.com/site/kenkompass/</url>.</p

Trabecular Meshwork Gene Expression after Selective Laser Trabeculoplasty

BACKGROUND: Trabecular meshwork and Schlemm's canal are the tissues appointed to modulate the aqueous humour outflow from the anterior chamber. The impairment of their functions drives to an intraocular pressure increase. The selective laser trabeculoplasty is a laser therapy of the trabecular meshwork able to decrease intraocular pressure. The exact response mechanism to this treatment has not been clearly delineated yet. The herein presented study is aimed at studying the gene expression changes induced in trabecular meshwork cells by selective laser trabeculoplasty (SLT) in order to better understand the mechanisms subtending its efficacy. METHODOLOGY/PRINCIPAL FINDINGS: Primary human trabecular meshwork cells cultured in fibroblast medium underwent selective laser trabeculoplasty treatment. RNA was extracted from a pool of cells 30 minutes after treatment while the remaining cells were further cultured and RNA was extracted respectively 2 and 6 hours after treatment. Control cells stored in incubator in absence of SLT treatment were used as reference samples. Gene expression was evaluated by hybridization on miRNA-microarray and laser scanner analysis. Scanning electron microscopic examination was performed on 2 Trabecular meshwork samples after SLT at 4(th) and 6(th) hour from treatment. On the whole, selective laser trabeculoplasty modulates in trabecular meshwork the expression of genes involved in cell motility, intercellular connections, extracellular matrix production, protein repair, DNA repair, membrane repair, reactive oxygen species production, glutamate toxicity, antioxidant activities, and inflammation. CONCLUSIONS/SIGNIFICANCE: SLT did not induce any phenotypic alteration in TM samples. TM is a complex tissue possessing a great variety of function pivotal for the active regulation of aqueous humour outflow from the anterior chamber. SLT is able to modulate these functions at the postgenomic molecular level without inducing damage either at molecular or phenotypic levels

Bioinformatic and statistical analysis of the optic nerve head in a primate model of ocular hypertension

<p>Abstract</p> <p>Background</p> <p>The nonhuman primate model of glaucomatous optic neuropathy most faithfully reproduces the human disease. We used high-density oligonucleotide arrays to investigate whole genome transcriptional changes occurring at the optic nerve head during primate experimental glaucoma.</p> <p>Results</p> <p>Laser scarification of the trabecular meshwork of cynomolgus macaques produced elevated intraocular pressure that was monitored over time and led to varying degrees of damage in different samples. The macaques were examined clinically before enucleation and the myelinated optic nerves were processed post-mortem to determine the degree of neuronal loss. Global gene expression was examined in dissected optic nerve heads with Affymetrix GeneChip microarrays. We validated a subset of differentially expressed genes using qRT-PCR, immunohistochemistry, and immuno-enriched astrocytes from healthy and glaucomatous human donors. These genes have previously defined roles in axonal outgrowth, immune response, cell motility, neuroprotection, and extracellular matrix remodeling.</p> <p>Conclusion</p> <p>Our findings show that glaucoma is associated with increased expression of genes that mediate axonal outgrowth, immune response, cell motility, neuroprotection, and ECM remodeling. These studies also reveal that, as glaucoma progresses, retinal ganglion cell axons may make a regenerative attempt to restore lost nerve cell contact.</p