SNPExpress: integrated visualization of genome-wide genotypes, copy numbers and gene expression levels

Candidate pathway-based genome-wide association studies identify novel associations of genomic variants in the complement system associated with coronary artery disease

BACKGROUND

Genomic variants identified by genome-wide association studies (GWAS) explain <20% of heritability of coronary artery disease (CAD), thus many risk variants remain missing for CAD. Identification of new variants may unravel new biological pathways and genetic mechanisms for CAD. To identify new variants associated with CAD, we developed a candidate pathway-based GWAS by integrating expression quantitative loci analysis and mining of GWAS data with variants in a candidate pathway.

METHODS AND RESULTS

Mining of GWAS data was performed to analyze variants in 32 complement system genes for positive association with CAD. Functional variants in genes showing positive association were then identified by searching existing expression quantitative loci databases and validated by real-time reverse transcription polymerase chain reaction. A follow-up case-control design was then used to determine whether the functional variants are associated with CAD in 2 independent GeneID Chinese populations. Candidate pathway-based GWAS identified positive association between variants in C3AR1 and C6 and CAD. Two functional variants, rs7842 in C3AR1 and rs4400166 in C6, were found to be associated with expression levels of C3AR1 and C6, respectively. Significant association was identified between rs7842 and CAD (P=3.99×10(-6); odds ratio, 1.47) and between rs4400166 and CAD (P=9.30×10(-3); odds ratio, 1.24) in the validation cohort. The significant findings were confirmed in the replication cohort (P=1.53×10(-5); odds ratio, 1.37 for rs7842; P=8.41×10(-3); odds ratio, 1.21 for rs4400166).

CONCLUSIONS

Integration of GWAS with biological pathways and expression quantitative loci is effective in identifying new risk variants for CAD. Functional variants increasing C3AR1 and C6 expression were shown to confer significant risk of CAD for the first time.

Integrated genome-wide genotyping and gene expression profiling reveals BCL11B as a putative oncogene in acute myeloid leukemia with 14q32 aberrations

Acute myeloid leukemia is a neoplasm characterized by recurrent molecular aberrations traditionally demonstrated by cytogenetic analyses. We used high density genome-wide genotyping and gene expression profiling to reveal acquired cryptic abnormalities in acute myeloid leukemia. By genome-wide genotyping of 137 cases of primary acute myeloid leukemia, we disclosed a recurrent focal amplification on chromosome 14q32, which included the genes BCL11B, CCNK, C14orf177 and SETD3, in two cases. In the affected cases, the BCL11B gene showed consistently high mRNA expression, whereas the expression of the other genes was unperturbed. Fluorescence in situ hybridization on 40 cases of acute myeloid leukemia with high BCL11B mRNA expression [2.5-fold above median; 40 out of 530 cases (7.5%)] revealed 14q32 abnormalities in two additional cases. In the four BCL11B-rearranged cases the 14q32 locus was fused to different partner chromosomes. In fact, in two cases, we demonstrated that the focal 14q32 amplifications were integrated into transcriptionally active loci. The translocations involving BCL11B result in increased expression of full-length BCL11B protein. The BCL11B-rearranged acute myeloid leukemias expressed both myeloid and T-cell markers. These biphenotypic acute leukemias all carried FLT3 internal tandem duplications, a characteristic marker of acute myeloid leukemia. BCL11B mRNA expression in acute myeloid leukemia appeared to be strongly associated with expression of other T-cell-specific genes. Myeloid 32D(GCSF-R) cells ectopically expressing Bcl11b showed decreased proliferation rate and less maturation. In conclusion, by an integrated approach involving high-throughput genome-wide genotyping and gene expression profiling we identified BCL11B as a candidate oncogene in acute myeloid leukemia.

Genome-wide gene expression profiling, genotyping, and copy number analyses of acute myeloid leukemia using Affymetrix GeneChips

With novel genome-wide technologies it is nowadays possible to perform detailed molecular analyses of normal and malignant tissues. Acute myeloid leukemia (AML) is a heterogeneous group of diseases with variable response to therapy. Gene expression profiling and genome-wide genotyping have recently been successfully applied to unravel the heterogeneity of AML. This chapter gives instructions and recommendations for genome-wide gene expression analyses, genotyping, and copy number analyses, as performed for AML using Affymetrix GeneChips.

Promoter polymorphisms in two overlapping 6p25 genes implicate mitochondrial proteins in cognitive deficit in schizophrenia

In a previous study, we detected a 6p25-p24 region linked to schizophrenia in families with high composite cognitive deficit (CD) scores, a quantitative trait integrating multiple cognitive measures. Association mapping of a 10 Mb interval identified a 260 kb region with a cluster of single-nucleotide polymorphisms (SNPs) significantly associated with CD scores and memory performance. The region contains two colocalising genes, LYRM4 and FARS2, both encoding mitochondrial proteins. The two tagging SNPs with strongest evidence of association were located around the overlapping putative promoters, with rs2224391 predicted to alter a transcription factor binding site (TFBS). Sequencing the promoter region identified 22 SNPs, many predicted to affect TFBSs, in a tight linkage disequilibrium block. Luciferase reporter assays confirmed promoter activity in the predicted promoter region, and demonstrated marked downregulation of expression in the LYRM4 direction under the haplotype comprising the minor alleles of promoter SNPs, which however is not driven by rs2224391. Experimental evidence from LYRM4 expression in lymphoblasts, gel-shift assays and modelling of DNA breathing dynamics pointed to two adjacent promoter SNPs, rs7752203-rs4141761, as the functional variants affecting expression. Their C-G alleles were associated with higher transcriptional activity and preferential binding of nuclear proteins, whereas the G-A combination had opposite effects and was associated with poor memory and high CD scores. LYRM4 is a eukaryote-specific component of the mitochondrial biogenesis of Fe-S clusters, essential cofactors in multiple processes, including oxidative phosphorylation. LYRM4 downregulation may be one of the mechanisms involved in inefficient oxidative phosphorylation and oxidative stress, increasingly recognised as contributors to schizophrenia pathogenesis.

XX ovarian dysgenesis is caused by a PSMC3IP/HOP2 mutation that abolishes coactivation of estrogen-driven transcription

XX female gonadal dysgenesis (XX-GD) is a rare, genetically heterogeneous disorder characterized by lack of spontaneous pubertal development, primary amenorrhea, uterine hypoplasia, and hypergonadotropic hypogonadism as a result of streak gonads. Most cases are unexplained but thought to be autosomal recessive. We elucidated the genetic basis of XX-GD in a highly consanguineous Palestinian family by using homozygosity mapping and candidate-gene and whole-exome sequencing. Affected females were homozygous for a 3 bp deletion (NM_016556.2, c.600_602del) in the PSMC3IP gene, leading to deletion of a glutamic acid residue (p.Glu201del) in the highly conserved C-terminal acidic domain. Proteasome 26S subunit, ATPase, 3-Interacting Protein (PSMC3IP)/Tat Binding Protein Interacting Protein (TBPIP) is a nuclear, tissue-specific protein with multiple functions. It is critical for meiotic recombination as indicated by the known role of its yeast ortholog, Hop2. Through the C terminus (not present in yeast), PSMC3IP also coactivates ligand-driven transcription mediated by estrogen, androgen, glucocorticoid, progesterone, and thyroid nuclear receptors. In cell lines, the p.Glu201del mutation abolished PSMC3IP activation of estrogen-driven transcription. Impaired estrogenic signaling can lead to ovarian dysgenesis both by affecting the size of the follicular pool created during fetal development and by failing to counteract follicular atresia during puberty. PSMC3IP joins previous genes known to be mutated in XX-GD, the FSH receptor, and BMP15, highlighting the importance of hormonal signaling in ovarian development and maintenance and suggesting a common pathway perturbed in isolated XX-GD. By analogy to other XX-GD genes, PSMC3IP is also a candidate gene for premature ovarian failure, and its role in folliculogenesis should be further investigated.

NUP98/NSD1 characterizes a novel poor prognostic group in acute myeloid leukemia with a distinct HOX gene expression pattern

Translocations involving nucleoporin 98kD (NUP98) on chromosome 11p15 occur at relatively low frequency in acute myeloid leukemia (AML) but can be missed with routine karyotyping. In this study, high-resolution genome-wide copy number analyses revealed cryptic NUP98/NSD1 translocations in 3 of 92 cytogenetically normal (CN)-AML cases. To determine their exact frequency, we screened > 1000 well-characterized pediatric and adult AML cases using a NUP98/NSD1-specific RT-PCR. Twenty-three cases harbored the NUP98/NSD1 fusion, representing 16.1% of pediatric and 2.3% of adult CN-AML patients. NUP98/NSD1-positive AML cases had significantly higher white blood cell counts (median, 147 × 10⁹/L), more frequent FAB-M4/M5 morphology (in 63%), and more CN-AML (in 78%), FLT3/internal tandem duplication (in 91%) and WT1 mutations (in 45%) than NUP98/NSD1-negative cases. NUP98/NSD1 was mutually exclusive with all recurrent type-II aberrations. Importantly, NUP98/NSD1 was an independent predictor for poor prognosis; 4-year event-free survival was < 10% for both pediatric and adult NUP98/NSD1-positive AML patients. NUP98/NSD1-positive AML showed a characteristic HOX-gene expression pattern, distinct from, for example, MLL-rearranged AML, and the fusion protein was aberrantly localized in nuclear aggregates, providing insight into the leukemogenic pathways of these AMLs. Taken together, NUP98/NSD1 identifies a previously unrecognized group of young AML patients, with distinct characteristics and dismal prognosis, for whom new treatment strategies are urgently needed.

Rheumatoid arthritis-associated polymorphisms are not protective against Alzheimer's disease

Background

Rheumatoid arthritis (RA) and Alzheimer's disease (AD) are inversely associated. To test the hypothesis that genetic elements associated with increased RA risk are associated with decreased AD risk, we evaluated RA genetic risk factors recently identified in genome-wide association studies (GWAS) for their association with AD in a two-stage, case-control analysis.

Results

In our Stage 1 analysis of ~800 AD and ~1,200 non-AD individuals, three of seventeen RA-associated SNPs were nominally associated with AD (p < 0.05) with one SNP, rs2837960, retaining significance after correction for multiple testing (p = 0.03). The rs2837960_G (minor) allele, which is associated with increased RA risk, was associated with increased AD risk. Analysis of these three SNPs in a Stage 2 population, consisting of ~1,100 AD and ~2,600 non-AD individuals, did not confirm their association with AD. Analysis of Stage 1 and 2 combined suggested that rs2837960 shows a trend for association with AD. When the Stage 2 population was age-matched for the Stage 1 population, rs2837960 exhibited a non-significant trend with AD. Combined analysis of Stage 1 and the age-matched Stage 2 subset showed a significant association of rs2837960 with AD (p = 0.002, OR 1.24) that retained significance following correction for age, sex and APOE (p = 0.02, OR = 1.20). Rs2837960 is near BACE2, which encodes an aspartic protease capable of processing the AD-associated amyloid precursor protein. Testing for an association between rs2837960 and the expression of BACE2 isoforms in human brain, we observed a trend between rs2837960 and the total expression of BACE2 and the expression of a BACE2 transcript lacking exon 7 (p = 0.07 and 0.10, respectively).

Conclusions

RA-associated SNPs are generally not associated with AD. Moreover, rs2837960_G is associated with increased risk of both RA and, in individuals less than 80 years of age, with AD. Overall, these results contest the hypothesis that genetic variants associated with RA confer protection against AD. Further investigation of rs2837960 is necessary to elucidate the mechanism by which rs2837960 contributes to both AD and RA risk, likely via modulation of BACE2 expression.

SIGMA2: a system for the integrative genomic multi-dimensional analysis of cancer genomes, epigenomes, and transcriptomes

Background

High throughput microarray technologies have afforded the investigation of genomes, epigenomes, and transcriptomes at unprecedented resolution. However, software packages to handle, analyze, and visualize data from these multiple 'omics disciplines have not been adequately developed.

Results

Here, we present SIGMA², a system for the integrative genomic multi-dimensional analysis of cancer genomes, epigenomes, and transcriptomes. Multi-dimensional datasets can be simultaneously visualized and analyzed with respect to each dimension, allowing combinatorial integration of the different assays belonging to the different 'omics.

Conclusion

The identification of genes altered at multiple levels such as copy number, loss of heterozygosity (LOH), DNA methylation and the detection of consequential changes in gene expression can be concertedly performed, establishing SIGMA² as a novel tool to facilitate the high throughput systems biology analysis of cancer.