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.

Genome wide interactions of wild-type and activator bypass forms of σ54

Enhancer-dependent transcription involving the promoter specificity factor σ54 is widely distributed amongst bacteria and commonly associated with cell envelope function. For transcription initiation, σ54-RNA polymerase yields open promoter complexes through its remodelling by cognate AAA+ ATPase activators. Since activators can be bypassed in vitro, bypass transcription in vivo could be a source of emergent gene expression along evolutionary pathways yielding new control networks and transcription patterns. At a single test promoter in vivo bypass transcription was not observed. We now use genome-wide transcription profiling, genome-wide mutagenesis and gene over-expression strategies in Escherichia coli, to (i) scope the range of bypass transcription in vivo and (ii) identify genes which might alter bypass transcription in vivo. We find little evidence for pervasive bypass transcription in vivo with only a small subset of σ54 promoters functioning without activators. Results also suggest no one gene limits bypass transcription in vivo, arguing bypass transcription is strongly kept in check. Promoter sequences subject to repression by σ54 were evident, indicating loss of rpoN (encoding σ54) rather than creating rpoN bypass alleles would be one evolutionary route for new gene expression patterns. Finally, cold-shock promoters showed unusual σ54-dependence in vivo not readily correlated with conventional σ54 binding-sites.

Genome-wide association study for sight-threatening diabetic retinopathy reveals association with genetic variation near the GRB2 gene

AIMS/HYPOTHESIS: Diabetic retinopathy is a serious complication of diabetes mellitus and can lead to blindness. A genetic component, in addition to traditional risk factors, has been well described although strong genetic factors have not yet been identified. Here, we aimed to identify novel genetic risk factors for sight-threatening diabetic retinopathy using a genome-wide association study.

METHODS: Retinopathy was assessed in white Australians with type 2 diabetes mellitus. Genome-wide association analysis was conducted for comparison of cases of sight-threatening diabetic retinopathy (n = 336) with diabetic controls with no retinopathy (n = 508). Top ranking single nucleotide polymorphisms were typed in a type 2 diabetes replication cohort, a type 1 diabetes cohort and an Indian type 2 cohort. A mouse model of proliferative retinopathy was used to assess differential expression of the nearby candidate gene GRB2 by immunohistochemistry and quantitative western blot.

RESULTS: The top ranked variant was rs3805931 with p = 2.66 × 10(-7), but no association was found in the replication cohort. Only rs9896052 (p = 6.55 × 10(-5)) was associated with sight-threatening diabetic retinopathy in both the type 2 (p = 0.035) and the type 1 (p = 0.041) replication cohorts, as well as in the Indian cohort (p = 0.016). The study-wide meta-analysis reached genome-wide significance (p = 4.15 × 10(-8)). The GRB2 gene is located downstream of this variant and a mouse model of retinopathy showed increased GRB2 expression in the retina.

CONCLUSIONS/INTERPRETATION: Genetic variation near GRB2 on chromosome 17q25.1 is associated with sight-threatening diabetic retinopathy. Several genes in this region are promising candidates and in particular GRB2 is upregulated during retinal stress and neovascularisation.

LD Score regression distinguishes confounding from polygenicity in genome-wide association studies

Both polygenicity (many small genetic effects) and confounding biases, such as cryptic relatedness and population stratification, can yield an inflated distribution of test statistics in genome-wide association studies (GWAS). However, current methods cannot distinguish between inflation from a true polygenic signal and bias. We have developed an approach, LD Score regression, that quantifies the contribution of each by examining the relationship between test statistics and linkage disequilibrium (LD). The LD Score regression intercept can be used to estimate a more powerful and accurate correction factor than genomic control. We find strong evidence that polygenicity accounts for the majority of the inflation in test statistics in many GWAS of large sample size.

Genome-wide association studies in oesophageal adenocarcinoma and Barrett's oesophagus: a large-scale meta-analysis

Background: Esophageal adenocarcinoma (EA) is one of the fastest rising cancers in western countries. Barrett’s Esophagus (BE) is the premalignant precursor of EA. However, only a subset of BE patients develop EA, which complicates the clinical management in the absence of valid predictors. Genetic risk factors for BE and EA are incompletely understood. This study aimed to identify novel genetic risk factors for BE and EA.Methods: Within an international consortium of groups involved in the genetics of BE/EA, we performed the first meta-analysis of all genome-wide association studies (GWAS) available, involving 6,167 BE patients, 4,112 EA patients, and 17,159 representative controls, all of European ancestry, genotyped on Illumina high-density SNP-arrays, collected from four separate studies within North America, Europe, and Australia. Meta-analysis was conducted using the fixed-effects inverse variance-weighting approach. We used the standard genome-wide significant threshold of 5×10-8 for this study. We also conducted an association analysis following reweighting of loci using an approach that investigates annotation enrichment among the genome-wide significant loci. The entire GWAS-data set was also analyzed using bioinformatics approaches including functional annotation databases as well as gene-based and pathway-based methods in order to identify pathophysiologically relevant cellular pathways.Findings: We identified eight new associated risk loci for BE and EA, within or near the CFTR (rs17451754, P=4·8×10-10), MSRA (rs17749155, P=5·2×10-10), BLK (rs10108511, P=2·1×10-9), KHDRBS2 (rs62423175, P=3·0×10-9), TPPP/CEP72 (rs9918259, P=3·2×10-9), TMOD1 (rs7852462, P=1·5×10-8), SATB2 (rs139606545, P=2·0×10-8), and HTR3C/ABCC5 genes (rs9823696, P=1·6×10-8). A further novel risk locus at LPA (rs12207195, posteriori probability=0·925) was identified after re-weighting using significantly enriched annotations. This study thereby doubled the number of known risk loci. The strongest disease pathways identified (P<10-6) belong to muscle cell differentiation and to mesenchyme development/differentiation, which fit with current pathophysiological BE/EA concepts. To our knowledge, this study identified for the first time an EA-specific association (rs9823696, P=1·6×10-8) near HTR3C/ABCC5 which is independent of BE development (P=0·45).Interpretation: The identified disease loci and pathways reveal new insights into the etiology of BE and EA. Furthermore, the EA-specific association at HTR3C/ABCC5 may constitute a novel genetic marker for the prediction of transition from BE to EA. Mutations in CFTR, one of the new risk loci identified in this study, cause cystic fibrosis (CF), the most common recessive disorder in Europeans. Gastroesophageal reflux (GER) belongs to the phenotypic CF-spectrum and represents the main risk factor for BE/EA. Thus, the CFTR locus may trigger a common GER-mediated pathophysiology.

Genome-wide association study reveals greater polygenic loading for schizophrenia in cases with a family history of illness

Genome-wide association studies (GWAS) of schizophrenia have yielded more than 100 common susceptibility variants, and strongly support a substantial polygenic contribution of a large number of small allelic effects. It has been hypothesized that familial schizophrenia is largely a consequence of inherited rather than environmental factors. We investigated the extent to which familiality of schizophrenia is associated with enrichment for common risk variants detectable in a large GWAS. We analyzed single nucleotide polymorphism (SNP) data for cases reporting a family history of psychotic illness (N = 978), cases reporting no such family history (N = 4,503), and unscreened controls (N = 8,285) from the Psychiatric Genomics Consortium (PGC1) study of schizophrenia. We used a multinomial logistic regression approach with model-fitting to detect allelic effects specific to either family history subgroup. We also considered a polygenic model, in which we tested whether family history positive subjects carried more schizophrenia risk alleles than family history negative subjects, on average. Several individual SNPs attained suggestive but not genome-wide significant association with either family history subgroup. Comparison of genome-wide polygenic risk scores based on GWAS summary statistics indicated a significant enrichment for SNP effects among family history positive compared to family history negative cases (Nagelkerke's R(2 ) = 0.0021; P = 0.00331; P-value threshold <0.4). Estimates of variability in disease liability attributable to the aggregate effect of genome-wide SNPs were significantly greater for family history positive compared to family history negative cases (0.32 and 0.22, respectively; P = 0.031). We found suggestive evidence of allelic effects detectable in large GWAS of schizophrenia that might be specific to particular family history subgroups. However, consideration of a polygenic risk score indicated a significant enrichment among family history positive cases for common allelic effects. Familial illness might, therefore, represent a more heritable form of schizophrenia, as suggested by previous epidemiological studies.

Design, recruitment, logistics, and data management of the GEHA (Genetics of Healthy Ageing) project

In 2004, the integrated European project GEHA (Genetics of Healthy Ageing) was initiated with the aim of identifying genes involved in healthy ageing and longevity. The first step in the project was the recruitment of more than 2500 pairs of siblings aged 90 years or more together with one younger control person from 15 areas in 11 European countries through a coordinated and standardised effort. A biological sample, preferably a blood sample, was collected from each participant, and basic physical and cognitive measures were obtained together with information about health, life style, and family composition. From 2004 to 2008 a total of 2535 families comprising 5319 nonagenarian siblings were identified and included in the project. In addition, 2548 younger control persons aged 50-75 years were recruited. A total of 2249 complete trios with blood samples from at least two old siblings and the younger control were formed and are available for genetic analyses (e.g. linkage studies and genome-wide association studies). Mortality follow-up improves the possibility of identifying families with the most extreme longevity phenotypes. With a mean follow-up time of 3.7 years the number of families with all participating siblings aged 95 years or more has increased by a factor of 5 to 750 families compared to when interviews were conducted. Thus, the GEHA project represents a unique source in the search for genes related to healthy ageing and longevity.

Association Study of 167 Candidate Genes for Schizophrenia Selected by a Multi-Domain Evidence-Based Prioritization Algorithm and Neurodevelopmental Hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this
ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

Background: Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

Results: First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

Conclusions: The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.

Genetics of Healthy Aging in Europe.The EU-Integrated Project GEHA (GEnetics of Healthy Aging)

The aim of the 5-year European Union (EU)-Integrated Project GEnetics of Healthy Aging (GEHA), constituted by 25 partners (24 from Europe plus the Beijing Genomics Institute from China), is to identify genes involved in healthy aging and longevity, which allow individuals to survive to advanced old age in good cognitive and physical function and in the absence of major age-related diseases. To achieve this aim a coherent, tightly integrated program of research that unites demographers, geriatricians, geneticists, genetic epidemiologists, molecular biologists, bioinfomaticians, and statisticians has been set up. The working plan is to: (a) collect DNA and information on the health status from an unprecedented number of long-lived 90+ sibpairs (n = 2650) and of younger ethnically matched controls (n = 2650) from 11 European countries; (b) perform a genome-wide linkage scannning in all the sibpairs (a total of 5300 individuals); this investigation will be followed by linkage disequilibrium mapping (LD mapping) of the candidate chromosomal regions; (c) study in cases (i.e., the 2650 probands of the sibpairs) and controls (2650 younger people), genomic regions (chromosome 4, D4S1564, chromosome 11, 11.p15.5) which were identified in previous studies as possible candidates to harbor longevity genes; (d) genotype all recruited subjects for apoE polymorphisms; and (e) genotype all recruited subjects for inherited as well as epigenetic variability of the mitochondrial DNA (mtDNA). The genetic analysis will be performed by 9 high-throughput platforms, within the framework of centralized databases for phenotypic, genetic, and mtDNA data. Additional advanced approaches (bioinformatics, advanced statistics, mathematical modeling, functional genomics and proteomics, molecular biology, molecular genetics) are envisaged to identify the gene variant(s) of interest. The experimental design will also allow (a) to identify gender-specific genes involved in healthy aging and longevity in women and men stratified for ethnic and geographic origin and apoE genotype; (b) to perform a longitudinal survival study to assess the impact of the identified genetic loci on 90+ people mortality; and (c) to develop mathematical and statistical models capable of combining genetic data with demographic characteristics, health status, socioeconomic factors, lifestyle habits.

Genetic analysis of coronary artery disease single nucleotide polymorphisms in diabetic nephropathy

Background. Diabetic nephropathy is a leading cause of end-stage renal disease. Premature mortality is common in patients with nephropathy, largely due to cardiovascular disease. Genetic variants implicated in macrovascular disease are therefore excellent candidates to assess for association with diabetic nephropathy. Recent genome-wide association studies have identified a total of 15 single-nucleotide polymorphisms (SNPs) that are reproducibly associated with cardiovascular disease.

Methods. We initially assessed these SNPs for association in UK type 1 diabetic patients with (cases; n = 597) and without (controls; n = 502) nephropathy using iPLEXTM and TaqMan® assays. Replication studies were performed with DNA genotyped in a total of 2668 individuals from the British Isles.

Results. One SNP (rs4420638) on chromosome 19q13 was found to be significantly associated with diabetic nephropathy before (P = 0.0002) and after correction for multiple testing (Pcorrected = 0.002). We replicated this finding in a phenotypically similar case–control collection comprising 709 individuals with type 1 diabetes (P = 0.002; combined P < 0.00001; OR = 1.54, 95% CI: 1.29–1.84).

Conclusions. Our case–control data suggest that rs4420638, or a functional SNP in linkage disequilibrium with this SNP, may be associated with diabetic nephropathy.

The impact of newly identified loci on coronary heart disease, stroke and total mortality in the MORGAM prospective cohorts

Recently, genome wide association studies (GWAS) have identified a number of single nucleotide polymorphisms (SNPs) as being associated with coronary heart disease (CHD). We estimated the effect of these SNPs on incident CHD, stroke and total mortality in the prospective cohorts of the MORGAM Project. We studied cohorts from Finland, Sweden, France and Northern Ireland (total N=33,282, including 1,436 incident CHD events and 571 incident stroke events). The lead SNPs at seven loci identified thus far and additional SNPs (in total 42) were genotyped using a case-cohort design.We estimated the effect of the SNPs on disease history at baseline, disease events during follow-up and classic risk factors. Multiple testing was taken into account using false discovery rate (FDR) analysis. SNP rs1333049 on chromosome 9p21.3 was associated with both CHD and stroke (HR5=.20, 95% CI 1.08-1.34 for incident CHD events and 1.15, 0.99-1.34 for incident stroke). SNP rs11670734 (19q12) was associated with total mortality and stroke. SNP rs2146807 (10q11.21) showed some association with the fatality of acute coronary event. SNP rs2943634 (2q36.3) was associated with high density lipoprotein (HDL) cholesterol and SNPs rs599839, rs4970834 (1p13.3) and rs17228212 (15q22.23) were associated with non-HDL cholesterol. SNPs rs2943634 (2q36.3) and rs12525353 (6q25.1) were associated with blood pressure. These findings underline the need for replication studies in prospective settings and confirm the candidacy of several SNPs that may play a role in the etiology of cardiovascular disease.

Molecular analysis of endothelial progenitor cell (EPC) subtypes reveals two distinct cell populations with different identities

BACKGROUND: The term endothelial progenitor cells (EPCs) is currently used to refer to cell populations which are quite dissimilar in terms of biological properties. This study provides a detailed molecular fingerprint for two EPC subtypes: early EPCs (eEPCs) and outgrowth endothelial cells (OECs). METHODS: Human blood-derived eEPCs and OECs were characterised by using genome-wide transcriptional profiling, 2D protein electrophoresis, and electron microscopy. Comparative analysis at the transcript and protein level included monocytes and mature endothelial cells as reference cell types. RESULTS: Our data show that eEPCs and OECs have strikingly different gene expression signatures. Many highly expressed transcripts in eEPCs are haematopoietic specific (RUNX1, WAS, LYN) with links to immunity and inflammation (TLRs, CD14, HLAs), whereas many transcripts involved in vascular development and angiogenesis-related signalling pathways (Tie2, eNOS, Ephrins) are highly expressed in OECs. Comparative analysis with monocytes and mature endothelial cells clusters eEPCs with monocytes, while OECs segment with endothelial cells. Similarly, proteomic analysis revealed that 90% of spots identified by 2-D gel analysis are common between OECs and endothelial cells while eEPCs share 77% with monocytes. In line with the expression pattern of caveolins and cadherins identified by microarray analysis, ultrastructural evaluation highlighted the presence of caveolae and adherens junctions only in OECs. CONCLUSIONS: This study provides evidence that eEPCs are haematopoietic cells with a molecular phenotype linked to monocytes; whereas OECs exhibit commitment to the endothelial lineage. These findings indicate that OECs might be an attractive cell candidate for inducing therapeutic angiogenesis, while eEPC should be used with caution because of their monocytic nature.