sscMap: An extensible Java application for connecting small-molecule drugs using gene-expression signatures

Background
Connectivity mapping is a process to recognize novel pharmacological and toxicological properties in small molecules by comparing their gene expression signatures with others in a database. A simple and robust method for connectivity mapping with increased specificity and sensitivity was recently developed, and its utility demonstrated using experimentally derived gene signatures.

Results
This paper introduces sscMap (statistically significant connections' map), a Java application designed to undertake connectivity mapping tasks using the recently published method. The software is bundled with a default collection of reference gene-expression profiles based on the publicly available dataset from the Broad Institute Connectivity Map 02, which includes data from over 7000 Affymetrix microarrays, for over 1000 small-molecule compounds, and 6100 treatment instances in 5 human cell lines. In addition, the application allows users to add their custom collections of reference profiles and is applicable to a wide range of other 'omics technologies.

Conclusion
The utility of sscMap is two fold. First, it serves to make statistically significant connections between a user-supplied gene signature and the 6100 core reference profiles based on the Broad Institute expanded dataset. Second, it allows users to apply the same improved method to custom-built reference profiles which can be added to the database for future referencing. The software can be freely downloaded from http://purl.oclc.org/NET/sscMap

Loss-of-function mutations in the cathepsin C gene result in periodontal disease and palmoplantar keratosis

Papillon-Lefevre syndrome, or keratosis palmoplantaris with periodontopathia (PLS, MIM 245000), is an autosomal recessive disorder that is mainly ascertained by dentists because of the severe periodontitis that afflicts patients(1,2). Both the deciduous and permanent dentitions are affected, resulting in premature tooth loss. Palmoplantar keratosis, varying from mild psoriasiform scaly skin to overt hyperkeratosis, typically develops within the first three years of life. Keratosis also affects other sites such as elbows and knees. Most PLS patients display both periodontitis and hyperkeratosis. some patients have only palmoplantar keratosis or periodontitis, and in rare individuals the periodontitis is mild and of late onset(3-6). The PLS locus has been mapped to chromosome 11q14-q21 (refs 7-9). Using homozygosity mapping in eight small consanguineous families, we have narrowed the candidate region to a 1.2-cM interval between D11S4082 and D11S931. The gene (CTSC) encoding the lysosomal protease cathepsin C (or dipeptidyl aminopeptidase I) lies within this interval. We defined the genomic structure of CTSC and found mutations in all eight families. In two of these families we used a functional assay to demonstrate an almost total loss of cathepsin C activity in PLS patients and reduced activity in obligate carriers.

Rice-arsenate interactions in hydroponics: a three-gene model for tolerance

In this study, the genetic mapping of the tolerance of root growth to 13.3 muM arsenate [As(V)] using the BalaxAzucena population is improved, and candidate genes for further study are identified. A remarkable three-gene model of tolerance is advanced, which appears to involve epistatic interaction between three major genes, two on chromosome 6 and one on chromosome 10. Any combination of two of these genes inherited from the tolerant parent leads to the plant having tolerance. Lists of potential positional candidate genes are presented. These are then refined using whole genome transcriptomics data and bioinformatics. Physiological evidence is also provided that genes related to phosphate transport are unlikely to be behind the genetic loci conferring tolerance. These results offer testable hypotheses for genes related to As(V) tolerance that might offer strategies for mitigating arsenic (As) accumulation in consumed rice.

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

Exclusion of linkage between schizophrenia and the D2 dopamine receptor gene region of chromosome 11q in 112 Irish multiplex families

A leading theory hypothesizes that schizophrenia arises from dysregulation of the dopamine system in certain brain regions. As this dysregulation could arise from abnormal expression of D2 dopamine receptors, the D2 receptor gene (DRD2) on chromosome 11q is a candidate locus for schizophrenia. We tested whether allelic variation at DRD2 and five surrounding loci cosegregated with schizophrenia in 112 small- to moderate-size Irish families containing two or more members affected with schizophrenia or schizoaffective disorder, defined by DSM-III-R. Evidence of linkage was assessed using varying definitions of illness and modes of transmission. Assuming genetic homogeneity, linkage between schizophrenia and large regions of 11q around DRD2 could be strongly excluded. Assuming genetic heterogeneity, variation at the DRD2 locus could be rejected as a major risk factor for schizophrenia in more than 50% of these families for all models tested and in as few as 25% of the families for certain models. The DRD2 linkage in fewer than 25% of these families could not be excluded under any of the models tested. Our results suggest that the major component of genetic susceptibility to schizophrenia is not due to allelic variation at the DRD2 locus or other genes in the surrounding chromosomal region.

Genetic variation in the 6p22.3 gene DTNBP1, the human ortholog of the mouse dysbindin gene, is associated with schizophrenia

Prior evidence has supported the existence of multiple susceptibility genes for schizophrenia. Multipoint linkage analysis of the 270 Irish high-density pedigrees that we have studied, as well as results from several other samples, suggest that at least one such gene is located in region 6p24-21. In the present study, family-based association analysis of 36 simple sequence-length-polymorphism markers and of 17 SNP markers implicated two regions, separated by approximately 7 Mb. The first region, and the focus of this report, is 6p22.3. In this region, single-nucleotide polymorphisms within the 140-kb gene DTNBP1 (dystrobrevin-binding protein 1, or dysbindin) are strongly associated with schizophrenia. Uncorrected, empirical P values produced by the program TRANSMIT were significant (P

Inactivating mutations in an SH2 domain-encoding gene in X-linked lymphoproliferative syndrome

X-linked lymphoproliferative syndrome (XLP) is an inherited immunodeficiency characterized by increased susceptibility to Epstein-Barr virus (EBV). In affected males, primary EBV infection leads to the uncontrolled proliferation of virus-containing B cells and reactive cytotoxic T cells, often culminating in the development of high-grade lymphoma. The XLP gene has been mapped to chromosome band Xq25 through linkage analysis and the discovery of patients harboring large constitutional genomic deletions. We describe here the presence of small deletions and intragenic mutations that specifically disrupt a gene named DSHP in 6 of 10 unrelated patients with XLP. This gene encodes a predicted protein of 128 amino acids composing a single SH2 domain with extensive homology to the SH2 domain of SHIP, an inositol polyphosphate 5-phosphatase that functions as a negative regulator of lymphocyte activation. DSHP is expressed in transformed T cell lines and is induced following in vitro activation of peripheral blood T lymphocytes. Expression of DSHP is restricted in vivo to lymphoid tissues, and RNA in situ hybridization demonstrates DSHP expression in activated T and B cell regions of reactive lymph nodes and in both T and B cell neoplasms. These observations confirm the identity of DSHP as the gene responsible for XLP, and suggest a role in the regulation of lymphocyte activation and proliferation. Induction of DSHP may sustain the immune response by interfering with SHIP-mediated inhibition of lymphocyte activation, while its inactivation in XLP patients results in a selective immunodeficiency to EBV.

An arsenate tolerance gene on chromosome 6 of rice

The genetics of arsenic tolerance in plants has not been extensively studied and no arsenic tolerance gene has been genetically mapped. Screening 20 diverse genotypes of rice for reduced root growth in 13.3 μM arsenate identified marked differences in tolerance. The most sensitive variety, Dawn, is known to be highly susceptible to straighthead, a condition linked to arsenic contamination of soil. Screening 108 recombinant inbred lines of the Bala x Azucena mapping population revealed the presence of a major gene, AsTol, which mapped between markers RZ516 and RG213 on chromosome 6. This gene is a good target for further characterisation. It should prove valuable for investigations into the physiological and molecular mechanism behind arsenic tolerance in plants and may lead to strategies aimed at breeding for arsenic contaminated regions. © New Phytologist (2004).

Autosomal dominant retinitis pigmentosa (ADRP):localization of an ADRP gene to the long arm of chromosome 3

Members of a large pedigree of Irish origin presenting with early onset Type I autosomal dominant retinitis pigmentosa (ADRP) have been typed for D3S47 (C17), a polymorphic marker from the long arm of chromosome 3. Significant, tight linkage of ADRP to D3S47, with a lod score of 14.7 maximizing at 0.00 recombination, has been obtained, hence localizing the ADRP gene (RP1) segregating in this pedigree to 3q.

Connectivity Mapping (ssCMap) To Predict A20 Inducing Drugs: Anti-inflammatory Action In Cystic Fibrosis

Cystic Fibrosis (CF) lung disease is characterised by a chronic and exaggerated inflammation in the airways. Despite recent developments to therapeutically overcome the underlying functional defect in CFTR (cystic fibrosis transmembrane conductance regulator), there is still an unmet need to also normalise the inflammatory response. The prolonged and heightened inflammatory response in CF is in part mediated by a lack of intrinsic downregulation of the pro-inflammatory NF-kB pathway. We have previously identified reduced expression of the NF-kB down-regulator A20 in CF as a key target to normalise the inflammatory response. Here we have used publically available gene array expression data together with sscMap (statistically significant connections’map)to successfully predict drugs already licensed for the use in humans to induce A20 mRNA and protein expression and thereby reduce inflammation. The effect of the predicted drugs on A20 and NFkB (p65) expression (mRNA) as well as pro-inflammatory cytokine release (IL-8) in the presence and absence of bacterial LPS was shown in bronchial epithelial cells lines (16HBE14o-, CFBE41o-) and in primary nasal epithelial cells (PNECs) from patients with CF (Phe508del homozygous) and non-CF controls. Additionally, the specificity of the drug action on A20 was confirmed using cell lines with TNFAIP3 (A20) knockdown (siRNA). We also show that the A20 inducing effect of ikarugamycin and quercetin is lower in CF derived airway epithelial cells than in non-CF cells.

Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma.

Multiple myeloma is characterized by genomic alterations frequently involving gains and losses of chromosomes. Single nucleotide polymorphism (SNP)-based mapping arrays allow the identification of copy number changes at the sub-megabase level and the identification of loss of heterozygosity (LOH) due to monosomy and uniparental disomy (UPD). We have found that SNP-based mapping array data and fluorescence in situ hybridization (FISH) copy number data correlated well, making the technique robust as a tool to investigate myeloma genomics. The most frequently identified alterations are located at 1p, 1q, 6q, 8p, 13, and 16q. LOH is found in these large regions and also in smaller regions throughout the genome with a median size of 1 Mb. We have identified that UPD is prevalent in myeloma and occurs through a number of mechanisms including mitotic nondisjunction and mitotic recombination. For the first time in myeloma, integration of mapping and expression data has allowed us to reduce the complexity of standard gene expression data and identify candidate genes important in both the transition from normal to monoclonal gammopathy of unknown significance (MGUS) to myeloma and in different subgroups within myeloma. We have documented these genes, providing a focus for further studies to identify and characterize those that are key in the pathogenesis of myeloma.