Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants.

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n=197, temporal cortex n=202) and with other brain pathologies (non-AD, cerebellar n=177, temporal cortex n=197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ± 100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non-ADs (q<0.05, p=7.70 × 10(-5)-1.67 × 10(-82)). Of these, 2,089 were also significant in the temporal cortex (p=1.85 × 10(-5)-1.70 × 10(-141)). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10(-6)). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9-3.3 fold enrichment (p<10(-6)) of significant cisSNPs with suggestive AD-risk association (p<10(-3)) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non-CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

Quantitative miRNA expression analysis using fluidigm microfluidics dynamic arrays.

BACKGROUND: MicroRNAs (miRNAs) represent a growing class of small non-coding RNAs that are important regulators of gene expression in both plants and animals. Studies have shown that miRNAs play a critical role in human cancer and they can influence the level of cell proliferation and apoptosis by modulating gene expression. Currently, methods for the detection and measurement of miRNA expression include small and moderate-throughput technologies, such as standard quantitative PCR and microarray based analysis. However, these methods have several limitations when used in large clinical studies where a high-throughput and highly quantitative technology needed for the efficient characterization of a large number of miRNA transcripts in clinical samples. Furthermore, archival formalin fixed, paraffin embedded (FFPE) samples are increasingly becoming the primary resource for gene expression studies because fresh frozen (FF) samples are often difficult to obtain and requires special storage conditions. In this study, we evaluated the miRNA expression levels in FFPE and FF samples as well as several lung cancer cell lines employing a high throughput qPCR-based microfluidic technology. The results were compared to standard qPCR and hybridization-based microarray platforms using the same samples. RESULTS: We demonstrated highly correlated Ct values between multiplex and singleplex RT reactions in standard qPCR assays for miRNA expression using total RNA from A549 (R = 0.98; p < 0.0001) and H1299 (R = 0.95; p < 0.0001) lung cancer cell lines. The Ct values generated by the microfluidic technology (Fluidigm 48.48 dynamic array systems) resulted in a left-shift toward lower Ct values compared to those observed by ABI 7900 HT (mean difference, 3.79), suggesting that the microfluidic technology exhibited a greater sensitivity. In addition, we show that as little as 10 ng total RNA can be used to reliably detect all 48 or 96 tested miRNAs using a 96-multiplexing RT reaction in both FFPE and FF samples. Finally, we compared miRNA expression measurements in both FFPE and FF samples by qPCR using the 96.96 dynamic array and Affymetrix microarrays. Fold change comparisons for comparable genes between the two platforms indicated that the overall correlation was R = 0.60. The maximum fold change detected by the Affymetrix microarray was 3.5 compared to 13 by the 96.96 dynamic array. CONCLUSION: The qPCR-array based microfluidic dynamic array platform can be used in conjunction with multiplexed RT reactions for miRNA gene expression profiling. We showed that this approach is highly reproducible and the results correlate closely with the existing singleplex qPCR platform at a throughput that is 5 to 20 times higher and a sample and reagent usage that was approximately 50-100 times lower than conventional assays. We established optimal conditions for using the Fluidigm microfluidic technology for rapid, cost effective, and customizable arrays for miRNA expression profiling and validation.

DNA methylation and RNA expression profiles in lung adenocarcinomas of never-smokers.

Lung cancer occurs in never-smokers. Epigenetic changes in lung cancer potentially represent important diagnostic, prognostic, and therapeutic targets. We compared DNA methylation profiles of 28 adenocarcinomas of the lungs of never-smokers with paired adjacent nonmalignant lung tissue. We correlated differential methylation changes with gene expression changes from the same 28 sample pairs. Using principal component analysis, we observed a distinct separation in methylation profiles between tumor and adjacent nonmalignant lung tissue. Tumors were generally hypomethylated compared with adjacent nonmalignant tissue. Of 1,906 CpG sites differentially methylated between tumor and nonmalignant tissue, 1,198 were within classically defined CpG islands where tumors were hypermethylated compared with nonmalignant tissue. A total of 708 sites were outside CpG islands where tumors were hypomethylated compared with nonmalignant tissue. There were significant differences in expression of 351 genes (23%) of the 1,522 genes matched to the differentially methylated CpG sites. Genes that were not significantly differentially expressed and were hypermethylated within CpG sites were enriched for homeobox genes. These results suggest that the methylation profiles of lung adenocarcinomas of never-smokers and adjacent nonmalignant lung tissue are significantly different. Despite the differential methylation of homeobox genes, no significant changes in expression of these genes were detected.

Methylmercury induces apoptosis in cultured rat dorsal root ganglion neurons.

Methylmercury is known to have devastating effects on the mammalian nervous system. In order to characterize the dose dependence of methylmercury-induced neurotoxicity, we first studied neurite outgrowth from rat dorsal root ganglia explants. In this model, methylmercury inhibited neurite outgrowth with a TD(50) of approximately 0.5 microM. We then used this relationship to optimize dosing for subsequent transcriptional profiling analyses in two independent neuronal model systems: dissociated sensory neurons and PC12 cells. As seen in previous studies, the expression of a number of genes associated with oxidative stress was altered following a 6h challenge with 1 microM methylmercury. When PC12 cells were subjected to a longer exposure (24h), a relative increase was noted in the representation of genes associated with cell cycling and apoptosis. To confirm the presence of apoptosis in cultured neurons, we then applied TUNEL staining and bis-benzimide staining techniques to primary cultures of dissociated sensory neurons. After 24h, 1 microM methylmercury increased both DNA end-labeling (P<0.01) and nuclear fragmentation (P<0.02). The latter effect appeared to be dose-dependent.

Association of MAPT haplotypes with Alzheimer's disease risk and MAPT brain gene expression levels.

INTRODUCTION: MAPT encodes for tau, the predominant component of neurofibrillary tangles that are neuropathological hallmarks of Alzheimer's disease (AD). Genetic association of MAPT variants with late-onset AD (LOAD) risk has been inconsistent, although insufficient power and incomplete assessment of MAPT haplotypes may account for this. METHODS: We examined the association of MAPT haplotypes with LOAD risk in more than 20,000 subjects (n-cases = 9,814, n-controls = 11,550) from Mayo Clinic (n-cases = 2,052, n-controls = 3,406) and the Alzheimer's Disease Genetics Consortium (ADGC, n-cases = 7,762, n-controls = 8,144). We also assessed associations with brain MAPT gene expression levels measured in the cerebellum (n = 197) and temporal cortex (n = 202) of LOAD subjects. Six single nucleotide polymorphisms (SNPs) which tag MAPT haplotypes with frequencies greater than 1% were evaluated. RESULTS: H2-haplotype tagging rs8070723-G allele associated with reduced risk of LOAD (odds ratio, OR = 0.90, 95% confidence interval, CI = 0.85-0.95, p = 5.2E-05) with consistent results in the Mayo (OR = 0.81, p = 7.0E-04) and ADGC (OR = 0.89, p = 1.26E-04) cohorts. rs3785883-A allele was also nominally significantly associated with LOAD risk (OR = 1.06, 95% CI = 1.01-1.13, p = 0.034). Haplotype analysis revealed significant global association with LOAD risk in the combined cohort (p = 0.033), with significant association of the H2 haplotype with reduced risk of LOAD as expected (p = 1.53E-04) and suggestive association with additional haplotypes. MAPT SNPs and haplotypes also associated with brain MAPT levels in the cerebellum and temporal cortex of AD subjects with the strongest associations observed for the H2 haplotype and reduced brain MAPT levels (ß = -0.16 to -0.20, p = 1.0E-03 to 3.0E-03). CONCLUSIONS: These results confirm the previously reported MAPT H2 associations with LOAD risk in two large series, that this haplotype has the strongest effect on brain MAPT expression amongst those tested and identify additional haplotypes with suggestive associations, which require replication in independent series. These biologically congruent results provide compelling evidence to screen the MAPT region for regulatory variants which confer LOAD risk by influencing its brain gene expression.

Multi-platform analysis of microRNA expression measurements in RNA from fresh frozen and FFPE tissues.

MicroRNAs play a role in regulating diverse biological processes and have considerable utility as molecular markers for diagnosis and monitoring of human disease. Several technologies are available commercially for measuring microRNA expression. However, cross-platform comparisons do not necessarily correlate well, making it difficult to determine which platform most closely represents the true microRNA expression level in a tissue. To address this issue, we have analyzed RNA derived from cell lines, as well as fresh frozen and formalin-fixed paraffin embedded tissues, using Affymetrix, Agilent, and Illumina microRNA arrays, NanoString counting, and Illumina Next Generation Sequencing. We compared the performance within- and between the different platforms, and then verified these results with those of quantitative PCR data. Our results demonstrate that the within-platform reproducibility for each method is consistently high and although the gene expression profiles from each platform show unique traits, comparison of genes that were commonly detectable showed that detection of microRNA transcripts was similar across multiple platforms.

Effects of joint contracture on the contralateral unoperated limb in a rabbit knee contracture model: a biomechanical and genetic study.

In most animal models, unoperated contralateral limbs are used as controls. However, in some experimental circumstances, the contralateral limb may represent a skewed control. The main purpose of this study was to determine if the unoperated contralateral limb could be used as a control, or if a different unoperated animal's limb should be used instead. Seventeen rabbits were divided into two groups. Group 1 rabbits (n = 12) underwent surgery on their right limbs to induce a contracture. Group 2 rabbits (n = 5) underwent no surgery. The left non-operated limbs of rabbits in group 1 were biomechanically and genetically compared to the limbs of unoperated rabbits in group 2 with the use of a validated joint measuring device and custom microarray, respectively. After 8 weeks of immobilization, there was a statistically greater flexion contracture in the unoperated contralateral limbs compared to the limbs of animals that received no surgery(8.4 ± 8.9° vs. 0 ± 0°; p-value = 0.03). When animals were remobilized for an additional 16 weeks, the significance between groups was lost (11.9 ± 21.4° vs. 8.9 ± 9.5°; p = 0.38). Similarly, there was a statistically significant increase in nine genes at 8 weeks (p < 0.001). However, at 24 weeks, only the PMCA 1 gene was statically increased (p < 0.001). In our rabbit model, the non-operated limb develops a small flexion contracture at 8 weeks. After 16 weeks of remobilization, there is no biomechanical or genetic difference between contralateral non-operated limbs and limbs of animals not undergoing any surgical intervention. Given the biomechanical and genetic findings, the contralateral non-operated limb can be used as a valid control.

Increased miR-708 expression in NSCLC and its association with poor survival in lung adenocarcinoma from never smokers.

PURPOSE: miRNA plays an important role in human disease and cancer. We seek to investigate the expression status, clinical relevance, and functional role of miRNA in non-small cell lung cancer. EXPERIMENTAL DESIGN: We conducted miRNA expression profiling in matched lung adenocarcinoma and uninvolved lung using 56 pairs of fresh-frozen (FF) and 47 pairs of formalin-fixed, paraffin-embedded (FFPE) samples from never smokers. The most differentially expressed miRNA genes were evaluated by Cox analysis and log-rank test. Among the best candidate, miR-708 was further examined for differential expression in two independent cohorts. Functional significance of miR-708 expression in lung cancer was examined by identifying its candidate mRNA target and through manipulating its expression levels in cultured cells. RESULTS: Among the 20 miRNAs most differentially expressed between tested tumor and normal samples, high expression level of miR-708 in the tumors was most strongly associated with an increased risk of death after adjustments for all clinically significant factors including age, sex, and tumor stage (FF cohort: HR, 1.90; 95% CI, 1.08-3.35; P = 0.025 and FFPE cohort: HR, 1.93; 95% CI, 1.02-3.63; P = 0.042). The transcript for TMEM88 gene has a miR-708 binding site in its 3' UTR and was significantly reduced in tumors high of miR-708. Forced miR-708 expression reduced TMEM88 transcript levels and increased the rate of cell proliferation, invasion, and migration in culture. CONCLUSIONS: miRNA-708 acts as an oncogene contributing to tumor growth and disease progression by directly downregulating TMEM88, a negative regulator of the Wnt signaling pathway in lung cancer.

Glutathione S-transferase omega genes in Alzheimer and Parkinson disease risk, age-at-diagnosis and brain gene expression: an association study with mechanistic implications.

BACKGROUND: Glutathione S-transferase omega-1 and 2 genes (GSTO1, GSTO2), residing within an Alzheimer and Parkinson disease (AD and PD) linkage region, have diverse functions including mitigation of oxidative stress and may underlie the pathophysiology of both diseases. GSTO polymorphisms were previously reported to associate with risk and age-at-onset of these diseases, although inconsistent follow-up study designs make interpretation of results difficult. We assessed two previously reported SNPs, GSTO1 rs4925 and GSTO2 rs156697, in AD (3,493 ADs vs. 4,617 controls) and PD (678 PDs vs. 712 controls) for association with disease risk (case-controls), age-at-diagnosis (cases) and brain gene expression levels (autopsied subjects). RESULTS: We found that rs156697 minor allele associates with significantly increased risk (odds ratio = 1.14, p = 0.038) in the older ADs with age-at-diagnosis > 80 years. The minor allele of GSTO1 rs4925 associates with decreased risk in familial PD (odds ratio = 0.78, p = 0.034). There was no other association with disease risk or age-at-diagnosis. The minor alleles of both GSTO SNPs associate with lower brain levels of GSTO2 (p = 4.7 × 10-11-1.9 × 10-27), but not GSTO1. Pathway analysis of significant genes in our brain expression GWAS, identified significant enrichment for glutathione metabolism genes (p = 0.003). CONCLUSION: These results suggest that GSTO locus variants may lower brain GSTO2 levels and consequently confer AD risk in older age. Other glutathione metabolism genes should be assessed for their effects on AD and other chronic, neurologic diseases.

Novel late-onset Alzheimer disease loci variants associate with brain gene expression.

OBJECTIVE: Recent genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) identified 9 novel risk loci. Discovery of functional variants within genes at these loci is required to confirm their role in Alzheimer disease (AD). Single nucleotide polymorphisms that influence gene expression (eSNPs) constitute an important class of functional variants. We therefore investigated the influence of the novel LOAD risk loci on human brain gene expression. METHODS: We measured gene expression levels in the cerebellum and temporal cortex of autopsied AD subjects and those with other brain pathologies (∼400 total subjects). To determine whether any of the novel LOAD risk variants are eSNPs, we tested their cis-association with expression of 6 nearby LOAD candidate genes detectable in human brain (ABCA7, BIN1, CLU, MS4A4A, MS4A6A, PICALM) and an additional 13 genes ±100 kb of these SNPs. To identify additional eSNPs that influence brain gene expression levels of the novel candidate LOAD genes, we identified SNPs ±100 kb of their location and tested for cis-associations. RESULTS: CLU rs11136000 (p = 7.81 × 10(-4)) and MS4A4A rs2304933/rs2304935 (p = 1.48 × 10(-4)-1.86 × 10(-4)) significantly influence temporal cortex expression levels of these genes. The LOAD-protective CLU and risky MS4A4A locus alleles associate with higher brain levels of these genes. There are other cis-variants that significantly influence brain expression of CLU and ABCA7 (p = 4.01 × 10(-5)-9.09 × 10(-9)), some of which also associate with AD risk (p = 2.64 × 10(-2)-6.25 × 10(-5)). CONCLUSIONS: CLU and MS4A4A eSNPs may at least partly explain the LOAD risk association at these loci. CLU and ABCA7 may harbor additional strong eSNPs. These results have implications in the search for functional variants at the novel LOAD risk loci.

Endothelial nitric oxide synthase gene variation associated with chronic kidney disease after liver transplant.

OBJECTIVE: To identify single nucleotide polymorphisms (SNPs) associated with risk of developing chronic kidney disease (CKD), a prevalent comorbidity, after liver transplant (LT). PATIENTS AND METHODS: This study consists of a cohort of adult (> or =18 years) primary-LT recipients who had normal renal function before LT and who survived 1 year or more after LT at a high-volume US LT program between January 1, 1990, and December 31, 2000. Patients with adequate renal function (estimated glomerular filtration rate, > or =40 mL/min per 1.73 m(2) during follow-up; n=308) and patients with incident CKD (estimated glomerular filtration rate, <40 mL/min per 1.73 m(2) after LT; n=92) were identified. To investigate the association of 6 candidate genes with post-LT CKD, we selected SNPs that have been associated with renal function in the literature. Hazard ratios were estimated using Cox regression, adjusted for potential confounding variables. RESULTS: The variant allele (298Asp) of the Glu298Asp SNP in the endothelial nitric oxide synthase gene (NOS3) was significantly associated with CKD after LT (P=.05; adjusted for multiple comparisons). The 5-year incidence of CKD was 70% among patients homozygous for the NOS3 variant allele (298Asp) compared with 42% among those not homozygous for the NOS3 variant allele. Specifically, homozygosity for the NOS3 variant allele conferred a 2.5-fold increased risk of developing CKD after LT (P=.005, adjusted for confounding variables). CONCLUSION: Homozygosity for the variant allele of NOS3 (298Asp) is associated with CKD after LT and may be useful for identifying recipients at higher risk of post-LT CKD.

Casp8p41 expression in primary T cells induces a proinflammatory response.

OBJECTIVE: HIV infection of CD4 T cells can lead to HIV protease-mediated cleavage of procaspase 8 generating a novel, HIV-specific peptide called Casp8p41. Casp8p41 has at least two biologic functions: induction of cell death via mitochondrial depolarization and release of cytochrome C, as well as activation of nuclear factor kappa B (NFkappaB). We have previously shown that Casp8p41-induced NFkappaB activation enhances HIV LTR transcription and consequently increases HIV replication. Herein, we questioned whether Casp8p41-induced NFkappaB activation impacts the cytokine profile of cells expressing Casp8p41. DESIGN: Analysis of cells expressing Casp8p41 and HIV-infected T cells. METHODS: We assessed whether host genes are transcriptionally activated following Casp8p41 production, using microarray analysis, cytokine quantification, followed by western blot and flow cytometry. RESULTS: Microarray analysis identified 259 genes significantly upregulated following expression of Casp8p41. Furthermore, Casp8p41 expression in primary CD4 T cells results in increased production of interleukin (IL)-2, IL-15 and tumor necrosis factor (TNF), as well as IL-1RA; whereas levels of granulocyte macrophage colony-stimulating factor and interferon (IFN)-gamma were reduced in the Casp8p41 expressing cells. Intracellular flow cytometry confirmed the co-association of Casp8p41 with elevated TNF in HIV-infected cells. CONCLUSION: These data indicate that the expression of Casp8p41 in HIV-infected CD4 T cells in addition to promoting apoptosis and enhancing HIV replication also promotes a proinflammatory cytokine milieu, which is characteristic of untreated HIV infection.