Ecological and genetic interactions between cyanobacteria and viruses in a low-oxygen mat community inferred through metagenomics and metatranscriptomics.

Metagenomic and metatranscriptomic sequencing was conducted on cyanobacterial mats of the Middle Island Sinkhole (MIS), Lake Huron. Metagenomic data from 14 samples collected over 5 years were used to reconstruct genomes of two genotypes of a novel virus, designated PhV1 type A and PhV1 type B. Both viral genotypes encode and express nblA, a gene involved in degrading phycobilisomes, which are complexes of pigmented proteins that harvest light for photosynthesis. Phylogenetic analysis indicated that the viral-encoded nblA is derived from the host cyanobacterium, Phormidium MIS-PhA. The cyanobacterial host also has two complete CRISPR (clustered regularly interspaced short palindromic repeats) systems that serve as defence mechanisms for bacteria and archaea against viruses and plasmids. One 45 bp CRISPR spacer from Phormidium had 100% nucleotide identity to PhV1 type B, but this region was absent from PhV1 type A. Transcripts from PhV1 and the Phormidium CRISPR loci were detected in all six metatranscriptomic data sets (three during the day and three at night), indicating that both are transcriptionally active in the environment. These results reveal ecological and genetic interactions between viruses and cyanobacteria at MIS, highlighting the value of parallel analysis of viruses and hosts in understanding ecological interactions in natural communities.

LINE-1 Mediated Insertion into Poc1a (Protein of Centriole 1 A) Causes Growth Insufficiency and Male Infertility in Mice.

Skeletal dysplasias are a common, genetically heterogeneous cause of short stature that can result from disruptions in many cellular processes. We report the identification of the lesion responsible for skeletal dysplasia and male infertility in the spontaneous, recessive mouse mutant chagun. We determined that Poc1a, encoding protein of the centriole 1a, is disrupted by the insertion of a processed Cenpw cDNA, which is flanked by target site duplications, suggestive of a LINE-1 retrotransposon-mediated event. Mutant fibroblasts have impaired cilia formation and multipolar spindles. Male infertility is caused by defective spermatogenesis early in meiosis and progressive germ cell loss. Spermatogonial stem cell transplantation studies revealed that Poc1a is essential for normal function of both Sertoli cells and germ cells. The proliferative zone of the growth plate is small and disorganized because chondrocytes fail to re-align after cell division and undergo increased apoptosis. Poc1a and several other genes associated with centrosome function can affect the skeleton and lead to skeletal dysplasias and primordial dwarfisms. This mouse mutant reveals how centrosome dysfunction contributes to defects in skeletal growth and male infertility.

Genome-Wide Identification of Klebsiella pneumoniae Fitness Genes during Lung Infection.

UNLABELLED: Klebsiella pneumoniae is an urgent public health threat because of resistance to carbapenems, antibiotics of last resort against Gram-negative bacterial infections. Despite the fact that K. pneumoniae is a leading cause of pneumonia in hospitalized patients, the bacterial factors required to cause disease are poorly understood. Insertion site sequencing combines transposon mutagenesis with high-throughput sequencing to simultaneously screen thousands of insertion mutants for fitness defects during infection. Using the recently sequenced K. pneumoniae strain KPPR1 in a well-established mouse model of pneumonia, insertion site sequencing was performed on a pool of >25,000 transposon mutants. The relative fitness requirement of each gene was ranked based on the ratio of lung to inoculum read counts and concordance between insertions in the same gene. This analysis revealed over 300 mutants with at least a 2-fold fitness defect and 69 with defects ranging from 10- to >2,000-fold. Construction of 6 isogenic mutants for use in competitive infections with the wild type confirmed their requirement for lung fitness. Critical fitness genes included those for the synthesis of branched-chain and aromatic amino acids that are essential in mice and humans, the transcriptional elongation factor RfaH, and the copper efflux pump CopA. The majority of fitness genes were conserved among reference strains representative of diverse pathotypes. These results indicate that regulation of outer membrane components and synthesis of amino acids that are essential to its host are critical for K. pneumoniae fitness in the lung. IMPORTANCE: Klebsiella pneumoniae is a bacterium that commonly causes pneumonia in patients after they are admitted to the hospital. K. pneumoniae is becoming resistant to all available antibiotics, and when these infections spread to the bloodstream, over half of patients die. Since currently available antibiotics are failing, we must discover new ways to treat these infections. In this study, we asked what genes the bacterium needs to cause an infection, since the proteins encoded by these genes could be targets for new antibiotics. We identified over 300 genes that K. pneumoniae requires to grow in a mouse model of pneumonia. Many of the genes that we identified are found in K. pneumoniae isolates from throughout the world, including antibiotic-resistant forms. If new antibiotics could be made against the proteins that these genes encode, they may be broadly effective against K. pneumoniae.

Complete Genome Sequence of Klebsiella pneumoniae Strain ATCC 43816 KPPR1, a Rifampin-Resistant Mutant Commonly Used in Animal, Genetic, and Molecular Biology Studies.

Klebsiella pneumoniae is an urgent public health threat due to the spread of carbapenem-resistant strains causing serious, and frequently fatal, infections. To facilitate genetic, molecular, and immunological studies of this pathogen, we report the complete chromosomal sequence of a genetically tractable, prototypical strain used in animal models.

CD1d-unrestricted NKT cells are endowed with a hybrid function far superior than that of iNKT cells.

Invariant natural killer T (iNKT) cells to date represent the best example of cells known to have a hybrid function, representing both innate and adaptive immunity. Shared phenotypic similarities with NK cells together with a rapid response to a cytokine stimulus and a productive TCR engagement are the features that underline the hybrid nature of iNKT cells. Using these criteria, we provide molecular and functional evidence demonstrating that CD1d-independent (CD1d(ind)) NKT cells, a population of CD1d-unrestricted NKT cells, are endowed with a hybrid function far superior to that of iNKT cells: (i) an extensive shared program with NK cells, (ii) a closer Euclidian distance with NK cells, and (iii) the ability to respond to innate stimuli (Poly:IC) with cytotoxic potential in the same manner as NK cells identify a hybrid feature in CD1d(ind)NKT cells that truly fulfills the dual function of an NK and a T cell. Our finding that CD1d(ind)NKT cells are programmed to act like NK cells in response to innate signals while being capable of adaptive responses is unprecedented, and thus might reemphasize CD1d-unrestricted NKT cells as a subset of lymphocytes that could affect biological processes of antimicrobial and tumor immunity in a unique way.

Novel bioinformatics method for identification of genome-wide non-canonical spliced regions using RNA-Seq data.

SETTING: During endoplasmic reticulum (ER) stress, the endoribonuclease (RNase) Ire1α initiates removal of a 26 nt region from the mRNA encoding the transcription factor Xbp1 via an unconventional mechanism (atypically within the cytosol). This causes an open reading frame-shift that leads to altered transcriptional regulation of numerous downstream genes in response to ER stress as part of the unfolded protein response (UPR). Strikingly, other examples of targeted, unconventional splicing of short mRNA regions have yet to be reported. OBJECTIVE: Our goal was to develop an approach to identify non-canonical, possibly very short, splicing regions using RNA-Seq data and apply it to ER stress-induced Ire1α heterozygous and knockout mouse embryonic fibroblast (MEF) cell lines to identify additional Ire1α targets. RESULTS: We developed a bioinformatics approach called the Read-Split-Walk (RSW) pipeline, and evaluated it using two Ire1α heterozygous and two Ire1α-null samples. The 26 nt non-canonical splice site in Xbp1 was detected as the top hit by our RSW pipeline in heterozygous samples but not in the negative control Ire1α knockout samples. We compared the Xbp1 results from our approach with results using the alignment program BWA, Bowtie2, STAR, Exonerate and the Unix "grep" command. We then applied our RSW pipeline to RNA-Seq data from the SKBR3 human breast cancer cell line. RSW reported a large number of non-canonical spliced regions for 108 genes in chromosome 17, which were identified by an independent study. CONCLUSIONS: We conclude that our RSW pipeline is a practical approach for identifying non-canonical splice junction sites on a genome-wide level. We demonstrate that our pipeline can detect novel splice sites in RNA-Seq data generated under similar conditions for multiple species, in our case mouse and human.

Regulation of the human endogenous retrovirus K (HML-2) transcriptome by the HIV-1 Tat protein.

UNLABELLED: Approximately 8% of the human genome is made up of endogenous retroviral sequences. As the HIV-1 Tat protein activates the overall expression of the human endogenous retrovirus type K (HERV-K) (HML-2), we used next-generation sequencing to determine which of the 91 currently annotated HERV-K (HML-2) proviruses are regulated by Tat. Transcriptome sequencing of total RNA isolated from Tat- and vehicle-treated peripheral blood lymphocytes from a healthy donor showed that Tat significantly activates expression of 26 unique HERV-K (HML-2) proviruses, silences 12, and does not significantly alter the expression of the remaining proviruses. Quantitative reverse transcription-PCR validation of the sequencing data was performed on Tat-treated PBLs of seven donors using provirus-specific primers and corroborated the results with a substantial degree of quantitative similarity. IMPORTANCE: The expression of HERV-K (HML-2) is tightly regulated but becomes markedly increased following infection with HIV-1, in part due to the HIV-1 Tat protein. The findings reported here demonstrate the complexity of the genome-wide regulation of HERV-K (HML-2) expression by Tat. This work also demonstrates that although HERV-K (HML-2) proviruses in the human genome are highly similar in terms of DNA sequence, modulation of the expression of specific proviruses in a given biological situation can be ascertained using next-generation sequencing and bioinformatics analysis.

A bioinformatics approach reveals novel interactions of the OVOL transcription factors in the regulation of epithelial - mesenchymal cell reprogramming and cancer progression.

BACKGROUND: Mesenchymal to Epithelial Transition (MET) plasticity is critical to cancer progression, and we recently showed that the OVOL transcription factors (TFs) are critical regulators of MET. Results of that work also posed the hypothesis that the OVOLs impact MET in a range of cancers. We now test this hypothesis by developing a model, OVOL Induced MET (OI-MET), and sub-model (OI-MET-TF), to characterize differential gene expression in MET common to prostate cancer (PC) and breast cancer (BC). RESULTS: In the OI-MET model, we identified 739 genes differentially expressed in both the PC and BC models. For this gene set, we found significant enrichment of annotation for BC, PC, cancer, and MET, as well as regulation of gene expression by AP1, STAT1, STAT3, and NFKB1. Focusing on the target genes for these four TFs plus the OVOLs, we produced the OI-MET-TF sub-model, which shows even greater enrichment for these annotations, plus significant evidence of cooperation among these five TFs. Based on known gene/drug interactions, we prioritized targets in the OI-MET-TF network for follow-on analysis, emphasizing the clinical relevance of this work. Reflecting these results back to the OI-MET model, we found that binding motifs for the TF pair AP1/MYC are more frequent than expected and that the AP1/MYC pair is significantly enriched in binding in cancer models, relative to non-cancer models, in these promoters. This effect is seen in both MET models (solid tumors) and in non-MET models (leukemia). These results are consistent with our hypothesis that the OVOLs impact cancer susceptibility by regulating MET, and extend the hypothesis to include mechanisms not specific to MET. CONCLUSIONS: We find significant evidence of the OVOL, AP1, STAT1, STAT3, and NFKB1 TFs having important roles in MET, and more broadly in cancer. We prioritize known gene/drug targets for follow-up in the clinic, and we show that the AP1/MYC TF pair is a strong candidate for intervention.

H2A.B facilitates transcription elongation at methylated CpG loci.

H2A.B is a unique histone H2A variant that only exists in mammals. Here we found that H2A.B is ubiquitously expressed in major organs. Genome-wide analysis of H2A.B in mouse ES cells shows that H2A.B is associated with methylated DNA in gene body regions. Moreover, H2A.B-enriched gene loci are actively transcribed. One typical example is that H2A.B is enriched in a set of differentially methylated regions at imprinted loci and facilitates transcription elongation. These results suggest that H2A.B positively regulates transcription elongation by overcoming DNA methylation in the transcribed region. It provides a novel mechanism by which transcription is regulated at DNA hypermethylated regions.

Community transcriptomic assembly reveals microbes that contribute to deep-sea carbon and nitrogen cycling.

The deep ocean is an important component of global biogeochemical cycles because it contains one of the largest pools of reactive carbon and nitrogen on earth. However, the microbial communities that drive deep-sea geochemistry are vastly unexplored. Metatranscriptomics offers new windows into these communities, but it has been hampered by reliance on genome databases for interpretation. We reconstructed the transcriptomes of microbial populations from Guaymas Basin, in the deep Gulf of California, through shotgun sequencing and de novo assembly of total community RNA. Many of the resulting messenger RNA (mRNA) contiguous sequences contain multiple genes, reflecting co-transcription of operons, including those from dominant members. Also prevalent were transcripts with only limited representation (2.8 times coverage) in a corresponding metagenome, including a considerable portion (1.2 Mb total assembled mRNA sequence) with similarity (96%) to a marine heterotroph, Alteromonas macleodii. This Alteromonas and euryarchaeal marine group II populations displayed abundant transcripts from amino-acid transporters, suggesting recycling of organic carbon and nitrogen from amino acids. Also among the most abundant mRNAs were catalytic subunits of the nitrite oxidoreductase complex and electron transfer components involved in nitrite oxidation. These and other novel genes are related to novel Nitrospirae and have limited representation in accompanying metagenomic data. High throughput sequencing of 16S ribosomal RNA (rRNA) genes and rRNA read counts confirmed that Nitrospirae are minor yet widespread members of deep-sea communities. These results implicate a novel bacterial group in deep-sea nitrite oxidation, the second step of nitrification. This study highlights metatranscriptomic assembly as a valuable approach to study microbial communities.

Next-generation sequencing identifies the Danforth's short tail mouse mutation as a retrotransposon insertion affecting Ptf1a expression.

The semidominant Danforth's short tail (Sd) mutation arose spontaneously in the 1920s. The homozygous Sd phenotype includes severe malformations of the axial skeleton with an absent tail, kidney agenesis, anal atresia, and persistent cloaca. The Sd mutant phenotype mirrors features seen in human caudal malformation syndromes including urorectal septum malformation, caudal regression, VACTERL association, and persistent cloaca. The Sd mutation was previously mapped to a 0.9 cM region on mouse chromosome 2qA3. We performed Sanger sequencing of exons and intron/exon boundaries mapping to the Sd critical region and did not identify any mutations. We then performed DNA enrichment/capture followed by next-generation sequencing (NGS) of the critical genomic region. Standard bioinformatic analysis of paired-end sequence data did not reveal any causative mutations. Interrogation of reads that had been discarded because only a single end mapped correctly to the Sd locus identified an early transposon (ETn) retroviral insertion at the Sd locus, located 12.5 kb upstream of the Ptf1a gene. We show that Ptf1a expression is significantly upregulated in Sd mutant embryos at E9.5. The identification of the Sd mutation will lead to improved understanding of the developmental pathways that are misregulated in human caudal malformation syndromes.

Identification of cross-species shared transcriptional networks of diabetic nephropathy in human and mouse glomeruli.

Murine models are valuable instruments in defining the pathogenesis of diabetic nephropathy (DN), but they only partially recapitulate disease manifestations of human DN, limiting their utility. To define the molecular similarities and differences between human and murine DN, we performed a cross-species comparison of glomerular transcriptional networks. Glomerular gene expression was profiled in patients with early type 2 DN and in three mouse models (streptozotocin DBA/2, C57BLKS db/db, and eNOS-deficient C57BLKS db/db mice). Species-specific transcriptional networks were generated and compared with a novel network-matching algorithm. Three shared human-mouse cross-species glomerular transcriptional networks containing 143 (Human-DBA STZ), 97 (Human-BKS db/db), and 162 (Human-BKS eNOS(-/-) db/db) gene nodes were generated. Shared nodes across all networks reflected established pathogenic mechanisms of diabetes complications, such as elements of Janus kinase (JAK)/signal transducer and activator of transcription (STAT) and vascular endothelial growth factor receptor (VEGFR) signaling pathways. In addition, novel pathways not previously associated with DN and cross-species gene nodes and pathways unique to each of the human-mouse networks were discovered. The human-mouse shared glomerular transcriptional networks will assist DN researchers in selecting mouse models most relevant to the human disease process of interest. Moreover, they will allow identification of new pathways shared between mice and humans.

Disruption of RAB40AL function leads to Martin--Probst syndrome, a rare X-linked multisystem neurodevelopmental human disorder.

BACKGROUND AND AIM: Martin--Probst syndrome (MPS) is a rare X-linked disorder characterised by deafness, cognitive impairment, short stature and distinct craniofacial dysmorphisms, among other features. The authors sought to identify the causative mutation for MPS. METHODS AND RESULTS: Massively parallel sequencing in two affected, related male subjects with MPS identified a RAB40AL (also called RLGP) missense mutation (chrX:102,079,078-102,079,079AC→GA p.D59G; hg18). RAB40AL encodes a small Ras-like GTPase protein with one suppressor of cytokine signalling box. The p.D59G variant is located in a highly conserved region of the GTPase domain between ß-2 and ß-3 strands. Using RT-PCR, the authors show that RAB40AL is expressed in human fetal and adult brain and kidney, and adult lung, heart, liver and skeletal muscle. RAB40AL appears to be a primate innovation, with no orthologues found in mouse, Xenopus or zebrafish. Western analysis and fluorescence microscopy of GFP-tagged RAB40AL constructs from transiently transfected COS7 cells show that the D59G missense change renders RAB40AL unstable and disrupts its cytoplasmic localisation. CONCLUSIONS: This is the first study to show that mutation of RAB40AL is associated with a human disorder. Identification of RAB40AL as the gene mutated in MPS allows for further investigations into the molecular mechanism(s) of RAB40AL and its roles in diverse processes such as cognition, hearing and skeletal development.

Convergence of genetic influences in comorbidity.

BACKGROUND: Predisposition to complex diseases is explained in part by genetic variation, and complex diseases are frequently comorbid, consistent with pleiotropic genetic variation influencing comorbidity. Genome Wide Association (GWA) studies typically assess association between SNPs and a single-disease phenotype. Fisher meta-analysis combines evidence of association from single-disease GWA studies, assuming that each study is an independent test of the same hypothesis. The Rank Product (RP) method overcomes limitations posed by Fisher assumptions, though RP was not designed for GWA data. METHODS: We modified RP to accommodate GWA data, and we call it modRP. Using p-values output from GWA studies, we aggregate evidence for association between SNPs and related phenotypes. To assess significance, RP randomly samples the observed ranks to develop the null distribution of the RP statistic, and then places the observed RPs into the null distribution. ModRP eliminates the effect of linkage disequilibrium and controls for differences in power at tested SNPs, to meet RP assumptions in application to GWA data. RESULTS: After validating modRP based on both positive and negative control studies, we searched for pleiotropic influences on comorbid substance use disorders in a novel study, and found two SNPs to be significantly associated with comorbid cocaine, opium, and nicotine dependence. Placing these SNPs into biological context, we developed a protein network modeling the interaction of cocaine, nicotine, and opium with these variants. CONCLUSIONS: ModRP is a novel approach to identifying pleiotropic genetic influences on comorbid complex diseases. It can be used to assess association for related phenotypes where raw data is unavailable or inappropriate for analysis using other approaches. The method is conceptually simple and produces statistically significant, biologically relevant results.

Decade-long bacterial community dynamics in cystic fibrosis airways.

The structure and dynamics of bacterial communities in the airways of persons with cystic fibrosis (CF) remain largely unknown. We characterized the bacterial communities in 126 sputum samples representing serial collections spanning 8-9 y from six age-matched male CF patients. Sputum DNA was analyzed by bar-coded pyrosequencing of the V3-V5 hypervariable region of the 16S rRNA gene, defining 662 operational taxonomic units (OTUs) from >633,000 sequences. Bacterial community diversity decreased significantly over time in patients with typically progressive lung disease but remained relatively stable in patients with a mild lung disease phenotype. Antibiotic use, rather than patient age or lung function, was the primary driver of decreasing diversity. Interpatient variability in community structure exceeded intrapatient variability in serial samples. Antibiotic treatment was associated with pronounced shifts in community structure, but communities showed both short- and long-term resilience after antibiotic perturbation. There was a positive correlation between OTU occurrence and relative abundance, with a small number of persistent OTUs accounting for the greatest abundance. Significant changes in community structure, diversity, or total bacterial density at the time of pulmonary exacerbation were not observed. Despite decreasing community diversity in patients with progressive disease, total bacterial density remained relatively stable over time. These findings show the critical relationship between airway bacterial community structure, disease stage, and clinical state at the time of sample collection. These features are the key parameters with which to assess the complex ecology of the CF airway.

The NIH National Center for Integrative Biomedical Informatics (NCIBI).

The National Center for Integrative and Biomedical Informatics (NCIBI) is one of the eight NCBCs. NCIBI supports information access and data analysis for biomedical researchers, enabling them to build computational and knowledge models of biological systems to address the Driving Biological Problems (DBPs). The NCIBI DBPs have included prostate cancer progression, organ-specific complications of type 1 and 2 diabetes, bipolar disorder, and metabolic analysis of obesity syndrome. Collaborating with these and other partners, NCIBI has developed a series of software tools for exploratory analysis, concept visualization, and literature searches, as well as core database and web services resources. Many of our training and outreach initiatives have been in collaboration with the Research Centers at Minority Institutions (RCMI), integrating NCIBI and RCMI faculty and students, culminating each year in an annual workshop. Our future directions include focusing on the TranSMART data sharing and analysis initiative.

A genetic network model of cellular responses to lithium treatment and cocaine abuse in bipolar disorder.

BACKGROUND: Lithium is an effective treatment for Bipolar Disorder (BD) and significantly reduces suicide risk, though the molecular basis of lithium's effectiveness is not well understood. We seek to improve our understanding of this effectiveness by posing hypotheses based on new experimental data as well as published data, testing these hypotheses in silico, and posing new hypotheses for validation in future studies. We initially hypothesized a gene-by-environment interaction where lithium, acting as an environmental influence, impacts signal transduction pathways leading to differential expression of genes important in the etiology of BD mania. RESULTS: Using microarray and rt-QPCR assays, we identified candidate genes that are differentially expressed with lithium treatment. We used a systems biology approach to identify interactions among these candidate genes and develop a network of genes that interact with the differentially expressed candidates. Notably, we also identified cocaine as having a potential influence on the network, consistent with the observed high rate of comorbidity for BD and cocaine abuse. The resulting network represents a novel hypothesis on how multiple genetic influences on bipolar disorder are impacted by both lithium treatment and cocaine use. Testing this network for association with BD and related phenotypes, we find that it is significantly over-represented for genes that participate in signal transduction, consistent with our hypothesized-gene-by environment interaction. In addition, it models related pharmacogenomic, psychiatric, and chemical dependence phenotypes. CONCLUSIONS: We offer a network model of gene-by-environment interaction associated with lithium's effectiveness in treating BD mania, as well as the observed high rate of comorbidity of BD and cocaine abuse. We identified drug targets within this network that represent immediate candidates for therapeutic drug testing. Posing novel hypotheses for validation in future work, we prioritized SNPs near genes in the network based on functional annotation. We also developed a "concept signature" for the genes in the network and identified additional candidate genes that may influence the system because they are significantly associated with the signature.

Candidate exome capture identifies mutation of SDCCAG8 as the cause of a retinal-renal ciliopathy.

Nephronophthisis-related ciliopathies (NPHP-RC) are recessive disorders that feature dysplasia or degeneration occurring preferentially in the kidney, retina and cerebellum. Here we combined homozygosity mapping with candidate gene analysis by performing 'ciliopathy candidate exome capture' followed by massively parallel sequencing. We identified 12 different truncating mutations of SDCCAG8 (serologically defined colon cancer antigen 8, also known as CCCAP) in 10 families affected by NPHP-RC. We show that SDCCAG8 is localized at both centrioles and interacts directly with OFD1 (oral-facial-digital syndrome 1), which is associated with NPHP-RC. Depletion of sdccag8 causes kidney cysts and a body axis defect in zebrafish and induces cell polarity defects in three-dimensional renal cell cultures. This work identifies loss of SDCCAG8 function as a cause of a retinal-renal ciliopathy and validates exome capture analysis for broadly heterogeneous single-gene disorders.

Metscape: a Cytoscape plug-in for visualizing and interpreting metabolomic data in the context of human metabolic networks.

SUMMARY: Metscape is a plug-in for Cytoscape, used to visualize and interpret metabolomic data in the context of human metabolic networks. We have developed a metabolite database by extracting and integrating information from several public sources. By querying this database, Metscape allows users to trace the connections between metabolites and genes, visualize compound networks and display compound structures as well as information for reactions, enzymes, genes and pathways. Applying the pathway filter, users can create subnetworks that consist of compounds and reactions from a given pathway. Metscape allows users to upload experimental data, and visualize and explore compound networks over time, or experimental conditions. Color and size of the nodes are used to visualize these dynamic changes. Metscape can display the entire metabolic network or any of the pathway-specific networks that exist in the database. AVAILABILITY: Metscape can be installed from within Cytoscape 2.6.x under 'Network and Attribute I/O' category. For more information, please visit http://metscape.ncibi.org/tryplugin.html.

Modeling complex genetic and environmental influences on comorbid bipolar disorder with tobacco use disorder.

BACKGROUND: Comorbidity of psychiatric and substance use disorders represents a significant complication in the clinical course of both disorders. Bipolar Disorder (BD) is a psychiatric disorder characterized by severe mood swings, ranging from mania to depression, and up to a 70% rate of comorbid Tobacco Use Disorder (TUD). We found epidemiological evidence consistent with a common underlying etiology for BD and TUD, as well as evidence of both genetic and environmental influences on BD and TUD. Therefore, we hypothesized a common underlying genetic etiology, interacting with nicotine exposure, influencing susceptibility to both BD and TUD. METHODS: Using meta-analysis, we compared TUD rates for BD patients and the general population. We identified candidate genes showing statistically significant, replicated, evidence of association with both BD and TUD. We assessed commonality among these candidate genes and hypothesized broader, multi-gene network influences on the comorbidity. Using Fisher Exact tests we tested our hypothesized genetic networks for association with the comorbidity, then compared the inferences drawn with those derived from the commonality assessment. Finally, we prioritized candidate SNPs for validation. RESULTS: We estimate risk for TUD among BD patients at 2.4 times that of the general population. We found three candidate genes associated with both BD and TUD (COMT, SLC6A3, and SLC6A4) and commonality analysis suggests that these genes interact in predisposing psychiatric and substance use disorders. We identified a 69 gene network that influences neurotransmitter signaling and shows significant over-representation of genes associated with BD and TUD, as well as genes differentially expressed with exposure to tobacco smoke. Twenty four of these genes are known drug targets. CONCLUSIONS: This work highlights novel bioinformatics resources and demonstrates the effectiveness of using an integrated bioinformatics approach to improve our understanding of complex disease etiology. We illustrate the development and testing of hypotheses for a comorbidity predisposed by both genetic and environmental influences. Consistent with our hypothesis, the selected network models multiple interacting genetic influences on comorbid BD with TUD, as well as the environmental influence of nicotine. This network nominates candidate genes for validation and drug testing, and we offer a panel of SNPs prioritized for follow-up.