Gene Sets Net Correlations Analysis (GSNCA): a multivariate differential coexpression test for gene sets

Motivation: To date, Gene Set Analysis (GSA) approaches primarily focus on identifying differentially expressed gene sets (pathways). Methods for identifying differentially coexpressed pathways also exist but are mostly based on aggregated pairwise correlations, or other pairwise measures of coexpression. Instead, we propose Gene Sets Net Correlations Analysis (GSNCA), a multivariate differential coexpression test that accounts for the complete correlation structure between genes.

Results: In GSNCA, weight factors are assigned to genes in proportion to the genes' cross-correlations (intergene correlations). The problem of finding the weight vectors is formulated as an eigenvector problem with a unique solution. GSNCA tests the null hypothesis that for a gene set there is no difference in the weight vectors of the genes between two conditions. In simulation studies and the analyses of experimental data, we demonstrate that GSNCA, indeed, captures changes in the structure of genes' cross-correlations rather than differences in the averaged pairwise correlations. Thus, GSNCA infers differences in coexpression networks, however, bypassing method-dependent steps of network inference. As an additional result from GSNCA, we define hub genes as genes with the largest weights and show that these genes correspond frequently to major and specific pathway regulators, as well as to genes that are most affected by the biological difference between two conditions. In summary, GSNCA is a new approach for the analysis of differentially coexpressed pathways that also evaluates the importance of the genes in the pathways, thus providing unique information that may result in the generation of novel biological hypotheses.

Rare and common variants in extracellular matrix gene Fibrillin 2 (FBN2) are associated with macular degeneration

Neurodegenerative diseases affecting the macula constitute a major cause of incurable vision loss and exhibit considerable clinical and genetic heterogeneity, from early-onset monogenic disease to multifactorial late-onset age-related macular degeneration (AMD). As part of our continued efforts to define genetic causes of macular degeneration, we performed whole exome sequencing in four individuals of a two-generation family with autosomal dominant maculopathy and identified a rare variant p.Glu1144Lys in Fibrillin 2 (FBN2), a glycoprotein of the elastin-rich extracellular matrix (ECM). Sanger sequencing validated the segregation of this variant in the complete pedigree, including two additional affected and one unaffected individual. Sequencing of 192 maculopathy patients revealed additional rare variants, predicted to disrupt FBN2 function. We then undertook additional studies to explore the relationship of FBN2 to macular disease. We show that FBN2 localizes to Bruch's membrane and its expression appears to be reduced in aging and AMD eyes, prompting us to examine its relationship with AMD. We detect suggestive association of a common FBN2 non-synonymous variant, rs154001 (p.Val965Ile) with AMD in 10,337 cases and 11,174 controls (OR=1.10; p-value=3.79×10(-5)). Thus, it appears that rare and common variants in a single gene - FBN2 - can contribute to Mendelian and complex forms of macular degeneration. Our studies provide genetic evidence for a key role of elastin microfibers and Bruch's membrane in maintaining blood-retina homeostasis and establish the importance of studying orphan diseases for understanding more common clinical phenotypes.

Mutational Analysis of the Rhodopsin Gene in Sector Retinitis Pigmentosa

Background: To determine the role of rhodopsin (RHO) gene mutations in patients with sector retinitis pigmentosa (RP) from Northern Ireland.

Design: A case series of sector RP in a tertiary ocular genetics clinic.

Participants: Four patients with sector RP were recruited from the Royal Victoria Hospital (Belfast, Northern Ireland) and Altnagelvin Hospital (Londonderry, Northern Ireland) following informed consent.

Methods: The diagnosis of sector RP was based on clinical examination, International Society for Clinical Electrophysiology of Vision (ISCEV) standard electrophysiology, and visual field analysis. DNA was extracted from peripheral blood leucocytes and the coding regions and adjacent flanking intronic sequences of the RHO gene were polymerase chain reaction (PCR) amplified and cycle sequenced.

Main Outcome Measure: Rhodopsin mutational status.

Results: A heterozygous missense mutation in RHO (c.173C > T) resulting in a non-conservative substitution of threonine to methionine (p. Thr58Met) was identified in one patient and was absent from 360 control individuals. This non-conservative substitution (p.Thr58Met) replaces a highly evolutionary conserved polar hydrophilic threonine residue with a non-polar hydrophobic methionine residue at position 58 near the cytoplasmic border of helix A of RHO.

Conclusions: The study identified a RHO gene mutation (p.Thr58Met) not previously reported in RP in a patient with sector RP. These findings outline the phenotypic variability associated with RHO mutations. It has been proposed that the regional effects of RHO mutations are likely to result from interplay between mutant alleles and other genetic, epigenetic and environmental factors.

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

Gene expression in normal Urothelium depends on location within the bladder:A possible link to bladder carcinogenesis

Objectives: Clinical studies have shown that more than 70% of primary bladder tumours arise in the area around the ureteric orifice. In this study a genomic approach was taken to explore the molecular mechanisms that may influence this phenomenon.

Methods: RNA was isolated from each individual normal ureteric orifice and the dome biopsy from 33 male patients. Equal amounts of the pooled ureteric orifice and dome mRNAs were labelled with Cy3 and Cy5, respectively before hybridising to the gene chip (UniGEM 2.0, Incyte Genomics Inc., Wilmington, Delaware, USA). Results: Significant changes (more than a twofold difference) in gene expression were observed in 3.1% (312) of the 10,176 gene array: 211 genes upregulated and 101 downregulated. Analysis of Cdc25B, TK1, PKM, and PDGFra with RT-PCR supported the reliability of the microarray result. Seladin-1 was the most upregulated gene in the ureteric orifice: 8.3-fold on the microarray and 11.4-fold by real time PCR.

Conclusions: Overall, this study suggests significant altered gene expression between these two anatomically distinct areas of the normal human bladder. Of particular note is Seladin-1, whose significance in cancer is yet to be clarified. Further studies of the genes discovered by this work will help clarify which of these differences influence primary bladder carcinogenesis. (c) 2006 European Association of Urology. Published by Elsevier B.V. All rights reserved.

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.

In Helicobacter pylori auto-inducer-2, but not LuxS/MccAB catalysed reverse transsulphuration, regulates motility through modulation of flagellar gene transcription

BACKGROUND: LuxS may function as a metabolic enzyme or as the synthase of a quorum sensing signalling molecule, auto-inducer-2 (AI-2); hence, the mechanism underlying phenotypic changes upon luxS inactivation is not always clear. In Helicobacter pylori, we have recently shown that, rather than functioning in recycling methionine as in most bacteria, LuxS (along with newly-characterised MccA and MccB), synthesises cysteine via reverse transsulphuration. In this study, we investigated whether and how LuxS controls motility of H. pylori, specifically if it has its effects via luxS-required cysteine metabolism or via AI-2 synthesis only.

RESULTS: We report that disruption of luxS renders H. pylori non-motile in soft agar and by microscopy, whereas disruption of mccAHp or mccBHp (other genes in the cysteine provision pathway) does not, implying that the lost phenotype is not due to disrupted cysteine provision. The motility defect of the DeltaluxSHp mutant was complemented genetically by luxSHp and also by addition of in vitro synthesised AI-2 or 4, 5-dihydroxy-2, 3-pentanedione (DPD, the precursor of AI-2). In contrast, exogenously added cysteine could not restore motility to the DeltaluxSHp mutant, confirming that AI-2 synthesis, but not the metabolic effect of LuxS was important. Microscopy showed reduced number and length of flagella in the DeltaluxSHp mutant. Immunoblotting identified decreased levels of FlaA and FlgE but not FlaB in the DeltaluxSHp mutant, and RT-PCR showed that the expression of flaA, flgE, motA, motB, flhA and fliI but not flaB was reduced. Addition of DPD but not cysteine to the DeltaluxSHp mutant restored flagellar gene transcription, and the number and length of flagella.

CONCLUSIONS: Our data show that as well as being a metabolic enzyme, H. pylori LuxS has an alternative role in regulation of motility by modulating flagellar transcripts and flagellar biosynthesis through production of the signalling molecule AI-2.

The Fasciola hepatica genome: gene duplication and polymorphism reveals adaptation to the host environment and the capacity for rapid evolution

BACKGROUND: The liver fluke Fasciola hepatica is a major pathogen of livestock worldwide, causing huge economic losses to agriculture, as well as 2.4 million human infections annually.

RESULTS: Here we provide a draft genome for F. hepatica, which we find to be among the largest known pathogen genomes at 1.3 Gb. This size cannot be explained by genome duplication or expansion of a single repeat element, and remains a paradox given the burden it may impose on egg production necessary to transmit infection. Despite the potential for inbreeding by facultative self-fertilisation, substantial levels of polymorphism were found, which highlights the evolutionary potential for rapid adaptation to changes in host availability, climate change or to drug or vaccine interventions. Non-synonymous polymorphisms were elevated in genes shared with parasitic taxa, which may be particularly relevant for the ability of the parasite to adapt to a broad range of definitive mammalian and intermediate molluscan hosts. Large-scale transcriptional changes, particularly within expanded protease and tubulin families, were found as the parasite migrated from the gut, across the peritoneum and through the liver to mature in the bile ducts. We identify novel members of anti-oxidant and detoxification pathways and defined their differential expression through infection, which may explain the stage-specific efficacy of different anthelmintic drugs.

CONCLUSIONS: The genome analysis described here provides new insights into the evolution of this important pathogen, its adaptation to the host environment and external selection pressures. This analysis also provides a platform for research into novel drugs and vaccines.

Genome-wide gene pathway analysis of psychotic illness symptom dimensions based on a new schizophrenia-specific model of the OPCRIT

Empirically derived phenotypic measurements have the potential to enhance gene-finding efforts in schizophrenia. Previous research based on factor analyses of symptoms has typically included schizoaffective cases. Deriving factor loadings from analysis of only narrowly defined schizophrenia cases could yield more sensitive factor scores for gene pathway and gene ontology analyses. Using an Irish family sample, this study 1) factor analyzed clinician-rated Operational Criteria Checklist items in cases with schizophrenia only, 2) scored the full sample based on these factor loadings, and 3) implemented genome-wide association, gene-based, and gene-pathway analysis of these SCZ-based symptom factors (final N= 507). Three factors emerged from the analysis of the schizophrenia cases: a manic, a depressive, and a positive symptom factor. In gene-based analyses of these factors, multiple genes had q<. 0.01. Of particular interest are findings for PTPRG and WBP1L, both of which were previously implicated by the Psychiatric Genomics Consortium study of SCZ; results from this study suggest that variants in these genes might also act as modifiers of SCZ symptoms. Gene pathway analyses of the first factor indicated over-representation of glutamatergic transmission, GABA-A receptor, and cyclic GMP pathways. Results suggest that these pathways may have differential influence on affective symptom presentation in schizophrenia.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

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).

Gene expression and epigenetic discovery screen reveal methylation of SFRP2 in prostate cancer

Aberrant activation of Wnts is common in human cancers, including prostate. Hypermethylation associated transcriptional silencing of Wnt antagonist genes SFRPs (Secreted Frizzled-Related Proteins) is a frequent oncogenic event. The significance of this is not known in prostate cancer. The objectives of our study were to (i) profile Wnt signaling related gene expression and (ii) investigate methylation of Wnt antagonist genes in prostate cancer. Using TaqMan Low Density Arrays, we identified 15 Wnt signaling related genes with significantly altered expression in prostate cancer; the majority of which were upregulated in tumors. Notably, histologically benign tissue from men with prostate cancer appeared more similar to tumor (r = 0.76) than to benign prostatic hyperplasia (BPH; r = 0.57, p < 0.001). Overall, the expression profile was highly similar between tumors of high (≥ 7) and low (≤ 6) Gleason scores. Pharmacological demethylation of PC-3 cells with 5-Aza-CdR reactivated 39 genes (≥ 2-fold); 40% of which inhibit Wnt signaling. Methylation frequencies in prostate cancer were 10% (2/20) (SFRP1), 64.86% (48/74) (SFRP2), 0% (0/20) (SFRP4) and 60% (12/20) (SFRP5). SFRP2 methylation was detected at significantly lower frequencies in high-grade prostatic intraepithelial neoplasia (HGPIN; 30%, (6/20), p = 0.0096), tumor adjacent benign areas (8.82%, (7/69), p < 0.0001) and BPH (11.43% (4/35), p < 0.0001). The quantitative level of SFRP2 methylation (normalized index of methylation) was also significantly higher in tumors (116) than in the other samples (HGPIN = 7.45, HB = 0.47, and BPH = 0.12). We show that SFRP2 hypermethylation is a common event in prostate cancer. SFRP2 methylation in combination with other epigenetic markers may be a useful biomarker of prostate cancer.

IGFBP7 promoter methylation and gene expression analysis in prostate cancer

PURPOSE: IGFBP7 belongs to a family of insulin-like growth factor-1 regulatory binding proteins. IGFBP7 hypermethylation is associated with its down-regulation in various carcinomas. In prostate cancer IGFBP7 down-regulation has been widely reported but to our knowledge the mechanisms behind this event are unknown. We performed a denaturing high performance liquid chromatography screening and validation strategy to profile the methylation status of IGFBP7 in prostate cancer.

MATERIALS AND METHODS: We combined denaturing high performance liquid chromatography and bisulfite sequencing to examine IGFBP7 methylation in a panel of prostate cancer cell lines. Quantitative methylation specific polymerase chain reaction was used to determine methylation levels in prostate tissue specimens of primary prostate cancer, histologically benign prostate adjacent to tumor, high grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. IGFBP7 gene expression was measured by quantitative methylation specific polymerase chain reaction in cell lines and tissue specimens.

RESULTS: IGFBP7 was methylated in the 4 prostate cancer cell lines DU145, LNCaP, PC-3 and 22RV1. Quantitative methylation specific polymerase chain reaction analysis revealed that promoter methylation was associated with decreased IGFBP7 expression. Quantitative methylation specific polymerase chain reaction showed that IGFBP7 methylation was more frequently detected in prostate cancer (60% (31/52)) and high grade prostatic intraepithelial neoplasia (40% (6/15)) samples compared to histologically benign prostate adjacent to tumor (10%) and benign prostatic hyperplasia (0%) samples.

CONCLUSIONS: To our knowledge this is the first report of aberrant IGFBP7 promoter hypermethylation and concurrent IGFBP7 gene silencing in prostate cancer cell lines. Results demonstrate that CpG methylation of IGFBP7 may represent a novel biomarker of prostate cancer and pre-invasive neoplasms. Thus, future examination of IGFBP7 methylation and expression in a larger patient cohort, including bodily fluids, is justified to further evaluate its role in a diagnostic and prognostic setting.

Hypoxia response element-driven cytosine deaminase/5-fluorocytosine gene therapy system:a highly effective approach to overcome the dynamics of tumour hypoxia and enhance the radiosensitivity of prostate cancer cells in vitro

BACKGROUND: We proposed to exploit hypoxia-inducible factor (HIF)-1alpha overexpression in prostate tumours and use this transcriptional machinery to control the expression of the suicide gene cytosine deaminase (CD) through binding of HIF-1alpha to arrangements of hypoxia response elements. CD is a prodrug activation enzyme, which converts inactive 5-fluorocytosine to active 5-fluorouracil (5-FU), allowing selective killing of vector containing cells.

METHODS: We developed a pair of vectors, containing either five or eight copies of the hypoxia response element (HRE) isolated from the vascular endothelial growth factor (pH5VCD) or glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (pH8GCD) gene, respectively. The kinetics of the hypoxic induction of the vectors and sensitization effects were evaluated in 22Rv1 and DU145 cells in vitro.

RESULTS: The CD protein as selectively detected in lysates of transiently transfected 22Rv1 and DU145 cells following hypoxic exposure. This is the first evidence of GAPDH HREs being used to control a suicide gene therapy strategy. Detectable CD levels were sustained upon reoxygenation and prolonged hypoxic exposures. Hypoxia-induced chemoresistance to 5-FU was overcome in both cell lines treated with this suicide gene therapy approach. Hypoxic transfectants were sensitized to prodrug concentrations that were ten-fold lower than those that are clinically relevant. Moreover, the surviving fraction of reoxygenated transfectants could be further reduced with the concomitant delivery of clinically relevant single radiation doses.

CONCLUSIONS: This strategy thus has the potential to sensitize the hypoxic compartment of prostate tumours and improve the outcome of current therapies.