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.

SORBS1 gene, a new candidate for diabetic nephropathy: results from a multi-stage genome-wide association study in patients with type 1 diabetes

Aims/hypothesis

A Candidate Gene Association Study Identifies DAPL1 as a Female-Specific Susceptibility Locus for Age-Related Macular Degeneration (AMD)

Age-related macular degeneration (AMD) is the leading cause of blindness among white caucasians over the age of 50 years with a prevalence rate expected to increase markedly with an anticipated increase in the life span of the world population. To further expand our knowledge of the genetic architecture of the disease, we pursued a candidate gene approach assessing 25 genes and a total of 109 variants. Of these, synonymous single nucleotide polymorphism (SNP) rs17810398 located in death-associated protein-like 1 (DAPL1) was found to be associated with AMD in a joint analysis of 3,229 cases and 2,835 controls from five studies [combined P ADJ = 1.15 × 10(-6), OR 1.332 (1.187-1.496)]. This association was characterized by a highly significant sex difference (P diff = 0.0032) in that it was clearly confined to females with genome-wide significance [P ADJ = 2.62 × 10(-8), OR 1.541 (1.324-1.796); males: P ADJ = 0.382, OR 1.084 (0.905-1.298)]. By targeted resequencing of risk and non-risk associated haplotypes in the DAPL1 locus, we identified additional potentially functional risk variants, namely a common 897-bp deletion and a SNP predicted to affect a putative binding site of an exonic splicing enhancer. We show that the risk haplotype correlates with a reduced retinal transcript level of two, less frequent, non-canonical DAPL1 isoforms. DAPL1 plays a role in epithelial differentiation and may be involved in apoptotic processes thereby suggesting a possible novel pathway in AMSaveD pathogenesis.

An siRNA-based functional genomics screen for the identification of regulators of ciliogenesis and ciliopathy genes

Defects in primary cilium biogenesis underlie the ciliopathies, a growing group of genetic disorders. We describe a whole-genome siRNA-based reverse genetics screen for defects in biogenesis and/or maintenance of the primary cilium, obtaining a global resource. We identify 112 candidate ciliogenesis and ciliopathy genes, including 44 components of the ubiquitin-proteasome system, 12 G-protein-coupled receptors, and 3 pre-mRNA processing factors (PRPF6, PRPF8 and PRPF31) mutated in autosomal dominant retinitis pigmentosa. The PRPFs localize to the connecting cilium, and PRPF8- and PRPF31-mutated cells have ciliary defects. Combining the screen with exome sequencing data identified recessive mutations in PIBF1, also known as CEP90, and C21orf2, also known as LRRC76, as causes of the ciliopathies Joubert and Jeune syndromes. Biochemical approaches place C21orf2 within key ciliopathy-associated protein modules, offering an explanation for the skeletal and retinal involvement observed in individuals with C21orf2 variants. Our global, unbiased approaches provide insights into ciliogenesis complexity and identify roles for unanticipated pathways in human genetic disease.

Connectivity mapping using a combined gene signature from multiple colorectal cancer datasets identified candidate drugs including existing chemotherapies

BACKGROUND: While the discovery of new drugs is a complex, lengthy and costly process, identifying new uses for existing drugs is a cost-effective approach to therapeutic discovery. Connectivity mapping integrates gene expression profiling with advanced algorithms to connect genes, diseases and small molecule compounds and has been applied in a large number of studies to identify potential drugs, particularly to facilitate drug repurposing. Colorectal cancer (CRC) is a commonly diagnosed cancer with high mortality rates, presenting a worldwide health problem. With the advancement of high throughput omics technologies, a number of large scale gene expression profiling studies have been conducted on CRCs, providing multiple datasets in gene expression data repositories. In this work, we systematically apply gene expression connectivity mapping to multiple CRC datasets to identify candidate therapeutics to this disease.

RESULTS: We developed a robust method to compile a combined gene signature for colorectal cancer across multiple datasets. Connectivity mapping analysis with this signature of 148 genes identified 10 candidate compounds, including irinotecan and etoposide, which are chemotherapy drugs currently used to treat CRCs. These results indicate that we have discovered high quality connections between the CRC disease state and the candidate compounds, and that the gene signature we created may be used as a potential therapeutic target in treating the disease. The method we proposed is highly effective in generating quality gene signature through multiple datasets; the publication of the combined CRC gene signature and the list of candidate compounds from this work will benefit both cancer and systems biology research communities for further development and investigations.

Molecular subtyping of primary prostate cancer reveals specific and shared target genes of different ETS rearrangements

This work aimed to evaluate whether ETS transcription factors frequently involved in rearrangements in prostate carcinomas (PCa), namely ERG and ETV1, regulate specific or shared target genes. We performed differential expression analysis on nine normal prostate tissues and 50 PCa enriched for different ETS rearrangements using exon-level expression microarrays, followed by in vitro validation using cell line models. We found specific deregulation of 57 genes in ERG-positive PCa and 15 genes in ETV1-positive PCa, whereas deregulation of 27 genes was shared in both tumor subtypes. We further showed that the expression of seven tumor-associated ERG target genes (PLA1A, CACNA1D, ATP8A2, HLA-DMB, PDE3B, TDRD1, and TMBIM1) and two tumor-associated ETV1 target genes (FKBP10 and GLYATL2) was significantly affected by specific ETS silencing in VCaP and LNCaP cell line models, respectively, whereas the expression of three candidate ERG and ETV1 shared targets (GRPR, KCNH8, and TMEM45B) was significantly affected by silencing of either ETS. Interestingly, we demonstrate that the expression of TDRD1, the topmost overexpressed gene of our list of ERG-specific candidate targets, is inversely correlated with the methylation levels of a CpG island found at -66 bp of the transcription start site in PCa and that TDRD1 expression is regulated by direct binding of ERG to the CpG island in VCaP cells. We conclude that ETS transcription factors regulate specific and shared target genes and that TDRD1, FKBP10, and GRPR are promising therapeutic targets and can serve as diagnostic markers for molecular subtypes of PCa harboring specific fusion gene rearrangements.

Selecting Relevant Genes with a Spectral Approach

Array technologies have made it possible to record simultaneously the expression pattern of thousands of genes. A fundamental problem in the analysis of gene expression data is the identification of highly relevant genes that either discriminate between phenotypic labels or are important with respect to the cellular process studied in the experiment: for example cell cycle or heat shock in yeast experiments, chemical or genetic perturbations of mammalian cell lines, and genes involved in class discovery for human tumors. In this paper we focus on the task of unsupervised gene selection. The problem of selecting a small subset of genes is particularly challenging as the datasets involved are typically characterized by a very small sample size ?? the order of few tens of tissue samples ??d by a very large feature space as the number of genes tend to be in the high thousands. We propose a model independent approach which scores candidate gene selections using spectral properties of the candidate affinity matrix. The algorithm is very straightforward to implement yet contains a number of remarkable properties which guarantee consistent sparse selections. To illustrate the value of our approach we applied our algorithm on five different datasets. The first consists of time course data from four well studied Hematopoietic cell lines (HL-60, Jurkat, NB4, and U937). The other four datasets include three well studied treatment outcomes (large cell lymphoma, childhood medulloblastomas, breast tumors) and one unpublished dataset (lymph status). We compared our approach both with other unsupervised methods (SOM,PCA,GS) and with supervised methods (SNR,RMB,RFE). The results clearly show that our approach considerably outperforms all the other unsupervised approaches in our study, is competitive with supervised methods and in some case even outperforms supervised approaches.

"Genetic Diversity and Selection in Three Plasmodium vivax Merozoite Surface Protein 7 (Pvmsp-7) Genes in a Colombian Population"

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Genetic diversity and selection in three plasmodium vivax merozoite surface protein 7 (pvmsp-7) genes in a colombian population

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Generation of candidate human influenza vaccine strains in cell culture: rehearsing the European response to an H7N1 pandemic threat

Background: Although H5N1 avian influenza viruses pose the most obvious imminent pandemic threat, there have been several recent zoonotic incidents involving transmission of H7 viruses to humans. Vaccines are the primary public health defense against pandemics, but reliance on embryonated chickens eggs to propagate vaccine and logistic problems posed by the use of new technology may slow our ability to respond rapidly in a pandemic situation. Objectives: We sought to generate an H7 candidate vaccine virus suitable for administration to humans whose generation and amplification avoided the use of eggs. Methods: We generated a suitable H7 vaccine virus by reverse genetics. This virus, known as RD3, comprises the internal genes of A/Puerto Rico/8/34 with surface antigens of the highly pathogenic avian strain A/Chicken/Italy/13474/99 (H7N1). The multi-basic amino acid site in the HA gene, associated with high pathogenicity in chickens, was removed. Results: The HA modification did not alter the antigenicity of the virus and the resultant single basic motif was stably retained following several passages in Vero and PER. C6 cells. RD3 was attenuated for growth in embryonated eggs, chickens, and ferrets. RD3 induced an antibody response in infected animals reactive against both the homologous virus and other H7 influenza viruses associated with recent infection by H7 viruses in humans. Conclusions: This is the first report of a candidate H7 vaccine virus for use in humans generated by reverse genetics and propagated entirely in mammalian tissue culture. The vaccine has potential use against a wide range of H7 strains.

Expression of target and reference genes in Daphnia magna exposed to ibuprofen

Background: Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e. g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application. Results: We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l(-1). Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was upregulated following 24 h exposure to 63-81 mg IB l(-1). Conclusions: As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.

Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

Identification and expression analyses of Cytosolic Glutamine Synthetase genes in barley (Hordeum vulgare L.)

Glutamine synthetase (GS) is a key enzyme in nitrogen (N) assimilation, particularly during seed development. Three cytosolic GS isoforms (HvGS1) were identified in barley (Hordeum vulgare L. cv Golden Promise). Quantitation of gene expression, localization and response to N supply revealed that each gene plays a non-redundant role in different tissues and during development. Localization of HvGS1_1 in vascular cells of different tissues, combined with its abundance in the stem and its response to changes in N supply, indicate that it is important in N transport and remobilization. HvGS1_1 is located on chromosome 6H at 72.54 cM, close to the marker HVM074 which is associated with a major quantitative trait locus (QTL) for grain protein content (GPC). HvGS1_1 may be a potential candidate gene to manipulate barley GPC. HvGS1_2 mRNA was localized to the leaf mesophyll cells, in the cortex and pericycle of roots, and was the dominant HvGS1 isoform in these tissues. HvGS1_2 expression increased in leaves with an increasing supply of N, suggesting its role in the primary assimilation of N. HvGS1_3 was specifically and predominantly localized in the grain, being highly expressed throughout grain development. HvGS1_3 expression increased specifically in the roots of plants grown on high NH+4, suggesting that it has a primary role in grain N assimilation and also in the protection against ammonium toxicity in roots. The expression of HvGS1 genes is directly correlated with protein and enzymatic activity, indicating that transcriptional regulation is of prime importance in the control of GS activity in barley.

Arabinose and protocatechuate catabolism genes are important for growth of Rhizobium leguminosarum biovar viciae in the pea rhizosphere

Background and aims: To form nitrogen-fixing nodules on pea roots, Rhizobium leguminosarum biovar viciae must be competitive in the rhizosphere. Our aim was to identify genes important for rhizosphere fitness. Methods: Signature-tagged mutants were screened using microarrays to identify mutants reduced for growth in pea rhizospheres. Candidate mutants were assessed relative to controls for growth in minimal medium, growth in pea rhizospheres and for infection of peas in mixed inoculants. Mutated genes were identified by DNA sequencing and confirmed by transduction. Results: Of 5508 signature-tagged mutants, microarrays implicated 50 as having decreased rhizosphere fitness. Growth tests identified six mutants with rhizosphere-specific phenotypes. The mutation in one of the genes (araE) was in an arabinose catabolism operon and blocked growth on arabinose. The mutation in another gene (pcaM), encoding a predicted solute binding protein for protocatechuate and hydroxybenzoate uptake, decreased growth on protocatechuate. Both mutants were decreased for nodule infection competitiveness with mixed inoculants, but nodulated peas normally when inoculated alone. Other mutants with similar phenotypes had mutations predicted to affect secondary metabolism. Conclusions: Catabolism of arabinose and protocatechuate in the pea rhizosphere is important for competitiveness of R.l. viciae. Other genes predicted to be involved in secondary metabolism are also important.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7 Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5` and 3` segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition. we uncovered two surprising patterns of evolution in repetitive sequence. Firsts in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi. (C) 2009 Elsevier B.V. All rights reserved.